Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,601 C→A 7.7% intergenic (+42/+274) ybaQ → / ← copA putative DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,6010CA 7.7% 169.7 / 4.8 58intergenic (+42/+274)ybaQ/copAputative DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base A (2/3);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CAGAGTTGCTGCATGTTCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTT  >  NC_000913/508309‑508888
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                               
cAGAGTTGCTGCATGTTCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCgca                                                                                                                                                                                                                                                                                         >  2:73443/1‑301 (MQ=255)
       gttggaTGTTCATCGTAATTGCGTCAGTGTACTGATCAATAACATTCGTTAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATGCGGCGTGAACGCCttatccggcc                                                                                                                                                                                                                                                                                  <  1:609677/296‑2 (MQ=255)
         tgcATGTTCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCACTTCATCTCGTAGACCGGAGAAGGCGTTACCGCCcaatccggaatc                                                                                                                                                                                                                                                                                >  1:120926/1‑301 (MQ=255)
              tttCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGGATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAga                                                                                                                                                                                                                                                                           <  1:20096/300‑2 (MQ=255)
                    acccgacgctcttccgatcTGATTAATAACAATCGTAAACTCACTTCTGAGTTGGCATTTCGTCTGGCGAAAGTTTTTGATACCCCAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAgatgaagc                                                                                                                                                                                                                                                                     <  1:432828/283‑8 (MQ=255)
                     gTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTg                                                                                                                                                                                                                                                                    >  2:201296/1‑301 (MQ=255)
                           gCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCAACCGGCAACCGAGCACCATTGCCTGAt                                                                                                                                                                                                                                                              >  1:461858/1‑301 (MQ=255)
                                    cACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGAAGCGACGCtt                                                                                                                                                                                                                                                     >  2:183824/1‑301 (MQ=255)
                                        gATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAgatgaagaagatcggaagagcacacgtc                                                                                                                                                                                                                                                 >  2:433392/1‑274 (MQ=255)
                                         aTCAATAACAATCGTAAACTCACTATTGGGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCg                                                                                                                                                                                                                                                <  1:73353/301‑1 (MQ=255)
                                            aaTAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCGTCATCTCGTAGGCCGGATAAGGCGTTCACACAGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCt                                                                                                                                                                                                                                             >  2:385249/1‑301 (MQ=255)
                                                 cAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGCTCTTTGGGAAGTTGCAAACAACATGCGCACCCAGGAAGACTTGGGACGGCTTGCAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCACAAAAGGTCGCGTAAGGCCAACACATCCCTTCATCTCGTCAGCCGGACAAGGCTTTCACGCCCCATCCAGCCTCCGCCAACCATTGCCCGAtcagcaccctacccgccgcaacc                                                                                                                                                                                                                                        >  1:373367/1‑279 (MQ=255)
                                                 cAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCa                                                                                                                                                                                                                                        >  2:442683/1‑301 (MQ=255)
                                                             cacttctttgaTGGCATTTCGTCTGGCGAAAGTTGTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                                                                                                                                                                                            <  1:201053/292‑1 (MQ=255)
                                                                tactGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                                                                                                                                                                                         >  2:639625/1‑301 (MQ=255)
                                                                 tctGAGATTGCATTTCGTTTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgtg                                                                                                                                                                                                                        <  2:546262/300‑1 (MQ=255)
                                                                   ggagaTGGCATTTCGTCTTGCGAAAGTTTTGTATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGcc                                                                                                                                                                                                                      <  1:488106/300‑1 (MQ=255)
                                                                           cATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCCGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCCCGCCGTGAAGAGCGTGCAACAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTCGGCCGGATAAGGCGTTCACGCCGCATCCGTCATCCGAGCCCCCTTGCCTGATGCGACGCTTACCCGTCTTATCCGGCCTACCAACTTATTCCTGAcccct                                                                                                                                                                                                              >  1:152706/1‑299 (MQ=255)
                                                                                  tCTGTCGAAAGTTTTTGATACCACTGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCgg                                                                                                                                                                                                       <  2:82186/301‑1 (MQ=255)
                                                                                    tGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGGTCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTATGCCTGAACCGTAGGCCAGaa                                                                                                                                                                                                     >  2:166242/1‑300 (MQ=255)
                                                                                      gCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                                                                                                                                                                                   >  2:56411/1‑301 (MQ=255)
                                                                                              gaggtttgtccgttctgATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAAGATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTa                                                                                                                                                                                           <  1:497714/285‑1 (MQ=255)
                                                                                                ttGATACCTCAGTCGATTTTTGGCTAACCCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                                                                         <  2:462422/301‑1 (MQ=255)
                                                                                                  gATACCACAGTCGATGTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGAGTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                                                       <  2:152857/301‑1 (MQ=255)
                                                                                                      ccACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcaa                                                                                                                                                                                   >  2:265907/1‑300 (MQ=255)
                                                                                                       ccctgTCGATTTTTGGCTAAACCTCCGTGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcacc                                                                                                                                                                                  <  1:183583/297‑3 (MQ=255)
                                                                                                       aacaGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatc                                                                                                                                                                                  <  1:56343/300‑1 (MQ=255)
                                                                                                            tCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggca                                                                                                                                                                             >  2:290180/1‑301 (MQ=255)
                                                                                                              gATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCCGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTaagatcggaagagcaca                                                                                                                                                                           >  2:498361/1‑285 (MQ=255)
                                                                                                               aTTTTTGGGTAAACGTCCAGGGGGCGGTTGATCTTTGGGAAGTTGAAAACAACAGGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggcatcc                                                                                                                                                                          <  1:289778/301‑1 (MQ=255)
                                                                                                                        tAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACACACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTCCGCCGCATCCGGCCTCCCAGcccaat                                                                                                                                                                 >  1:447633/1‑296 (MQ=255)
                                                                                                                             cTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAAGCCGCATCCGGCATCCGAGAAACATTAccc                                                                                                                                                            >  2:331451/1‑300 (MQ=255)
                                                                                                                                    gggcggTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGGGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGa                                                                                                                                                     <  1:442119/300‑1 (MQ=255)
                                                                                                                                               tcttGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGt                                                                                                                                          <  2:77954/299‑1 (MQ=255)
                                                                                                                                                  tGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCATCCGAGCCACATTGCCCGATACGACGCTTACGCGtaga                                                                                                                                       >  1:77864/1‑298 (MQ=255)
                                                                                                                                                   gggAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa                                                                                                                                      >  2:527268/1‑301 (MQ=255)
                                                                                                                                                          tGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCCGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAACAGGCGTTTACGCCGCATACGGCATCCGAGCACCAATGCCTGATCCGACGCTTAACCCTCTTAACAGccc                                                                                                                               >  1:49334/1‑301 (MQ=37)
                                                                                                                                                                          ctccAGGAATAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTg                                                                                                               <  2:121050/299‑1 (MQ=255)
                                                                                                                                                                              aGGAAGTATTGGGACGGTGTTAATCTGTGGCGGAATGTTTGGCACGCTGTGAAGATCGGGCAAAAAAGGTCGCGTAAGGCTAAAACATCGCTTCATCGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACGTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAAcc                                                                                                           <  1:166046/301‑1 (MQ=25)
                                                                                                                                                                                      ttGGGACGTATTGAAACAGTGGCTGATTATTTGGCACGCCGTGAAGAGCGTGAAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCg                                                                                                   <  2:49402/301‑1 (MQ=255)
                                                                                                                                                                                                        gTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgccg                                                                                 <  2:448197/301‑1 (MQ=255)
                                                                                                                                                                                                             tGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTATCCCGGATAAACCGACTACGCCCCATcc                                                                            >  1:110241/1‑301 (MQ=255)
                                                                                                                                                                                                                  atTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCACCCCGCATCCGAGCACCATTGCCTGATGCGACGCTTCAGCCTCCTATCACGCCTACAAACCTGTGCCTGAACCCCCAGCAGGAAAAAGCGATAACGCAGCACCCGGCAg                                                                       >  1:353678/1‑301 (MQ=25)
                                                                                                                                                                                                                     ggggaCGCCGTTAAGGGCGGGCAAAAAAGGTCGCGTAAGGCATAAAAATCGCTTCTGCTCGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATCCGAGCCCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCa                                                                    <  1:301499/297‑1 (MQ=25)
                                                                                                                                                                                                                        cACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTTCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAACCCTGGGACTTAACCCTCAGCCGGaaaaagagatacacacgcataccggaaccacag                                                                 >  1:412733/1‑272 (MQ=25)
                                                                                                                                                                                                                                    agCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGGCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGATTACGCCGCATCCGGCaaccatccctcacaacacaa                                                     >  1:146402/1‑282 (MQ=255)
                                                                                                                                                                                                                                     gCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAGAAGGCGTTTACGCCGCATACGGCATCCGAGCACCATTGACTGATGCGACGCTTACGCGTCCTAACAGGCCCAACAACCTGTACCTGAAACGTAGACCGGATAAGGCGTTTACGCCCACACCGGCCGTAATGCGCCGATCCCacaa                                                    >  1:353133/1‑298 (MQ=25)
                                                                                                                                                                                                                                            aaaacGGTTGCGTAAGGCCAATACGTCGAGTTATCTCGTAGGCCGGTGAAGGCGTTCACGCCACTGCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTGATCAGGCCTACAAACGTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCAc                                             <  1:384753/301‑1 (MQ=25)
                                                                                                                                                                                                                                                 ggTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTaa                                        >  2:553779/1‑301 (MQ=255)
                                                                                                                                                                                                                                                  gTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTaaa                                       >  2:477452/1‑301 (MQ=255)
                                                                                                                                                                                                                                                     gcgTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCAACCCTAAaaca                                    >  2:583843/1‑301 (MQ=255)
                                                                                                                                                                                                                                                              cAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCg                           <  2:146528/301‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaatcgtcggttgatctgtgaGGCGGGATAAGGCGTGCAGTCCTCATCGGGCAGCCGGGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCAGCCGGCATCCGAGCACGATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGAGAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGc                          <  1:638748/281‑1 (MQ=37)
                                                                                                                                                                                                                                                                    aTCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAatga                     >  2:330891/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                             tctcGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCATCCGATCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTAt            <  2:353596/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGCCGCTTACGCGTCTTATCAGGCCTACCAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCCGCCATGCGTCGACGCCCAATGCCCCACCCTAAaccacccccacccccatcacgccccactgcc         >  1:243070/1‑274 (MQ=25)
                                                                                                                                                                                                                                                                                   gggCCGGATAAGGCGTTCATGCCGCATCCGGCATCCGAGGACCTGTGCCTGATGCGACTCTTAGCCGTCTTATCAGGCCTACAAACTTGGGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTc      <  1:476851/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       cGGATTAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCTTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGtt  <  2:413296/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                               
CAGAGTTGCTGCATGTTCATCGTAATAGCGTCAGTGCACTGATCAATAACAATCGTAAACTCACTACTGAGATGGCATTTCGTCTGGCGAAAGTTTTTGATACCACAGTCGATTTTTGGCTAAACCTCCAGGCGGCGGTTGATCTTTGGGAAGTTGAAAACAACATGCGCACCCAGGAAGAATTGGGACGGATTGAAACAGTGGCTGAATATTTGGCACGCCGTGAAGAGCGTGCAAAAAAGGTCGCGTAAGGCCAAAACATCGCTTCATCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACCTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTT  >  NC_000913/508309‑508888

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑