Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,331,850 T→C 9.4% V268A (GTG→GCG)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,331,8500TC 9.4% 119.1 / 3.5 47V268A (GTG→GCG) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base T (16/25);  new base C (2/3);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGC  >  NC_000913/1331559‑1332147
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                          
ttgttctcgTGTTGGGGTTTATTGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTTGTCTGTTTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGTTGGTCGGGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAATTCGATGCCAGCACGACCATGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGAAAAACCGTTACGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCACGCTGTAGCGCGCtttttcc                                                                                                                                                                                                                                                                                                   <  2:608326/294‑7 (MQ=255)
tGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCAACAGCGTGCCGGAACg                                                                                                                                                                                                                                                                                                   >  1:164457/1‑301 (MQ=255)
ggggccGCGTTGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGGCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCACGCTGTAGCGCGCtttttcc                                                                                                                                                                                                                                                                                                   <  1:378423/297‑7 (MQ=255)
      gcgTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCACGCTGTAGCGCGCtttttccagcaga                                                                                                                                                                                                                                                                                             <  2:346135/301‑13 (MQ=255)
            tttccgatctGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAcc                                                                                                                                                                                                                                                                                       <  2:8540/292‑1 (MQ=255)
            tggGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAACTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGACCGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTACGGTAAGTCTGCTGGAAAAAACACGCTACAGCGCGCTGGAACGGAAAGAAAAAcc                                                                                                                                                                                                                                                                                       >  1:124505/1‑301 (MQ=255)
            gggggTTTTTGGTTTGGCCGCTGCTTTGGAAAAAGATCGTTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAAcc                                                                                                                                                                                                                                                                                       <  1:69032/300‑1 (MQ=255)
               ggTATATGTTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGAc                                                                                                                                                                                                                                                                                    <  2:133783/301‑1 (MQ=255)
                      ggTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGCAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGATACAGCGTGCTGGAACGTGAAGACaaaacaaaacggaag                                                                                                                                                                                                                                                                             >  1:8540/1‑289 (MQ=255)
                          gtgcgggtggtATGGAAAAAGATCTGTCGTGGCCTGTCTGCCGGTCGTGTGGAGTCGGGGGCGGTTCGCCTGGTGGGCGAGCGTGAGCGTGAAATTAAAGCGTGCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAAc                                                                                                                                                                                                                                                                         <  1:279529/291‑1 (MQ=255)
                            gccgcTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCAGGAACGTGAAGACAAACCGACAACCAGaaaaaca                                                                                                                                                                                                                                                                       >  2:108936/1‑298 (MQ=255)
                                     aTGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGGAAGTCTGCTGGAAAAAGCGCGATCACGCCGACTGGAACGTGAAGACAAACCGACAACAAAAAAACCCggaccacac                                                                                                                                                                                                                                                              >  1:251185/1‑294 (MQ=255)
                                            aaGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTttatta                                                                                                                                                                                                                                                       <  2:142959/301‑1 (MQ=255)
                                              ggtggtttgtgtctgtgcTGCCGGTTGGTTGCAGGTGGGGGCGGTTCGCCTGGTGGTCGAGGGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTGCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTAcc                                                                                                                                                                                                                                                     <  1:184754/284‑1 (MQ=255)
                                                 cGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAAACTGGCAATCCTTTTATTACctaa                                                                                                                                                                                                                                                  >  2:268939/1‑299 (MQ=255)
                                                  cgctcgtccgttcttgctgtgtgtgtGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTa                                                                                                                                                                                                                                                 <  1:528914/279‑1 (MQ=255)
                                                          ccTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACCGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGACGCACCTTTAATTACCTCTACCCAGcaa                                                                                                                                                                                                                                         >  1:224897/1‑301 (MQ=255)
                                                                tGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTaagatcggaagagca                                                                                                                                                                                                                                   >  2:529587/1‑287 (MQ=255)
                                                                    ggTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGGGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCa                                                                                                                                                                                                                               <  2:453433/301‑1 (MQ=255)
                                                                            gCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCga                                                                                                                                                                                                                       >  2:458382/1‑300 (MQ=255)
                                                                                                     ggcgggcgtgggggtgtatgtAAAGCGTTCGTGCCGGAAGAGTGCTGGGAAGTTGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCa                                                                                                                                                                                             <  1:268641/281‑1 (MQ=255)
                                                                                                      gcgtgcgtgagcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTGGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCat                                                                                                                                                                                            <  1:457818/297‑1 (MQ=255)
                                                                                                                gcgtgaAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGgcg                                                                                                                                                                                  <  2:124629/301‑1 (MQ=255)
                                                                                                                    gtaattaagaggttagttcgggAAGAGTTCTGGGAAGGCGATGCAGGCACGACGACGCCATCGGGTGAAGCGTTGGCGTTACAGAGGACTAATCAGAACGACAAACCGTTCCGTCCGGTAAACAAATAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATGACCGCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGc                                                                                                                                                                              <  1:262750/281‑1 (MQ=255)
                                                                                                                       aTTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGtt                                                                                                                                                                           >  2:290907/1‑301 (MQ=255)
                                                                                                                        ttAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGAGAAAAAAACCATGATGATGGAGCAGCGttt                                                                                                                                                                          >  2:163959/1‑301 (MQ=255)
                                                                                                                         tAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTg                                                                                                                                                                         <  2:164609/301‑1 (MQ=255)
                                                                                                                             gagttcagacgtgtgctcttccgatctgTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATg                                                                                                                                                                     <  2:91409/274‑1 (MQ=255)
                                                                                                                                 tCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACTCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGc                                                                                                                                                                 <  1:108812/301‑1 (MQ=255)
                                                                                                                                  tgTGCCGGAAGAGTTTTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCa                                                                                                                                                                <  2:251452/300‑1 (MQ=255)
                                                                                                                                   gggCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAg                                                                                                                                                               <  1:290505/299‑1 (MQ=255)
                                                                                                                                     gCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCAACCGTCTTGGATTTGGCGTGAAAAAAACCATGATGAAGGCCCACCGTTTGTATGAAACAGGc                                                                                                                                                             >  1:546047/1‑301 (MQ=255)
                                                                                                                                        ggAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCtat                                                                                                                                                          <  2:555250/301‑1 (MQ=255)
                                                                                                                                              gTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGAGTTTGTATGAAGCAGGCTATAACACt                                                                                                                                                    >  1:162945/1‑301 (MQ=255)
                                                                                                                                                  tGGGAAGTCGATGTCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACa                                                                                                                                                <  2:586459/301‑1 (MQ=255)
                                                                                                                                                        gTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGATTGTATGagatcggaagagcgtcgtgtagggaaa                                                                                                                                          >  1:91285/1‑275 (MQ=255)
                                                                                                                                                                     cGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCAATTAAATTAGTAACGACTCCAATaa                                                                                                                             >  1:154822/1‑301 (MQ=255)
                                                                                                                                                                                                  ttACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTc                                                                                                <  2:154973/301‑1 (MQ=255)
                                                                                                                                                                                                              catcaGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCa                                                                                    >  2:149305/1‑301 (MQ=255)
                                                                                                                                                                                                                      cGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGAGGTTAAATCAGCGaaaaa                                                                            >  2:514870/1‑300 (MQ=255)
                                                                                                                                                                                                                                       cgGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATc                                                            <  2:225140/301‑1 (MQ=255)
                                                                                                                                                                                                                                                     cAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCCCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGa                                             >  2:577157/1‑301 (MQ=255)
                                                                                                                                                                                                                                                      aaaCTCATGCTGCGGTAATTCTTCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGaa                                            <  2:419315/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                tGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGcc                                  <  2:163097/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                           tgcgggaaataggggtgtACATCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCGTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGaaaa                       <  1:172773/284‑1 (MQ=255)
                                                                                                                                                                                                                                                                               cgttccaGCACGCTGTAGCGCGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAAACAGTACGCCAGCAAAAAAAAcaac                   >  1:475662/7‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                gCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGc  <  1:514198/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                          
TGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCACAGGAAGCGCACGAAGC  >  NC_000913/1331559‑1332147

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