Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,828,010 C→G 6.1% intergenic (+156/+43) 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,0100CG 6.1% 475.5 / 15.6 147intergenic (+156/+43)mtgA/rpoHmonofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
Reads supporting (aligned to +/- strand):  ref base C (71/66);  new base G (3/6);  total (75/72)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

Only 100 of 147 total aligned reads displayed.

CGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGG  >  NC_002947/5827864‑5828156
                                                                                                                                                   |                                                                                                                                                  
aggcccTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCAcc                                                                                                                                                 >  2:1656054/2‑150 (MQ=255)
    ccTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCTACCGCATCCCCTGTAGAACAGGGCCAAGGCGCTAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCAt                                                                                                                                             >  1:2342885/1‑150 (MQ=255)
     cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                            <  2:282339/150‑1 (MQ=255)
     cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                            <  2:2167900/150‑1 (MQ=255)
     cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGAAGGAGCGGGCTTACCCGCAAAGAGGCCGGGAAAGCCAAACAAAAAAGACGCTCACCCCATa                                                                                                                                            >  2:1898190/1‑150 (MQ=255)
     cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCACTGCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                            <  2:675077/150‑1 (MQ=255)
             ggTAAACCCGCTCCTACGGGAAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                    <  2:1493336/150‑1 (MQ=255)
              gTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGCAGCGGGCTTACCCGCGAAGAGGCCGGCACAGCCAAACAAAAAAGCCGCTCACCCCCTAGGGCCcgcg                                                                                                                                   >  1:2188120/1‑150 (MQ=255)
                   cccGCTCGTACGGCCAACCGCATACCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGAAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGTCGGCAAAGCCAAACATAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttt                                                                                                                              <  2:896741/150‑1 (MQ=255)
                   cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCttt                                                                                                                              >  2:2377660/1‑150 (MQ=255)
                    ccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtttt                                                                                                                             >  2:1001304/1‑150 (MQ=255)
                    ccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtttt                                                                                                                             >  2:1655911/1‑150 (MQ=255)
                       cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                          <  1:1655911/150‑1 (MQ=255)
                       cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                          <  1:2377660/150‑1 (MQ=255)
                        tgcTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCttt                                                                                                                         <  2:1127452/148‑1 (MQ=255)
                         ccTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAGAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgtt                                                                                                                        <  1:1921830/150‑6 (MQ=37)
                         ccTACGGGCAACCGAATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                        <  1:288035/150‑1 (MQ=255)
                          cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCTGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                       <  1:291552/150‑7 (MQ=37)
                          cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                       <  1:1668253/150‑7 (MQ=38)
                              gggCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAcc                                                                                                                   >  2:1020773/1‑150 (MQ=255)
                               ggCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTGCAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAccc                                                                                                                  <  1:214422/150‑1 (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)
                                        gtCCGCGGCGTCCGCGCGGGCGAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAAt                                                                                                         <  2:2188120/149‑1 (MQ=255)
                                        aTCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAAt                                                                                                         <  2:101295/150‑1 (MQ=255)
                                         tCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAAAGCATCACCTCCACGCGCGCCATCACCAGGGAACCCACCTGCAATGTCTAAGAAAAAAGCCGCTCACCCCAGAGGTCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                        <  2:2079817/150‑1 (MQ=255)
                                             cGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGc                                                                                                    >  1:481378/1‑150 (MQ=255)
                                                 ctcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTc                                                                                                                            >  1:945388/1‑122 (MQ=255)
                                                  tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                               >  1:2058368/1‑150 (MQ=255)
                                                  tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                               <  2:1750805/150‑1 (MQ=255)
                                                  tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                               >  1:340725/1‑150 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                            >  1:945998/1‑102 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                            <  2:945998/102‑1 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCaa                                                                                                          <  2:163958/136‑1 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCaa                                                                                                          >  1:163958/1‑136 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCa                                                                                            >  1:1552695/1‑150 (MQ=255)
                                                     tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCa                                                                                            <  2:758360/150‑1 (MQ=255)
                                                        cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                         <  1:1015968/150‑1 (MQ=255)
                                                              aaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAgg                                                                                                                                          >  2:2237229/1‑95 (MQ=255)
                                                              aaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAgg                                                                                                                                          <  1:2237229/95‑1 (MQ=255)
                                                                 ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                        <  1:1452693/110‑1 (MQ=255)
                                                                 ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                        >  2:1452693/1‑110 (MQ=255)
                                                                 ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                >  1:1507857/1‑150 (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)
                                                                           aGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAACCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCGAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGc                                                                      >  2:1737653/1‑150 (MQ=255)
                                                                             gCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAt                                                                    >  1:641149/1‑150 (MQ=255)
                                                                             gCAACAGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCCAAGCCAAACAAAACAGCCGCTCACCCCACAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGCGCCTTGATCTTGTTCATCGcct                                                                    >  2:95102/1‑148 (MQ=255)
                                                                              cAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAtt                                                                   <  1:1638764/150‑1 (MQ=255)
                                                                               aaCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttc                                                                  <  1:95102/150‑1 (MQ=255)
                                                                               aaCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAACCAGACTCAATCAGGCAGAGACCAGTGCCTTGACCTTGTTCACCGCCttc                                                                  >  1:1704380/1‑150 (MQ=255)
                                                                                  cgaaTCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAGGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                               >  1:398706/4‑150 (MQ=255)
                                                                                   gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCTGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                <  1:119048/132‑1 (MQ=255)
                                                                                   gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCTGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                >  2:119048/1‑132 (MQ=255)
                                                                                   gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATtcttct                                                              <  2:2380858/150‑1 (MQ=255)
                                                                                   gCATCCCCTGAAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCACTCAGGCAGCGATCAGTGCCTTGACCTGGTTCATCGCATtcttct                                                              <  2:216278/150‑1 (MQ=255)
                                                                                       ccccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCACCGGCGGCCGGCCCCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCaa                                                          <  2:1268904/150‑1 (MQ=255)
                                                                                         ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                        <  1:1389048/86‑1 (MQ=255)
                                                                                         ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                             <  2:69047/145‑1 (MQ=255)
                                                                                         ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACt                                                        >  1:319262/1‑150 (MQ=255)
                                                                                          cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                       >  2:2276985/1‑150 (MQ=255)
                                                                                            gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttt                                                                                                                             <  2:604601/78‑1 (MQ=18)
                                                                                            gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttt                                                                                                                             >  1:604601/1‑78 (MQ=18)
                                                                                                gAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAAt                                                 >  1:1067328/1‑150 (MQ=255)
                                                                                                      ggCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                            >  1:1661503/1‑150 (MQ=255)
                                                                                                       gtttACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                            <  2:404833/147‑1 (MQ=255)
                                                                                                         ttACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCt                                                                                      >  2:669679/1‑105 (MQ=255)
                                                                                                         ttACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCt                                                                                      <  1:669679/105‑1 (MQ=255)
                                                                                                             ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGa                                     >  1:2173398/1‑150 (MQ=255)
                                                                                                                    gagGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTg                              <  1:1176032/150‑1 (MQ=255)
                                                                                                                         cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTAc                           >  2:712170/1‑148 (MQ=255)
                                                                                                                         cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTAc                           <  1:712170/148‑1 (MQ=255)
                                                                                                                         cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                         <  1:2276985/150‑1 (MQ=255)
                                                                                                                         cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                         <  1:107488/150‑1 (MQ=255)
                                                                                                                          ggCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTg                        >  2:1296866/1‑150 (MQ=255)
                                                                                                                            cAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTc                      >  1:610355/1‑150 (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)
                                                                                                                                gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                   <  2:1617152/85‑1 (MQ=255)
                                                                                                                                 ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                             >  2:2206258/1‑90 (MQ=255)
                                                                                                                                 ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCAATACGCTGTACTTGTCTGCCa                 >  2:858074/1‑150 (MQ=255)
                                                                                                                                  caaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                         >  2:2266073/1‑141 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCTCTCACCCCAGGTGGCCAGCGGCTATTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGGATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt               >  2:762370/1‑150 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                     >  2:574048/1‑112 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt               >  1:2127879/1‑150 (MQ=255)
                                                                                                                                   aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt               >  1:1020321/1‑150 (MQ=255)
                                                                                                                                        aaaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg          <  2:1067328/150‑1 (MQ=255)
                                                                                                                                         aaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGc         >  1:144774/1‑150 (MQ=255)
                                                                                                                                          aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCa        >  1:1927530/1‑150 (MQ=255)
                                                                                                                                             gCCGCTAACCCCATGGGGCCATCGGTTTTGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCg     >  2:1956215/1‑150 (MQ=255)
                                                                                                                                              ccGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTAATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGt    <  2:2058368/150‑1 (MQ=255)
                                                                                                                                                gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATa                                   <  1:2235425/117‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGG  >  NC_002947/5827864‑5828156

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: