Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,331,839 G→A 9.9% A264A (GCG→GCA 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,8390GA 9.9% 129.4 / 3.6 46A264A (GCG→GCAtopADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base G (17/24);  new base A (2/3);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GCAGGCGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCA  >  NC_000913/1331542‑1332130
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    
gggggggggggtgttgaaggggcgcgGGGGGTGGGATGTGGTTTGGTCGCTGTTATGGGAAAAGGTCGGTCGGGGGCTGTTTGGCGGTCGGGTGCAGTCGGTGGGGGGTCGGCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGGTTGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAgcgcgc                                                                                                                                                                                                                                                                                                   <  1:331559/279‑1 (MQ=255)
   ggCGCGTCGCTTTATTGACCGCGTGGTGGGGTATATGGTTTCGCGGCTGCTATGGAAAAAGATCGCTCGTGGGCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGAATCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTAc                                                                                                                                                                                                                                                                                                <  2:220397/301‑1 (MQ=255)
                 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)
                                                                                                                                                                                                                                                                                                cgttccaGCACGCTGTAGCGCGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAAACAGTACGCCAGCAAAAAAAAcaac  >  1:475662/7‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    
GCAGGCGCGTCGCTTTATGGACCGCGTGGTGGGGTATATGGTTTCGCCGCTGCTATGGAAAAAGATCGCTCGTGGCCTGTCTGCCGGTCGTGTGCAGTCGGTGGCGGTTCGCCTGGTGGTCGAGCGTGAGCGTGAAATTAAAGCGTTCGTGCCGGAAGAGTTCTGGGAAGTCGATGCCAGCACGACCACGCCATCTGGTGAAGCGTTGGCGTTACAGGTGACTCATCAGAACGACAAACCGTTCCGTCCGGTCAACAAAGAACAAACTCAGGCTGCGGTAAGTCTGCTGGAAAAAGCGCGCTACAGCGTGCTGGAACGTGAAGACAAACCGACAACCAGTAAACCTGGCGCTCCTTTTATTACCTCTACGCTGCAACAAGCTGCCAGCACCCGTCTTGGATTTGGCGTGAAAAAAACCATGATGATGGCGCAGCGTTTGTATGAAGCAGGCTATATCACTTACATGCGTACCGACTCCACTAACCTGAGTCAGGACGCGGTAAATATGGTTCGCGGTTATATCAGCGATAATTTTGGTAAGAAATATCTGCCGGAAAGTCCGAATCAGTACGCCAGCAAAGAAAACTCA  >  NC_000913/1331542‑1332130

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: