Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,236,235 G→A intergenic (+87/‑27) yjbE → / → yjbF uncharacterized protein YjbE/lipoprotein YjbF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,236,2350GA98.3% 202.2 / NA 60intergenic (+87/‑27)yjbE/yjbFuncharacterized protein YjbE/lipoprotein YjbF
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (37/22);  total (38/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACAT  >  NC_000913/4235956‑4236525
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  
ttCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGgaga                                                                                                                                                                                                                                                                                <  1:383112/300‑1 (MQ=255)
  cTGCGTGGGCTGCACCTTTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGagaag                                                                                                                                                                                                                                                                              <  2:381548/300‑1 (MQ=255)
     cGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGNTTCGGCGGGGAGTTCCTCCGCGACCATGGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGTAGAAGAGt                                                                                                                                                                                                                                                                           <  2:526879/300‑1 (MQ=255)
         ggCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCTGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTg                                                                                                                                                                                                                                                                       >  1:57262/1‑300 (MQ=255)
         ggCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTg                                                                                                                                                                                                                                                                       >  1:473726/1‑300 (MQ=255)
         ggCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCAGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCTGTGGCGGTGATGGTTCTAATACCGGGACCACCAAAAACACGACCACCAGTACCCAGTAATAAAGTATGTATCCACAAAATAATTCGAGTCATTGCATCTGTGGCTAGACGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTcatg                                                                                                                                                                                                                                                                       >  1:436392/1‑297 (MQ=255)
             gCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCAGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGc                                                                                                                                                                                                                                                                   >  1:672561/1‑300 (MQ=255)
              cACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGCGAAGAGTCGTGAAGCg                                                                                                                                                                                                                                                                  >  1:335075/1‑300 (MQ=255)
              cACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCACCTATGGAGAAGAGTCGTGAAGCg                                                                                                                                                                                                                                                                  >  1:217854/1‑300 (MQ=255)
                ccTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGAc                                                                                                                                                                                                                                                                >  1:175822/1‑300 (MQ=255)
                 cTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTCATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACATTCAGGTGTGATTATTCTGCCCATCTGGAGAAGAATCGTGAAGCggcc                                                                                                                                                                                                                                                               >  2:228534/1‑297 (MQ=255)
                 cTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCACTCTGGAGAAGAGTCGTGAAGCGAcc                                                                                                                                                                                                                                                               >  2:549004/1‑300 (MQ=255)
                 cTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTACGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACTGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCACATCTGGAGAAGAGTCGTGAAGCGaac                                                                                                                                                                                                                                                               >  2:288019/1‑298 (MQ=255)
                                 aGCGGCAGGGTCGTCAGCAACGTCGGTTTCGTCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGTTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATGAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGATAAGAGTCGTGAAGCGACCTGCACTCATTCTTATc                                                                                                                                                                                                                                               <  2:701670/300‑1 (MQ=255)
                                                 gcaACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAAACACAATCCGGTGTGGTTATTCTGCCCCTCTGGAGGAGAGTCGTTAAGCGACCTGCACTCATTCCTATCTGCCTGTTATTAAAgg                                                                                                                                                                                                                               >  2:35336/1‑300 (MQ=255)
                                                   aaCGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGcc                                                                                                                                                                                                                             <  1:312595/300‑1 (MQ=255)
                                                        cGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTc                                                                                                                                                                                                                        <  1:619145/300‑1 (MQ=255)
                                                             tCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGTTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCa                                                                                                                                                                                                                   <  1:82542/300‑1 (MQ=255)
                                                                gcggGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAGTTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAANCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAACCACGa                                                                                                                                                                                                                >  2:624764/1‑300 (MQ=255)
                                                                     gAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGGAGTATGAAGGGATTAACCATAACCACACTCAGGTGTGATTATTCTGCCCCTCTGGAGGAGAGTCGTGAAGCGAACTGCACTCATTCTTATCTGCCCGTTATTACAGGCCTATTCAACCCCGACCaaa                                                                                                                                                                                                           >  2:104310/1‑300 (MQ=255)
                                                                     gAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTaaa                                                                                                                                                                                                           >  1:236436/1‑300 (MQ=255)
                                                                                     aGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCAc                                                                                                                                                                                           >  1:77524/1‑300 (MQ=255)
                                                                                     aGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCCGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCAc                                                                                                                                                                                           <  2:672326/300‑1 (MQ=255)
                                                                                       cGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTCTCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACtg                                                                                                                                                                                         >  1:554227/1‑300 (MQ=255)
                                                                                       cGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCACCTCTGGAGAAGAGTAGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCAccg                                                                                                                                                                                         >  2:452118/1‑298 (MQ=255)
                                                                                             cACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGTTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGTAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGTGCAATTCACTGTGGGAc                                                                                                                                                                                   <  2:386226/300‑1 (MQ=255)
                                                                                               ccGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTCTCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAg                                                                                                                                                                                 >  1:392984/1‑300 (MQ=255)
                                                                                               ccGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCCGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAg                                                                                                                                                                                 >  1:466279/1‑300 (MQ=255)
                                                                                                         ggcggTGGGTGTCGCGCTTGCGGCAACCGGTGGCGTTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGc                                                                                                                                                                       <  1:700231/300‑1 (MQ=255)
                                                                                                          gcggTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACACAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCa                                                                                                                                                                      >  2:399030/1‑300 (MQ=255)
                                                                                                                        gcTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCNTCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAg                                                                                                                                                        <  1:499064/300‑1 (MQ=255)
                                                                                                                         cTTGCGGCAACCGTTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGc                                                                                                                                                       <  2:694810/300‑1 (MQ=255)
                                                                                                                                                 ggTTCTAATACCGGGACCACCACAACCACGACCAACAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGAATAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAata                                                                                                                               >  1:190188/1‑300 (MQ=255)
                                                                                                                                                  gTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAatat                                                                                                                              >  1:223188/1‑300 (MQ=255)
                                                                                                                                                   ttCTAATACCGGGACCACCACAGCCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATACCCACCCTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGGACTGTACTCATTCTTATCTGTCTGTTATTAAAGGCCTGTTCAGCCACGACTAAAGAGATGGGCAATTCACTGTGGGACAGTCTTTTTGGCACGCCAGACGTATAGCTGACGGCTGGTGACGTTCAGAATATg                                                                                                                             >  2:9733/1‑300 (MQ=255)
                                                                                                                                                   ttCTAATACCGGGACCACCACAACCACGACCGCCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTAAAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGGTGATATTCAAAATATg                                                                                                                             >  2:8863/1‑300 (MQ=255)
                                                                                                                                                   ttCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGAAGGATGATGATATTCAAAATATg                                                                                                                             >  2:690829/1‑300 (MQ=255)
                                                                                                                                                           ccGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACgc                                                                                                                     <  1:473310/300‑1 (MQ=255)
                                                                                                                                                                   accacAACCACGACCACCAGTACCCCGTAATAACGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGAAATTCCATGTGGGACAGTCTGTTCGGCACGCCAGGCGTCCAGCTGGCGGATGATGATATTCAAAATATGCCCTACgccaactagg                                                                                                             >  2:185461/1‑294 (MQ=255)
                                                                                                                                                                   accacAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACCACTGACGGATGATGATCTTTAAAATATGCCCTACGCCAGCCAGt                                                                                                             >  1:317343/1‑300 (MQ=255)
                                                                                                                                                                         aCCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGATTCATTGCATCTGTGGCTAGAAGTATGAAGCGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCATGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGc                                                                                                       <  2:430369/300‑1 (MQ=255)
                                                                                                                                                                         aCCACGACCACCAGTACCCAGGAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTTTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTTGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGc                                                                                                       <  2:57131/300‑1 (MQ=255)
                                                                                                                                                                                             gTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCCGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCCAAATATGCCCTAAGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAg                                                                                   >  1:680570/1‑300 (MQ=255)
                                                                                                                                                                                             gTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAg                                                                                   >  1:549296/1‑300 (MQ=255)
                                                                                                                                                                                                   aaGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCACTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTTGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTAttt                                                                             <  1:415673/300‑1 (MQ=255)
                                                                                                                                                                                                     gtatgtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCCGTACATGCAGCTTAATGGCGGGCCGCAGTTATTtgt                                                                           >  1:45550/1‑300 (MQ=255)
                                                                                                                                                                                                                  aaaTAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTc                                                              >  1:32673/1‑300 (MQ=255)
                                                                                                                                                                                                                         tCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTACCCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGCAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTTAAg                                                       <  2:613529/300‑1 (MQ=255)
                                                                                                                                                                                                                           gAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTGCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGAc                                                     <  1:644386/300‑1 (MQ=255)
                                                                                                                                                                                                                                         gtgGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATg                                       >  1:399255/1‑300 (MQ=255)
                                                                                                                                                                                                                                                gAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTTAGCCACGACTAAAGAGTTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATACGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGAAGGACAACAAAAATGGGTCACt                                >  2:311342/1‑300 (MQ=255)
                                                                                                                                                                                                                                                           ggATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGTACTCATTCTTATCTGCCTGTTATTACAGGCCCGTTCAACCACGAATAAAGAGCTGGGCAATTCAATGTGGGGCAGTCTGTTCGGCAAGACAGGCGTACNGCTGACGGATGATGATATTCAAAATATGCCCTACACCAGCCAGTACATGCAACTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGATGAAAACGGACAACAAAAATGGGTCACTCAGGATCAGGc                     >  2:174600/1‑300 (MQ=255)
                                                                                                                                                                                                                                                           ggATTAACCATAACCACACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGc                     <  1:568068/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaccataaccaCACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCCGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATNCCCTCCGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTActc                >  1:482079/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                 accataaccaCACTCCGGTGTGATTATTCTGCCCCTCTGGATAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTActct               <  2:415942/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ccataaccaCACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCACTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGAATAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCCGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTCATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGCAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTCCtctc              >  1:706643/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                       accaCACTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGNACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTNCGCTGAAGACGGAAAACAAAAATGGGTCACTCAGGATCAGGCTCATCTCGTTAc         >  2:663327/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                          acacTCCGGTGTGATTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACCGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCa      >  2:675755/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                              tCCGGTGTGATTATTCTGCCCCTCTGTAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTTGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATTGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACAt  <  1:238665/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                              tCCGGTGTGATTATTCTGCCCCTCTGTAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACAt  <  1:302709/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  
TTCTGCGTGGGCTGCACCTGTACAGGTGGGCGAAGCGGCAGGGTCGGCAGCAACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACAGCTGACGGATGATGATATTCAAAATATGCCCTACGCCAGCCAGTACATGCAGCTTAATGGCGGGCCGCAGTTATTTGTGGTGCTGGCCTTCGCTGAAGACGGACAACAAAAATGGGTCACTCAGGATCAGGCTACTCTCGTTACCCAACAT  >  NC_000913/4235956‑4236525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: