Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,792,240 T→A 9.5% intergenic (‑220/+27) 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,2400TA 9.5% 181.0 / 8.7 63intergenic (‑220/+27)ydiV/nlpCanti‑FlhD4C2 factor, inactive EAL family phosphodiesterase/putative C40 clan peptidase lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base A (4/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAA  >  NC_000913/1791942‑1792539
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
cGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAACCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCAATGTaa                                                                                                                                                                                                                                                                                                           >  2:1238/1‑300 (MQ=255)
         cTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGAATTAAATGAACGACGATGTATCATCGCCgc                                                                                                                                                                                                                                                                                                  >  1:33877/1‑300 (MQ=255)
          ttGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGATACATCGCCGTc                                                                                                                                                                                                                                                                                                 <  2:319459/301‑1 (MQ=255)
           tgttGGTTGTCTTTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTcc                                                                                                                                                                                                                                                                                                <  1:149246/298‑1 (MQ=255)
            gCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTTTTTGCAAAATCCTCATTAATTCAGCATTAAATGAACGGCACTGTacaaacaccagcct                                                                                                                                                                                                                                                                                               >  1:482353/1‑288 (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)
                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                           
CGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTCAGTATCCGGCTAAAAGGGATGAGGCCATAAGACATAACAATTACAGAAGGAGTAACTTTCATTTGTTCCATGTTAACCACTTTTTCAGGGGTCCTTTTTTTAGATTATCCTGATTATAAACGAATAATCCATTTACGGAATTTTTGTCTGCAAAATACTACTGTATTCAGGAGTAAATGGACGGCGATGTATCATCGCCGTCCATTTTTCGACATACTAGATTCGTCTCGCCTGCCAGAAATTTTTCTGCCAATAGACATTATCAAGTGAGGAACGCATCACTCCCTTGCTGGTAGAGGCGTGGATAAATTGGTTGTTGGTATCGTAAATACCTACATGCAAACCATTTTGTCCGGACCCCGTTTTGAAAAAGACCAGGTCACCTGGCAGCAACTCGTCTTTATCAATTTGCGTGCCAATAGAGGCTTGTTGTTTGGTTTCACGGGGCAGCTGCAAATCGAAACGATCGCGCATCGTCACAACCACAAA  >  NC_000913/1791942‑1792539

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