Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,582,852 G→C 6.9% intergenic (‑75/+51) ypfG ← / ← nudK DUF1176 domain‑containing protein YpfG/GDP‑mannose hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,582,8520GC 6.9% 558.2 / 22.9 171intergenic (‑75/+51)ypfG/nudKDUF1176 domain‑containing protein YpfG/GDP‑mannose hydrolase
Reads supporting (aligned to +/- strand):  ref base G (85/74);  new base C (3/9);  total (88/83)
Fisher's exact test for biased strand distribution p-value = 7.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

Only 100 of 174 total aligned reads displayed.

TTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAA  >  NC_000913/2582703‑2583001
                                                                                                                                                     |                                                                                                                                                     
ttGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCg                                                                                                                                                       <  2:323662/150‑1 (MQ=255)
     cTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGCAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGAttt                                                                                                                                                  >  2:1391718/1‑150 (MQ=255)
     cTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGAttt                                                                                                                                                  <  2:1950732/150‑1 (MQ=255)
      tGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTg                                                                                                                                                 <  1:1226975/150‑1 (MQ=255)
         tGCCGCCCACACCAAAGACGTCGGCAACAGAGCGAAAAAAAAGAGAACGAAAATGCGATAGCGCATGGACGGTCCGACCCCAGAAGTTCAAACAAACAGCACAAGTCGGATCTTCGTTGCGAAGGTAGCGCGGCTCAATCCGATTGGTCg                                                                                                                                              >  1:2249277/1‑150 (MQ=255)
          gccgccCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTcag                                                                                                                                             <  2:74065/150‑3 (MQ=255)
          gccgccCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCgg                                                                                                                                             <  2:1655962/150‑1 (MQ=255)
          gccgccCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCgg                                                                                                                                             >  1:1404876/1‑150 (MQ=255)
           ccgccCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGa                                                                                                                                            >  1:271866/1‑150 (MQ=255)
            cgccCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGAt                                                                                                                                           >  1:8995/1‑150 (MQ=255)
            cgccCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGAt                                                                                                                                           <  2:624505/150‑1 (MQ=255)
            cgccCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGAt                                                                                                                                           >  1:720254/1‑150 (MQ=255)
              cccACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTg                                                                                                                                         <  1:832859/150‑1 (MQ=255)
                cacaCCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAg                                                                                                                                       <  2:1423732/150‑1 (MQ=255)
                 acacCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGc                                                                                                                                      <  2:15567/150‑1 (MQ=255)
                  cacCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGcc                                                                                                                                     <  1:651605/150‑1 (MQ=255)
                  cacCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCCGACAAATCCGATTGAGcc                                                                                                                                     <  1:399264/150‑1 (MQ=255)
                  cacCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCCGACAAATCCGATTGAGcc                                                                                                                                     <  2:1176300/150‑1 (MQ=255)
                  cacCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCCGACAAATCCGATTGAGcc                                                                                                                                     <  1:563116/150‑1 (MQ=255)
                  cacCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACACCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGcc                                                                                                                                     <  2:1609408/150‑1 (MQ=255)
                   acCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCCGACAAATCCGATTGAGccg                                                                                                                                    <  1:226979/150‑2 (MQ=255)
                   acCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCCGACAAATCCGATTGAGccg                                                                                                                                    <  1:218336/150‑2 (MQ=255)
                   acCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCCGACAAATCCGATTGAGccg                                                                                                                                    <  1:1707753/150‑2 (MQ=255)
                    ccAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCaa                                                                                                                                   >  1:949508/1‑150 (MQ=255)
                    ccAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCaa                                                                                                                                   >  2:1018453/1‑150 (MQ=255)
                      aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCaata                                                                                                                                 >  1:1621801/1‑150 (MQ=255)
                      aaaGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCaata                                                                                                                                 <  1:1141362/150‑1 (MQ=255)
                           cGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATaaa                                                                                                                            >  1:122762/1‑150 (MQ=255)
                           cGTCGGCAACAAAGCGAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAaac                                                                                                                           <  2:1817049/150‑1 (MQ=255)
                            ggcGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAaac                                                                                                                           <  2:232116/148‑1 (MQ=255)
                              cGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCGATAATAaacaa                                                                                                                         <  2:369354/150‑1 (MQ=255)
                              cGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAaacaa                                                                                                                         >  2:158417/1‑150 (MQ=255)
                              cGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAaacaa                                                                                                                         >  1:231857/1‑150 (MQ=255)
                               ggCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAaacaac                                                                                                                        >  2:798624/1‑150 (MQ=255)
                                gCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACg                                                                                                                       <  2:233283/150‑1 (MQ=255)
                                    caaAGCGAAAAAAAAGAGAAGGAAAATGAGATAGACCATGGACGGTACGCCCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTGAATCGGATTTGTCGGAATGAGGCAATAATAAACAACGaaaa                                                                                                                   >  2:1742491/1‑149 (MQ=255)
                                      aaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAatat                                                                                                                 >  1:2334191/1‑150 (MQ=255)
                                      aaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAatat                                                                                                                 >  1:2381318/1‑150 (MQ=255)
                                      aaGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGATAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAatat                                                                                                                 >  1:1893663/1‑150 (MQ=255)
                                       aGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAAtatt                                                                                                                >  1:1010151/1‑150 (MQ=255)
                                        gCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAAtatta                                                                                                               >  2:1805972/1‑150 (MQ=255)
                                             aaaaaaGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGtt                                                                                                          >  1:109766/1‑150 (MQ=255)
                                              aaaaaGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTa                                                                                                         <  2:1053409/150‑1 (MQ=255)
                                                   gagaAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTAttttt                                                                                                    >  1:841357/1‑150 (MQ=255)
                                                    agaAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTAttttta                                                                                                   >  2:1354786/1‑149 (MQ=255)
                                                    agaAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTc                                                                                                   <  2:1432534/150‑1 (MQ=255)
                                                    agaAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTc                                                                                                   <  2:820651/150‑1 (MQ=255)
                                                      aaGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAg                                                                                                 <  2:334250/150‑1 (MQ=255)
                                                      aaGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAg                                                                                                 <  2:8995/150‑1 (MQ=255)
                                                        cgtccaTGCGCTATCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTc                                                                                               >  2:1260177/6‑150 (MQ=255)
                                                         gAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTcc                                                                                              <  1:1306627/150‑1 (MQ=255)
                                                             aTGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGCCCATTa                                                                                          >  2:58070/1‑150 (MQ=255)
                                                                 gATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAAtg                                                                                      <  1:1360460/150‑1 (MQ=255)
                                                                   tAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAAtgtg                                                                                    <  2:2055041/150‑1 (MQ=255)
                                                                   aaGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAAtgtg                                                                                    <  2:582819/149‑1 (MQ=255)
                                                                    aGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGa                                                                                   >  1:2001732/1‑150 (MQ=255)
                                                                          tGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTg                                                                             >  1:221772/1‑150 (MQ=255)
                                                                           ggACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGc                                                                            >  1:1288058/1‑150 (MQ=255)
                                                                            gACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCa                                                                           <  1:1162241/150‑1 (MQ=255)
                                                                            gACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCa                                                                           <  1:1980513/150‑1 (MQ=255)
                                                                                  aCGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGt                                                                     >  2:465100/1‑150 (MQ=255)
                                                                                   cGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGtt                                                                    >  1:582976/1‑150 (MQ=255)
                                                                                    gACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTa                                                                   <  1:1846360/150‑1 (MQ=255)
                                                                                     aCCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTaa                                                                  >  1:1176891/1‑150 (MQ=255)
                                                                                      ccccAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAg                                                                 <  2:272803/150‑1 (MQ=255)
                                                                                        ccAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCa                                                               >  2:1101040/1‑150 (MQ=255)
                                                                                         cAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCaa                                                              >  1:2048135/1‑150 (MQ=255)
                                                                                             aGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAac                                                          >  2:1086327/1‑150 (MQ=255)
                                                                                              gTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGAAAATAGTTAAGCAATAaca                                                         >  1:379591/1‑150 (MQ=255)
                                                                                               ttCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAacac                                                        <  2:122762/150‑1 (MQ=255)
                                                                                                tCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACAcc                                                       >  2:1524923/1‑150 (MQ=255)
                                                                                                 caaacaaaCACCACAAGTCCTATCTTCGTAGCCCACGAAGCCCCGCTCAATCGGATATGTAGTATTTAGCCAATAATAAACAACGAAATATTAATGTTATATTTCAGACAATTAAATGTGACGTTTGCAAAAAGTTAAACAATAACACCg                                                      >  2:2359228/1‑150 (MQ=255)
                                                                                                       aacaGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATTTTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTAc                                                <  2:2309196/150‑1 (MQ=255)
                                                                                                         caGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCg                                              <  2:1804910/150‑1 (MQ=255)
                                                                                                         caGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCg                                              >  2:130449/1‑150 (MQ=255)
                                                                                                           gCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTc                                            >  1:1655969/1‑150 (MQ=255)
                                                                                                            cacaAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCa                                           <  1:1610313/150‑1 (MQ=255)
                                                                                                                aGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTa                                       <  2:2488920/150‑1 (MQ=255)
                                                                                                                 gTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAt                                      >  1:1020342/1‑150 (MQ=255)
                                                                                                                   cGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtct                                    <  1:47630/150‑1 (MQ=255)
                                                                                                                         ttCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCgg                              <  2:260567/150‑1 (MQ=255)
                                                                                                                              tGCGCAGATAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTg                         <  2:1893663/150‑1 (MQ=255)
                                                                                                                                cgcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGat                       <  1:158417/150‑1 (MQ=255)
                                                                                                                                 gcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGAGGTTTGCAACTAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGatc                      <  1:1550944/150‑1 (MQ=255)
                                                                                                                                      tAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCAtctc                 >  2:2406644/1‑150 (MQ=255)
                                                                                                                                       aGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTcc                <  1:1010971/150‑1 (MQ=255)
                                                                                                                                          gcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAgc             >  2:2137539/1‑150 (MQ=255)
                                                                                                                                          gcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAgc             >  1:892273/1‑150 (MQ=255)
                                                                                                                                          gcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAgc             <  1:1827314/150‑1 (MQ=255)
                                                                                                                                           tggCTCAATCCGACAAATCCGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACTTATCTCGCCGGTTTTGATCATCTCCAgcg            >  1:442809/2‑150 (MQ=255)
                                                                                                                                           cgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAgcg            <  1:788716/150‑1 (MQ=255)
                                                                                                                                             gCTCAATCGAATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGcc          >  2:1935237/1‑150 (MQ=255)
                                                                                                                                               tCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTg        <  2:1937388/150‑1 (MQ=255)
                                                                                                                                               tCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTg        <  1:798624/150‑1 (MQ=255)
                                                                                                                                                 aaTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGc      <  2:2242487/150‑1 (MQ=255)
                                                                                                                                                 aaTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGc      <  1:2082808/150‑1 (MQ=255)
                                                                                                                                                  aTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCt     >  2:289763/1‑150 (MQ=255)
                                                                                                                                                  aTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCt     <  2:841357/150‑1 (MQ=255)
                                                                                                                                                   tCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTg    >  1:2164096/1‑150 (MQ=255)
                                                                                                                                                     ggATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGaa  >  1:2483624/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGCGCTGGTGCCGCCCACACCAAAGACGTCGGCAACAAAGCGAAAAAAAAGAGAAGGAAAATGCGATAGCGCATGGACGGTACGACCCCAGAAGTTCAAACAAACAGCACAAGTCGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAA  >  NC_000913/2582703‑2583001

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