Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,828,011 A→G 6.1% intergenic (+157/+42) mtgA → / ← rpoH monofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,828,0110AG 6.1% 439.3 / 11.5 147intergenic (+157/+42)mtgA/rpoHmonofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
Reads supporting (aligned to +/- strand):  ref base A (72/66);  new base G (3/6);  total (75/72)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 154 total aligned reads displayed.

CGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTT  >  NC_002947/5827864‑5828160
                                                                                                                                                    |                                                                                                                                                     
aggcccTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCAcc                                                                                                                                                     >  2:1656054/2‑150 (MQ=255)
    ccTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCTACCGCATCCCCTGTAGAACAGGGCCAAGGCGCTAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCAt                                                                                                                                                 >  1:2342885/1‑150 (MQ=255)
     cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                                <  2:282339/150‑1 (MQ=255)
            gggTAAACCCGCTCCTACGGGAAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAg                                                                                                                                         >  1:1493336/1‑150 (MQ=255)
             ggTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                        >  2:1015968/1‑150 (MQ=255)
             ggTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCACCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGAGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                        <  2:764866/150‑1 (MQ=255)
             ggTAAACCCGCTCCTACGGGAAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                        <  2:1493336/150‑1 (MQ=255)
                   cccGCTCGTACGGCCAACCGCATACCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGAAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGTCGGCAAAGCCAAACATAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttt                                                                                                                                  <  2:896741/150‑1 (MQ=255)
                   cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCttt                                                                                                                                  >  1:1439981/1‑150 (MQ=255)
                    ccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtttt                                                                                                                                 >  2:1001304/1‑150 (MQ=255)
                      gCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttttt                                                                                                                               >  2:1668253/1‑148 (MQ=39)
                       cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                              <  1:1655911/150‑1 (MQ=255)
                       cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                              <  1:2377660/150‑1 (MQ=255)
                         ccTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAGAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgtt                                                                                                                            <  1:1921830/150‑6 (MQ=37)
                         ccTACGGGCAACCGAATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                            <  1:288035/150‑1 (MQ=255)
                          cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCTGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                           <  1:291552/150‑7 (MQ=37)
                          cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAGAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                           <  1:1517162/150‑7 (MQ=37)
                          cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                           <  1:1668253/150‑7 (MQ=38)
                            aCGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACaa                                                                                                                         >  2:11468/1‑150 (MQ=255)
                              gggCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAcc                                                                                                                       >  2:1020773/1‑150 (MQ=255)
                                gCAACCGCATCCCCGGCCTCCTCGCGGGCCAAACACCGGCTGGCCATAACCTCAGCCCCTGTAGGAGCGGGAAACCTCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCa                                                                                                                     <  1:1739473/150‑1 (MQ=255)
                                gCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCa                                                                                                                     >  2:1176032/1‑150 (MQ=255)
                                gCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCACCCCCAGCCGGAGAGCCACCTGAACCGAATAGGCCGGCTGAGTCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCa                                                                                                                     <  1:1497440/150‑1 (MQ=255)
                                  aaCCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGa                                                                                                                   >  2:1567176/1‑150 (MQ=255)
                                    ccGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGTCGACCTCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACt                                                                                                                 >  2:2269360/1‑150 (MQ=255)
                                        gtCCGCGGCGTCCGCGCGGGCGAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAAt                                                                                                             <  2:2188120/149‑1 (MQ=255)
                                         tCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                            <  2:2344992/150‑1 (MQ=255)
                                         tCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAAAGCATCACCTCCACGCGCGCCATCACCAGGGAACCCACCTGCAATGTCTAAGAAAAAAGCCGCTCACCCCAGAGGTCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                            <  2:2079817/150‑1 (MQ=255)
                                            ccGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCAGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAgg                                                                                                         >  1:758360/1‑150 (MQ=255)
                                                cctcctCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGc                                                                                                     <  1:624525/150‑1 (MQ=255)
                                                 ctcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTc                                                                                                                                <  2:945388/122‑1 (MQ=255)
                                                  tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                   >  1:340725/1‑150 (MQ=255)
                                                  tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                   >  1:2058368/1‑150 (MQ=255)
                                                  tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                   <  2:1750805/150‑1 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                                >  1:945998/1‑102 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCaa                                                                                                              <  2:163958/136‑1 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCaa                                                                                                              >  1:163958/1‑136 (MQ=255)
                                                        cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                             <  1:1020773/150‑1 (MQ=255)
                                                        cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                             <  1:1015968/150‑1 (MQ=255)
                                                              aaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAgg                                                                                                                                              >  2:2237229/1‑95 (MQ=255)
                                                                    cgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                 <  1:1001304/150‑1 (MQ=255)
                                                                      ccTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTc                                                                               >  1:1240782/1‑150 (MQ=255)
                                                                      ccTACAGGCAACCGCAGCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTc                                                                               <  2:2341995/150‑1 (MQ=255)
                                                                         aCAGGCAACCGCATCCGCTGTAGGAGCGTGCTTGCCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                            <  1:2054126/150‑1 (MQ=255)
                                                                         aCAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                            >  1:369912/1‑150 (MQ=255)
                                                                          cAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGACGAGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                           >  2:2131660/1‑150 (MQ=255)
                                                                             gCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAt                                                                        >  1:641149/1‑150 (MQ=255)
                                                                             gCAACAGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCCAAGCCAAACAAAACAGCCGCTCACCCCACAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGCGCCTTGATCTTGTTCATCGcct                                                                        >  2:95102/1‑148 (MQ=255)
                                                                               aaCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAACCAGACTCAATCAGGCAGAGACCAGTGCCTTGACCTTGTTCACCGCCttc                                                                      >  1:1704380/1‑150 (MQ=255)
                                                                                   gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATtcttct                                                                  <  2:2380858/150‑1 (MQ=255)
                                                                                   gCATCCCCTGAAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCACTCAGGCAGCGATCAGTGCCTTGACCTGGTTCATCGCATtcttct                                                                  <  2:216278/150‑1 (MQ=255)
                                                                                    cATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                                   >  2:2037523/1‑148 (MQ=255)
                                                                                       ccccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCACCGGCGGCCGGCCCCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCaa                                                              <  2:1268904/150‑1 (MQ=255)
                                                                                        cccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAAc                                                             <  1:11468/150‑1 (MQ=255)
                                                                                        cccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAAc                                                             <  2:641149/150‑1 (MQ=255)
                                                                                         ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                            >  2:1389048/1‑86 (MQ=255)
                                                                                         ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                                 <  2:69047/145‑1 (MQ=255)
                                                                                         ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                                 >  1:69047/1‑145 (MQ=255)
                                                                                         ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACt                                                            >  1:319262/1‑150 (MQ=255)
                                                                                          cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                           >  2:107488/1‑150 (MQ=255)
                                                                                            gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttt                                                                                                                                 >  1:604601/1‑78 (MQ=18)
                                                                                                gAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAAt                                                     >  1:1067328/1‑150 (MQ=255)
                                                                                                      ggCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                >  1:1661503/1‑150 (MQ=255)
                                                                                                       gtttACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                >  1:404833/3‑149 (MQ=255)
                                                                                                       gtttACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                <  2:404833/147‑1 (MQ=255)
                                                                                                         ttACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCt                                                                                          >  2:669679/1‑105 (MQ=255)
                                                                                                             ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGaa                                                      >  2:1631884/1‑137 (MQ=255)
                                                                                                             ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGaa                                                      <  1:1631884/137‑1 (MQ=255)
                                                                                                             ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGa                                         >  1:2173398/1‑150 (MQ=255)
                                                                                                                  aaGAGGCCGGCAAAGCCAAAAAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGc                                    <  1:2170254/150‑1 (MQ=255)
                                                                                                                    gagGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTg                                  <  1:1176032/150‑1 (MQ=255)
                                                                                                                      ggCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTa                                <  1:1228511/150‑1 (MQ=255)
                                                                                                                         cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTAc                               <  1:712170/148‑1 (MQ=255)
                                                                                                                           gCAAAGCCAAACAAAAAAGCCTCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGACTACGCCCAGCCGCTACGCTGTACTTGt                           <  2:1240782/150‑1 (MQ=255)
                                                                                                                            cAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTc                          <  2:369912/150‑1 (MQ=255)
                                                                                                                             aaaGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCt                         >  1:729936/1‑150 (MQ=255)
                                                                                                                               aGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGc                       >  2:1104481/1‑150 (MQ=255)
                                                                                                                                gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                            <  1:433262/64‑1 (MQ=255)
                                                                                                                                gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                       >  1:1617152/1‑85 (MQ=255)
                                                                                                                                 ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                 <  1:2206258/90‑1 (MQ=255)
                                                                                                                                 ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                     <  2:2173398/150‑1 (MQ=255)
                                                                                                                                 ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                     >  1:978828/1‑150 (MQ=255)
                                                                                                                                  caaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                             >  2:2266073/1‑141 (MQ=255)
                                                                                                                                  caaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                             <  1:2266073/141‑1 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCTCTCACCCCAGGTGGCCAGCGGCTATTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGGATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt                   >  2:762370/1‑150 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                         <  1:574048/112‑1 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                         >  2:574048/1‑112 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt                   >  2:177312/1‑150 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt                   >  1:2127879/1‑150 (MQ=255)
                                                                                                                                         aaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGc             >  1:144774/1‑150 (MQ=255)
                                                                                                                                          aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTCCCACCCAGACTCAATCAGCCAGCGACAAGTGCCTTAACCTTGTTCATCGCATTCTTCTCCAACTGCCGAATAAGCTCCGCCGCTACGCTCTACTTGTATGCCAGTTCATGcc            >  2:680052/1‑149 (MQ=255)
                                                                                                                                          aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCa            >  1:1927530/1‑150 (MQ=255)
                                                                                                                                           aaGCCGCTGGCCCTATGGGGTGAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAg           >  2:739088/1‑150 (MQ=255)
                                                                                                                                             gCCGCTAACCCCATGGGGCCATCGGTTTTGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCg         >  2:1956215/1‑150 (MQ=255)
                                                                                                                                                gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATa                                       <  1:2235425/117‑1 (MQ=255)
                                                                                                                                                gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATa                                       >  2:2235425/1‑117 (MQ=255)
                                                                                                                                                    aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGt                                                                                              >  1:2237929/1‑58 (MQ=255)
                                                                                                                                                    aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGc                                                 <  2:728453/103‑1 (MQ=255)
                                                                                                                                                    aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCtt  <  2:1439981/150‑1 (MQ=255)
                                                                                                                                                    aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCtt  <  2:1555475/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTT  >  NC_002947/5827864‑5828160

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: