Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,379,259 T→C K247E (AAA→GAA)  zapE ← divisome ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,379,2590TC98.2% 358.3 / ‑2.8 110K247E (AAA→GAA) zapEdivisome ATPase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (60/47);  total (62/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 113 total aligned reads displayed.

CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCA  >  NC_000913/3379162‑3379354
                                                                                                  |                                                                                                
 acaGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑c‑ccc                                                                                              >  2:552195/1‑100 (MQ=255)
    gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTGTAACGTCGGTGAATTTTCTCGTT‑C‑cccccc                                                                                           >  1:2518525/1‑100 (MQ=255)
    gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑cccccc                                                                                           >  1:772975/1‑100 (MQ=255)
    gCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑cccccc                                                                                           >  2:2007106/1‑100 (MQ=255)
      gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCg                                                                                         <  2:1636645/100‑1 (MQ=255)
         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCa                                                                                      >  1:2472459/1‑100 (MQ=255)
         gtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCa                                                                                      <  2:1335850/100‑1 (MQ=255)
         gtaAAAGAGAACGCCAGCGTCTGGTTATCAACGCCAATTGTTGCTAATGGCCGAAGGTTGATTGCTAAAATCGGCGAATTTTCTCCCT‑C‑CCCCCCCGCCa                                                                                      >  2:2199869/1‑100 (MQ=255)
           aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTCC‑CCCCCCCGCCAg                                                                                     >  2:1510614/1‑100 (MQ=255)
           aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAgc                                                                                    >  1:2521901/1‑100 (MQ=255)
            aaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAgcg                                                                                   <  1:2384202/100‑1 (MQ=255)
              agagACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGcc                                                                                 >  2:713920/1‑100 (MQ=255)
                 gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAAc                                                                              <  2:1553875/100‑1 (MQ=255)
                  aCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAAcc                                                                             <  1:1936519/100‑1 (MQ=255)
                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCGCCCCGCCAGCGCCAACCa                                                                            >  2:1029276/1‑100 (MQ=255)
                   ccGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCa                                                                            >  1:1427070/1‑100 (MQ=255)
                    cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCAt                                                                           >  1:718352/1‑100 (MQ=255)
                    cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCAt                                                                           <  2:187845/100‑1 (MQ=255)
                    cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCAt                                                                           >  2:2062255/1‑100 (MQ=255)
                      ccAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAg                                                                         <  2:2464604/100‑1 (MQ=255)
                       naGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCCCCAACCATAGt                                                                        >  1:2110509/2‑100 (MQ=255)
                       gaGCGTCTGGTTCACGACGCCCAGTGCTGCAAATGGCCGATGGTTGATGTCTAACGCCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGt                                                                        <  2:834637/99‑1 (MQ=255)
                         gCGTCTGGTTCTCGACGCCGATTGTTGGGAATGGCCCGTGGTTGATGTCTAAGGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGttt                                                                      <  2:497139/100‑1 (MQ=255)
                         gCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGttt                                                                      >  2:1179772/1‑100 (MQ=255)
                           gTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTAt                                                                    <  1:1197545/100‑1 (MQ=255)
                               ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                                >  1:1807212/1‑100 (MQ=255)
                               ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                                >  1:1615903/1‑100 (MQ=255)
                               ggTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAt                                                                >  2:375545/1‑100 (MQ=255)
                                gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCAtt                                                               >  1:1487793/1‑100 (MQ=255)
                                  tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                             <  1:714295/100‑1 (MQ=255)
                                  tctcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTg                                                             <  1:1482915/100‑1 (MQ=255)
                                    tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                           >  2:776698/1‑100 (MQ=255)
                                    tcGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcg                                                           <  1:2130315/100‑1 (MQ=255)
                                     cGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGGCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                          >  2:200870/1‑100 (MQ=255)
                                       aCGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGccc                                                        <  1:1564124/100‑1 (MQ=255)
                                         gCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCgg                                                      >  2:2458538/1‑100 (MQ=255)
                                          cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                     >  2:2025525/1‑100 (MQ=255)
                                          cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCggg                                                     >  2:1422009/1‑100 (MQ=255)
                                           ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGt                                                    <  2:340604/100‑1 (MQ=255)
                                             aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                  >  2:2249519/1‑100 (MQ=255)
                                             aTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGttt                                                  <  2:799201/100‑1 (MQ=255)
                                               tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCCACCATAGTTTATCCATTTGCGCCCGGGGTTc                                                 >  1:2565264/1‑100 (MQ=255)
                                               tgttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCa                                                >  1:1652451/1‑100 (MQ=255)
                                                gttgGTAGTGGCCGATGGTTGATTTGTAACGTCGGTCAATATTGGCGTT‑G‑CCCCGCCGCCAACGCCCACTATAGTTTATCCATTTGCGCCCGGATTTCAt                                               >  2:2440728/1‑100 (MQ=255)
                                                gttgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCGCCCCGCCAGCGCCAACCATAGTTTATCCATTTTCGCCCGGGTTTCAt                                               >  2:70368/1‑100 (MQ=255)
                                                  tgCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                             >  1:1887846/1‑100 (MQ=255)
                                                   gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCg                                             <  2:33176/100‑1 (MQ=255)
                                                   gCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGt                                            >  2:1830087/1‑100 (MQ=255)
                                                    cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                           >  1:2144406/1‑100 (MQ=255)
                                                    cTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTg                                           >  1:2234047/1‑100 (MQ=255)
                                                    cTAATGGCCGATGGTTGATTTCT‑ACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGATTATCCATTTGCGCCCGGGTTTCATCGtaa                                          >  2:1480927/1‑98 (MQ=255)
                                                     tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                          >  2:1013688/1‑100 (MQ=255)
                                                      aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTCTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGaa                                         <  1:1410404/100‑1 (MQ=255)
                                                        tGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTAAGCATTAGCGCCCGGGTTTCATCGTGAAg                                        >  1:2404942/1‑100 (MQ=255)
                                                         ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTg                                      <  2:2269714/100‑1 (MQ=255)
                                                           ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                    >  1:1266537/1‑100 (MQ=255)
                                                           ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGc                                    >  1:2008103/1‑100 (MQ=255)
                                                            cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                   <  1:1837155/100‑1 (MQ=255)
                                                            cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                   <  1:1565263/100‑1 (MQ=255)
                                                            cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                   <  1:1837121/100‑1 (MQ=255)
                                                             gATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCg                                   <  1:1531755/100‑1 (MQ=255)
                                                               tGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGaaa                                <  1:1075242/100‑1 (MQ=255)
                                                                ggTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                               >  2:432024/1‑100 (MQ=255)
                                                                ggTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                               <  1:1525478/100‑1 (MQ=255)
                                                                ggTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                               <  1:686636/100‑1 (MQ=255)
                                                                 gTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGc                              <  1:582242/100‑1 (MQ=255)
                                                                  ttGATTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGcc                             >  1:259713/1‑100 (MQ=255)
                                                                   tGATTTCTACCCCCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                            <  2:450243/100‑1 (MQ=255)
                                                                   tGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                            <  1:523545/100‑1 (MQ=255)
                                                                    gATTTCTAACGTCGGTGAATTTTCTCGTT‑CCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCca                            >  1:151326/1‑100 (MQ=255)
                                                                    gATTTCTAACGTCGGTGAATTGTCTCGTC‑C‑CCCCCCCTCCAGCGCCAACCATAGTTTATCCATTTGCGCCTGGGTATCATCGTGAAGTGGCGAAAGCcac                           <  1:2462879/100‑1 (MQ=255)
                                                                     aTTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                          <  1:2350977/100‑1 (MQ=255)
                                                                     aTTTCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCcaca                          <  1:2319581/100‑1 (MQ=255)
                                                                      tttCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                         <  2:1137875/100‑1 (MQ=255)
                                                                      tttCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                         >  2:1003524/1‑100 (MQ=255)
                                                                       ttCTAACGTCGGTGAATTTTCTCGTT‑C‑CACCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGa                        <  1:2398436/100‑1 (MQ=255)
                                                                        tCTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                       >  2:1061312/1‑100 (MQ=255)
                                                                         cTAACGTCGGTGAATTTTCTCGTTCC‑CCCCCCCGCCAGCGCCAACCCTAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAATGGCGAAAGCCACAAAt                       >  1:1644852/1‑100 (MQ=255)
                                                                         cTAACGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATg                      <  2:1022587/100‑1 (MQ=255)
                                                                           aaCGTCGGTGCATTTTCTCGTT‑C‑CCCCCCCGCCAGAGCCCACCATAGTTTATCCATTCGCGCCCGGGTTTCGACGGGAAGTGGCGAAAGGCACAGATgcg                    >  2:1559244/1‑100 (MQ=255)
                                                                             cGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCTAACAATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                  >  2:2436532/1‑100 (MQ=255)
                                                                             cGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGGAAGCCACAGATGCGcc                  >  1:2001335/1‑100 (MQ=255)
                                                                             cGTCGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGcc                  >  2:2464243/1‑100 (MQ=255)
                                                                                cGGTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga               <  2:756767/100‑1 (MQ=255)
                                                                                 ggTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag              >  2:1553654/1‑100 (MQ=255)
                                                                                 ggTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag              <  2:157632/100‑1 (MQ=255)
                                                                                 ggTGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag              <  1:889426/100‑1 (MQ=255)
                                                                                   tGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg            >  1:28850/1‑100 (MQ=255)
                                                                                   tGAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtg            <  1:2431755/100‑1 (MQ=255)
                                                                                    gAATTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga           <  1:239220/100‑1 (MQ=255)
                                                                                     aaTTTTCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtgag          <  1:1155743/100‑1 (MQ=255)
                                                                                      aTTTTCTCGTT‑C‑CCCGCCCGCCAGCGCCCACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGCTGCGCCTgagtgagt         >  2:1969279/1‑100 (MQ=255)
                                                                                       ttttCTCGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGagtgagtg        <  2:631607/100‑1 (MQ=255)
                                                                                          tctcGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc     <  1:1835366/100‑1 (MQ=255)
                                                                                          tctcGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTAc     >  1:1343213/1‑100 (MQ=255)
                                                                                           gtcGTT‑C‑CCCCCCCGCCAGCGCCAACGATAGTTTAGCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACGGATGCGCCTGAGTGAGTGTACg    <  2:54359/99‑1 (MQ=255)
                                                                                           ctcGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg    >  1:680241/1‑100 (MQ=255)
                                                                                           ctcGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg    >  1:680222/1‑100 (MQ=255)
                                                                                           ctcGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg    <  1:372978/100‑1 (MQ=255)
                                                                                             cGTT‑C‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCa  >  1:1432198/1‑100 (MQ=255)
                                                                                                  |                                                                                                
CACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑T‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCA  >  NC_000913/3379162‑3379354

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