Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,828,015 C→T 5.9% intergenic (+161/+38) 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,0150CT 5.9% 468.7 / 12.8 153intergenic (+161/+38)mtgA/rpoHmonofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
Reads supporting (aligned to +/- strand):  ref base C (74/70);  new base T (3/6);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

Only 100 of 158 total aligned reads displayed.

CCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCT  >  NC_002947/5827868‑5828162
                                                                                                                                                   |                                                                                                                                                   
ccTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCTACCGCATCCCCTGTAGAACAGGGCCAAGGCGCTAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCAt                                                                                                                                                   >  1:2342885/1‑150 (MQ=255)
 cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                                  <  2:2167900/150‑1 (MQ=255)
 cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                                  <  2:282339/150‑1 (MQ=255)
 cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGAAGGAGCGGGCTTACCCGCAAAGAGGCCGGGAAAGCCAAACAAAAAAGACGCTCACCCCATa                                                                                                                                                  >  2:1898190/1‑150 (MQ=255)
 cTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCACTGCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                                  <  2:675077/150‑1 (MQ=255)
     cgcgGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGc                                                                                                                                              <  1:36483/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)
                ccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtttt                                                                                                                                   >  2:1655911/1‑150 (MQ=255)
                  gCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttttt                                                                                                                                 >  2:1668253/1‑148 (MQ=39)
                     ccTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAGAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgtt                                                                                                                              <  1:1921830/150‑6 (MQ=37)
                     ccTACGGGCAACCGAATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                              <  1:288035/150‑1 (MQ=255)
                      cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAGAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                             <  1:1517162/150‑7 (MQ=37)
                        aCGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACaa                                                                                                                           >  2:11468/1‑150 (MQ=255)
                          gggCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAcc                                                                                                                         >  2:1020773/1‑150 (MQ=255)
                            gCAACCGCATCCCCGGCCTCCTCGCGGGCCAAACACCGGCTGGCCATAACCTCAGCCCCTGTAGGAGCGGGAAACCTCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCa                                                                                                                       <  1:1739473/150‑1 (MQ=255)
                                ccGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGTCGACCTCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACt                                                                                                                   >  2:2269360/1‑150 (MQ=255)
                                    gtCCGCGGCGTCCGCGCGGGCGAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAAt                                                                                                               <  2:2188120/149‑1 (MQ=255)
                                     tCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                              <  2:2344992/150‑1 (MQ=255)
                                        ccGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCAGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAgg                                                                                                           >  1:758360/1‑150 (MQ=255)
                                         cGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGc                                                                                                          >  1:481378/1‑150 (MQ=255)
                                             ctcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTc                                                                                                                                  <  2:945388/122‑1 (MQ=255)
                                              tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAgg                                                                                                           <  1:29562/144‑1 (MQ=255)
                                              tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                     >  1:2058368/1‑150 (MQ=255)
                                                 tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATa                                                                                                                                                  >  1:945998/1‑102 (MQ=255)
                                                 tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCaa                                                                                                                >  1:163958/1‑136 (MQ=255)
                                                 tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCa                                                                                                  <  2:758360/150‑1 (MQ=255)
                                                    cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                               <  1:1015968/150‑1 (MQ=255)
                                                    cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                               <  1:1020773/150‑1 (MQ=255)
                                                          aaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAgg                                                                                                                                                >  2:2237229/1‑95 (MQ=255)
                                                             ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                              >  2:1452693/1‑110 (MQ=255)
                                                             ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                              <  1:1452693/110‑1 (MQ=255)
                                                             ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                      >  1:1507857/1‑150 (MQ=255)
                                                                cgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                   <  1:1001304/150‑1 (MQ=255)
                                                                     aCAGGCAACCGCATCCGCTGTAGGAGCGTGCTTGCCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                              <  1:2054126/150‑1 (MQ=255)
                                                                      cAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGACGAGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                             >  2:2131660/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)
                                                                              cgaaTCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAGGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                                     >  1:398706/4‑150 (MQ=255)
                                                                               gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCTGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                      >  2:119048/1‑132 (MQ=255)
                                                                               gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATtcttct                                                                    <  2:2380858/150‑1 (MQ=255)
                                                                                cATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                                     >  2:2037523/1‑148 (MQ=255)
                                                                                    cccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAAc                                                               <  1:11468/150‑1 (MQ=255)
                                                                                    cccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAAc                                                               <  2:641149/150‑1 (MQ=255)
                                                                                     ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                              <  1:1389048/86‑1 (MQ=255)
                                                                                     ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                              >  2:1389048/1‑86 (MQ=255)
                                                                                     ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                                   >  1:69047/1‑145 (MQ=255)
                                                                                     ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACt                                                              >  1:319262/1‑150 (MQ=255)
                                                                                      cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                             >  2:107488/1‑150 (MQ=255)
                                                                                      cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                             >  2:2276985/1‑150 (MQ=255)
                                                                                        gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                           <  1:1656054/150‑1 (MQ=255)
                                                                                            gAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAAt                                                       >  1:1067328/1‑150 (MQ=255)
                                                                                                 ggCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                  >  1:1661503/1‑150 (MQ=255)
                                                                                                  gtttACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                  <  2:404833/147‑1 (MQ=255)
                                                                                                  gtttACCCGCGAAAGGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                  >  1:404833/3‑149 (MQ=255)
                                                                                                        ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGaa                                                        >  2:1631884/1‑137 (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)
                                                                                                                    cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTAc                                 >  2:712170/1‑148 (MQ=255)
                                                                                                                     ggCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTg                              >  2:1296866/1‑150 (MQ=255)
                                                                                                                      gCAAAGCCAAACAAAAAAGCCTCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGACTACGCCCAGCCGCTACGCTGTACTTGt                             <  2:1240782/150‑1 (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)
                                                                                                                           gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                         >  1:1617152/1‑85 (MQ=255)
                                                                                                                           gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                         <  2:1617152/85‑1 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                   <  1:2206258/90‑1 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                   >  2:2206258/1‑90 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                       >  1:2399558/1‑150 (MQ=255)
                                                                                                                            ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                       >  1:978828/1‑150 (MQ=255)
                                                                                                                             caaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                               <  1:2266073/141‑1 (MQ=255)
                                                                                                                              aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                           >  2:574048/1‑112 (MQ=255)
                                                                                                                              aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                           <  1:574048/112‑1 (MQ=255)
                                                                                                                              aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt                     >  1:2127879/1‑150 (MQ=255)
                                                                                                                              aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt                     >  2:177312/1‑150 (MQ=255)
                                                                                                                                     aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTCCCACCCAGACTCAATCAGCCAGCGACAAGTGCCTTAACCTTGTTCATCGCATTCTTCTCCAACTGCCGAATAAGCTCCGCCGCTACGCTCTACTTGTATGCCAGTTCATGcc              >  2:680052/1‑149 (MQ=255)
                                                                                                                                     aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCa              >  1:1927530/1‑150 (MQ=255)
                                                                                                                                         ccGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTAATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGt          <  2:2058368/150‑1 (MQ=255)
                                                                                                                                           gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATa                                         <  1:2235425/117‑1 (MQ=255)
                                                                                                                                           gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTgg        >  1:1690247/1‑150 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGt                                                                                                <  2:2237929/58‑1 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGt                                                                                                >  1:2237929/1‑58 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGc                                                   <  2:728453/103‑1 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCtt    >  2:2134604/1‑150 (MQ=255)
                                                                                                                                               aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCtt    <  2:1555475/150‑1 (MQ=255)
                                                                                                                                                ccccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTGCGTCGCATTCTTGTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCGtc   >  1:171860/1‑150 (MQ=255)
                                                                                                                                                 cccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtct  >  2:836281/1‑150 (MQ=255)
                                                                                                                                                  ccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGGTACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtc   <  2:622961/148‑1 (MQ=255)
                                                                                                                                                  ccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGGTACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtc   >  1:622961/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCT  >  NC_002947/5827868‑5828162

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