Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,242 T→A 12.1% intergenic (‑222/+25) ydiV ← / ← nlpC anti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,792,2420TA12.1% 163.5 / 10.4 61intergenic (‑222/+25)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (26/27);  new base A (6/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAA  >  NC_000913/1791951‑1792539
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         
cTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGAATTAAATGAACGACGATGTATCATCGCCgc                                                                                                                                                                                                                                                                                                  >  1:33877/1‑300 (MQ=255)
 ttGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                                                                                                                                                                                                                                 <  2:319459/301‑1 (MQ=255)
  tgttGGTTGTCTTTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTcc                                                                                                                                                                                                                                                                                                <  1:149246/298‑1 (MQ=255)
     tGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTCCAtt                                                                                                                                                                                                                                                                                             <  2:552511/301‑1 (MQ=255)
       gTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGAACGGCGATGTATCATCGCCGTCCAtttt                                                                                                                                                                                                                                                                                           >  1:162854/1‑301 (MQ=255)
              ctTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACAt                                                                                                                                                                                                                                                                                    >  1:237030/1‑301 (MQ=255)
              ctTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGAATTCAGGAGTAAATGGACGGCGATGTAACATCGCAGGCCATTTTCCGAAAt                                                                                                                                                                                                                                                                                    >  1:189097/1‑301 (MQ=255)
                   cGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTAACATCGACGTCCATTTTTCGACATACTaa                                                                                                                                                                                                                                                                               >  1:140851/1‑300 (MQ=255)
                       aGGAAATAGCTATCAGAATGAGTAATATTTTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTc                                                                                                                                                                                                                                                                           <  2:452128/301‑1 (MQ=255)
                        ggAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCg                                                                                                                                                                                                                                                                          <  2:140976/301‑1 (MQ=255)
                            aTAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGtctc                                                                                                                                                                                                                                                                      <  1:162467/301‑1 (MQ=255)
                               gCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGcc                                                                                                                                                                                                                                                                   >  2:296424/1‑301 (MQ=255)
                                            tAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCGGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGAGAtcggaagagcacacgtctgaactccagtcac                                                                                                                                                                                                                                                      >  2:37540/1‑271 (MQ=255)
                                                               tCTTCATGGTAACCGTCCCTTTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGTGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTGTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTTCCAGAAATTTTTCTGCCAATAGACATTATCa                                                                                                                                                                                                                                  <  2:317839/301‑1 (MQ=255)
                                                                cTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTAGCCTGCCAGAAATTTTTCTGCAAATAGACATTATCa                                                                                                                                                                                                                                  >  2:545467/1‑301 (MQ=255)
                                                                   tattgttaCCGTCCCTTTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGt                                                                                                                                                                                                                               <  2:33928/294‑1 (MQ=255)
                                                                    aTGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGta                                                                                                                                                                                                                              >  2:251472/1‑300 (MQ=255)
                                                                       ttAACCGTTCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAgg                                                                                                                                                                                                                           <  2:446250/300‑1 (MQ=255)
                                                                               cccTCTCGAAGGATTTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATc                                                                                                                                                                                                                   <  1:296020/301‑1 (MQ=255)
                                                                                   ctcGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACta                                                                                                                                                                                                               >  2:417330/1‑300 (MQ=255)
                                                                                   ctcGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTc                                                                                                                                                                                                               >  2:464131/1‑301 (MQ=255)
                                                                                            tttttttCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTgg                                                                                                                                                                                                      <  2:163006/299‑1 (MQ=255)
                                                                                                    aGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACACCCTTGCCGGTAAAGGCg                                                                                                                                                                                              >  1:164027/1‑301 (MQ=255)
                                                                                                       aTCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACCTTATCAAGTGAGGAACGCATCACCCCCTTCCTGGTAACGGCGttg                                                                                                                                                                                           >  1:575038/1‑299 (MQ=255)
                                                                                                           ggCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGCAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCAACACTCCCTTGCTGGTAACGGCGTGGATaa                                                                                                                                                                                       >  1:564606/1‑301 (MQ=255)
                                                                                                               aaaaGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAAAGGAGTGGATAACTTg                                                                                                                                                                                   >  2:351962/1‑301 (MQ=255)
                                                                                                                   gggATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTTTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAACCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGgttg                                                                                                                                                                               <  1:1238/301‑1 (MQ=255)
                                                                                                                       tGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTgg                                                                                                                                                                           <  2:237296/301‑1 (MQ=255)
                                                                                                                          ggCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAAAGGCCTGGATAAATTGGTTGTTGGTAt                                                                                                                                                                        >  1:280132/1‑301 (MQ=255)
                                                                                                                            ccATTAGACATAACAATTACCGAATGATTTACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCg                                                                                                                                                                      <  1:416767/301‑1 (MQ=255)
                                                                                                                                       aaCAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTAGGTAtttacgataccaac                                                                                                                                                           >  1:284455/1‑288 (MQ=255)
                                                                                                                                              tcAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTAGGTAtttacgataccaacaaccaat                                                                                                                                                    <  2:284758/300‑21 (MQ=255)
                                                                                                                                                    ggAGTAACTTTCTTTTTTTCCATGTTTACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCAttt                                                                                                                                              <  1:251205/301‑1 (MQ=255)
                                                                                                                                                      aGTTACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACGTGCAAACCATTTTg                                                                                                                                            <  2:329604/301‑1 (MQ=255)
                                                                                                                                                      aGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCACTAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTTGTATCGTAAATACCTACATGCAAACCAttttt                                                                                                                                            >  1:383892/1‑300 (MQ=255)
                                                                                                                                                          aCTTTCATTTTTTCCATTTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCg                                                                                                                                        <  2:384388/301‑1 (MQ=255)
                                                                                                                                                                    gTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACCTGCAAACCATTTTATCCGGACCCCGttt                                                                                                                              >  1:62485/1‑301 (MQ=255)
                                                                                                                                                                        cATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGaa                                                                                                                          <  2:164179/301‑1 (MQ=255)
                                                                                                                                                                                ccACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATATTGTCCGGACCCCGATTTGAAAAAGACCa                                                                                                                  >  2:128546/1‑301 (MQ=255)
                                                                                                                                                                                        tCAGGGGTCCTGTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCt                                                                                                          <  2:565320/301‑1 (MQ=255)
                                                                                                                                                                                                            ttttCCTTATTATTAACTAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTCCTGTATTAGGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATc                                                                                      <  1:35693/298‑1 (MQ=255)
                                                                                                                                                                                                             tATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATTGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCa                                                                                     <  2:189339/301‑1 (MQ=255)
                                                                                                                                                                                                                      ttcttccgatttatCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCGCCTGGCAGCAACTCGTCTTTATCAATTTGCGTg                                                                            <  2:112950/289‑1 (MQ=255)
                                                                                                                                                                                                                      tataAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTg                                                                            <  2:62553/301‑1 (MQ=255)
                                                                                                                                                                                                                          aaCGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGTACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCCTCCCaa                                                                        >  1:348885/1‑301 (MQ=255)
                                                                                                                                                                                                                                  atCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTGTGAAAAAGACCAGGTCGCCTGGCAGCAACTCGTCTTTATAAATTTGCGTgagatcggaagag                                                                >  1:112826/1‑289 (MQ=255)
                                                                                                                                                                                                                                             taaTTTTTGTCTGCTATATTCTTCTGTATTCAGGAGTAAATGTACGGCTATGTATCATCGTCGTCCATTTTTCGAAATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCGGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGt                                                     <  1:129653/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gtatgtcgaaaaaTGGACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATAACTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGATTCACAGTGCAACTACAAAACGaaaaa                          >  2:418615/11‑299 (MQ=255)
                                                                                                                                                                                                                                                                              tgtgtgctcttccgatctGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATcgcg                    <  2:377057/284‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGCTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTCAGTTTCACGGGGCAGCTGCAAATCGAAACGaacccg                    >  1:22312/1‑296 (MQ=255)
                                                                                                                                                                                                                                                                               aGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGCCATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGCCCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGAAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCAAAACAATCacgc                   >  2:352636/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                   tttGGACGGCGATGTATCATTGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCg               <  2:172991/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     tGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGATTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGAACgcgccacgag             >  1:172794/1‑295 (MQ=255)
                                                                                                                                                                                                                                                                                      ggACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTca            >  2:400375/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                      ggACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTTCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAAAGAACGCGCATCGtac            >  1:443210/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                      ggACGGCGATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGAACGCGCATCGTcc            >  1:498638/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      ggACGGCAATGATACATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAAAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAAACGAAACGATCGCGCAACGTcc            >  1:274383/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                        aCGGCGATGTATCATTGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTcaca          <  2:575753/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             taTGTATCATCGCCGTTCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACcac     <  2:280435/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                             gATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCCGCAACTCGTCTTTATCAATTTGCGTGCCCATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCACGCATCGTCAAcaccac     >  1:486703/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                gTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCAAAACGATcgcgagataggaaaagcggcgg  >  1:376561/1‑283 (MQ=255)
                                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                                         
CTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAA  >  NC_000913/1791951‑1792539

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