Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,828,018 A→G 6.8% intergenic (+164/+35) 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,0180AG 6.8% 459.2 / 13.3 153intergenic (+164/+35)mtgA/rpoHmonofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
Reads supporting (aligned to +/- strand):  ref base A (72/70);  new base G (5/6);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 160 total aligned reads displayed.

CTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCC  >  NC_002947/5827869‑5828167
                                                                                                                                                      |                                                                                                                                                     
    cgcgGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGc                                                                                                                                                   <  1:36483/150‑1 (MQ=255)
        ggTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                               >  2:1015968/1‑150 (MQ=255)
        ggTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCACCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGAGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                               <  2:764866/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)
              cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCttt                                                                                                                                         >  1:1439981/1‑150 (MQ=255)
                 gCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttttt                                                                                                                                      >  2:1668253/1‑148 (MQ=39)
                  cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                                     <  1:1655911/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)
                     cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                                  <  1:1668253/150‑7 (MQ=38)
                          ggCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTGCAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAccc                                                                                                                             <  1:214422/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)
                                           cctcctCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGc                                                                                                            <  1:624525/150‑1 (MQ=255)
                                            ctcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTc                                                                                                                                       <  2:945388/122‑1 (MQ=255)
                                            ctcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTc                                                                                                                                       >  1:945388/1‑122 (MQ=255)
                                             tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAgg                                                                                                                <  1:29562/144‑1 (MQ=255)
                                             tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAgg                                                                                                                >  2:29562/1‑144 (MQ=255)
                                             tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                          <  2:1750805/150‑1 (MQ=255)
                                             tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                          >  1:340725/1‑150 (MQ=255)
                                                tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                                       >  1:945998/1‑102 (MQ=255)
                                                tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCaa                                                                                                                     <  2:163958/136‑1 (MQ=255)
                                                tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCa                                                                                                       <  2:758360/150‑1 (MQ=255)
                                                tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCa                                                                                                       >  1:1552695/1‑150 (MQ=255)
                                                   cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                                    <  1:1015968/150‑1 (MQ=255)
                                                   cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                                    <  1:1020773/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)
                                                            ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                           >  1:1507857/1‑150 (MQ=255)
                                                               cgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                        <  1:1001304/150‑1 (MQ=255)
                                                               cgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                        <  2:1552695/150‑1 (MQ=255)
                                                                 ccTACAGGCAACCGCAGCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTc                                                                                      <  2:2341995/150‑1 (MQ=255)
                                                                    aCAGGCAACCGCATCCGCTGTAGGAGCGTGCTTGCCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                                   <  1:2054126/150‑1 (MQ=255)
                                                                     cAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGACGAGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                                  >  2:2131660/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)
                                                                             cgaaTCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAGGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                                          >  1:398706/4‑150 (MQ=255)
                                                                              gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCTGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                           >  2:119048/1‑132 (MQ=255)
                                                                              gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCTGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                           <  1:119048/132‑1 (MQ=255)
                                                                              gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATtcttct                                                                         <  2:2380858/150‑1 (MQ=255)
                                                                              gCATCCCCTGAAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCACTCAGGCAGCGATCAGTGCCTTGACCTGGTTCATCGCATtcttct                                                                         <  2:216278/150‑1 (MQ=255)
                                                                               cATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                                          <  1:2037523/148‑1 (MQ=255)
                                                                                  ccccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCACCGGCGGCCGGCCCCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCaa                                                                     <  2:1268904/150‑1 (MQ=255)
                                                                                   cccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAAc                                                                    <  1:11468/150‑1 (MQ=255)
                                                                                    ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                                   <  1:1389048/86‑1 (MQ=255)
                                                                                    ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                                        <  2:69047/145‑1 (MQ=255)
                                                                                    ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                                        >  1:69047/1‑145 (MQ=255)
                                                                                     cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                                  >  2:107488/1‑150 (MQ=255)
                                                                                     cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                                  >  2:2276985/1‑150 (MQ=255)
                                                                                       gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttt                                                                                                                                        >  1:604601/1‑78 (MQ=18)
                                                                                       gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttt                                                                                                                                        <  2:604601/78‑1 (MQ=18)
                                                                                       gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                <  1:1656054/150‑1 (MQ=255)
                                                                                       gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                <  1:1898190/150‑1 (MQ=255)
                                                                                           gAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAAt                                                            >  1:1067328/1‑150 (MQ=255)
                                                                                                  gtttACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                       <  2:404833/147‑1 (MQ=255)
                                                                                                  gtttACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                       >  1:404833/3‑149 (MQ=255)
                                                                                                    ttACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCt                                                                                                 >  2:669679/1‑105 (MQ=255)
                                                                                                             aaGAGGCCGGCAAAGCCAAAAAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGc                                           <  1:2170254/150‑1 (MQ=255)
                                                                                                                 ggCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTa                                       <  1:1228511/150‑1 (MQ=255)
                                                                                                                    cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCCGGGAACGATTAGTGTTTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                    >  1:2097442/1‑150 (MQ=255)
                                                                                                                    cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                    <  1:2276985/150‑1 (MQ=255)
                                                                                                                     ggCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTg                                   >  2:1296866/1‑150 (MQ=255)
                                                                                                                       cAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTc                                 >  1:610355/1‑150 (MQ=255)
                                                                                                                        aaaGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCt                                >  1:729936/1‑150 (MQ=255)
                                                                                                                           gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                                   >  2:433262/1‑64 (MQ=255)
                                                                                                                           gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                              >  1:1617152/1‑85 (MQ=255)
                                                                                                                           gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                              <  2:1617152/85‑1 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                        <  1:2206258/90‑1 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                            >  1:2399558/1‑150 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                            >  1:978828/1‑150 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                            <  2:2173398/150‑1 (MQ=255)
                                                                                                                             caaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                    >  2:2266073/1‑141 (MQ=255)
                                                                                                                             caaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                    <  1:2266073/141‑1 (MQ=255)
                                                                                                                              aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt                          >  2:177312/1‑150 (MQ=255)
                                                                                                                                   aaaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg                     <  2:1067328/150‑1 (MQ=255)
                                                                                                                                     aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTCCCACCCAGACTCAATCAGCCAGCGACAAGTGCCTTAACCTTGTTCATCGCATTCTTCTCCAACTGCCGAATAAGCTCCGCCGCTACGCTCTACTTGTATGCCAGTTCATGcc                   >  2:680052/1‑149 (MQ=255)
                                                                                                                                        gCCGCTAACCCCATGGGGCCATCGGTTTTGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCg                >  2:1956215/1‑150 (MQ=255)
                                                                                                                                         ccGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTAATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGt               <  2:2058368/150‑1 (MQ=255)
                                                                                                                                           gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATa                                              >  2:2235425/1‑117 (MQ=255)
                                                                                                                                           gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTgg             >  1:1690247/1‑150 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGc                                                        >  1:728453/1‑103 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGc                                                        <  2:728453/103‑1 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCtt         <  2:1439981/150‑1 (MQ=255)
                                                                                                                                                ccccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTGCGTCGCATTCTTGTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCGtc        >  1:171860/1‑150 (MQ=255)
                                                                                                                                                 cccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtct       <  2:1712074/150‑1 (MQ=255)
                                                                                                                                                 cccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtct       >  2:836281/1‑150 (MQ=255)
                                                                                                                                                  ccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGGTACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtc        <  2:622961/148‑1 (MQ=255)
                                                                                                                                                    aTAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCtt    <  2:2127879/150‑1 (MQ=255)
                                                                                                                                                     tAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTc   <  2:729936/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCC  >  NC_002947/5827869‑5828167

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: