Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,827,966 C→T 5.0% intergenic (+168/+31) 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,827,9660CT 5.0% 310.4 / 3.5 99intergenic (+168/+31)mtgA/rpoHmonofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
Reads supporting (aligned to +/- strand):  ref base C (41/53);  new base T (3/2);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCC  >  NC_002947/5827822‑5828111
                                                                                                                                                |                                                                                                                                                  
gTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                                <  6:73515/149‑1 (MQ=255)
    aCCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCt                                                                                                                                            >  3:365982/1‑149 (MQ=255)
     cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtt                                                                                                                                           <  5:358408/149‑1 (MQ=255)
         cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCt                                                                                                                                       <  5:336129/149‑1 (MQ=255)
            caaCGGGCAAACGCATCCCCGGCCTCCTCGCGGGCAAACCCGCCCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                                      <  4:237609/145‑1 (MQ=255)
            cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCttt                                                                                                                                     <  2:364669/148‑1 (MQ=255)
            cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCttt                                                                                                                                     >  3:237609/1‑148 (MQ=255)
            cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                                    >  2:60533/1‑149 (MQ=255)
             tACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCCCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                                    >  6:149844/1‑148 (MQ=255)
               cGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACaa                                                                                                                                 <  4:365982/149‑1 (MQ=255)
                   cAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCa                                                                                                                             <  5:218819/149‑1 (MQ=255)
                    aaCCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                                   >  7:263629/1‑136 (MQ=37)
                    aaCCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                                   <  8:263629/142‑7 (MQ=37)
                    aaCCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAcc                                                                                                                               >  5:200055/1‑146 (MQ=255)
                    aaCCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAcc                                                                                                                               >  6:408946/1‑146 (MQ=255)
                     ctatcttCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCTGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGa                                                                                                                           <  3:24778/143‑1 (MQ=255)
                      ccGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGTCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGAc                                                                                                                          <  6:220831/149‑1 (MQ=255)
                      ccGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGAc                                                                                                                          <  3:333139/149‑1 (MQ=255)
                        gCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTCCCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCCTAGGGCCAGCGGCGTTTCTGTACAAccccgcct                                                                                                                         >  3:245750/1‑143 (MQ=39)
                            ccccGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AAt                                                                                                                    >  1:119219/1‑148 (MQ=255)
                            ccccGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTAACAAACCaca                                                                                                                           >  4:24778/1‑140 (MQ=255)
                             cccGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AAt                                                                                                                    >  8:212293/1‑147 (MQ=255)
                             cccGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AAt                                                                                                                    >  7:77149/1‑147 (MQ=255)
                             cccGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCa                                                                                                                  <  4:235222/149‑1 (MQ=255)
                                    tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttt                                                                                                                                        >  8:22377/1‑120 (MQ=38)
                                    tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttt                                                                                                                                        <  7:22377/121‑2 (MQ=37)
                                      ctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGAt                                                                                                         >  2:51338/1‑149 (MQ=255)
                                       tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGAt                                                                                                         >  6:283081/1‑148 (MQ=255)
                                       tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGAt                                                                                                         >  1:263520/1‑148 (MQ=255)
                                       tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGAt                                                                                                         >  1:81753/1‑148 (MQ=255)
                                       tcGCGGACAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGAt                                                                                                         >  6:223401/1‑148 (MQ=255)
                                         gcgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCa                                                                                                       <  1:60533/148‑1 (MQ=255)
                                               aaaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAg                                                                                                                 >  1:412119/1‑132 (MQ=255)
                                               aaaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAg                                                                                                                 <  2:412119/132‑1 (MQ=255)
                                               aaaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCa                                                                                                              <  4:229517/135‑1 (MQ=255)
                                               aaaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCa                                                                                                              >  3:229517/1‑135 (MQ=255)
                                               aaaCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCtt                                                                                                <  7:23420/149‑1 (MQ=255)
                                                 aCCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGa                                                                                              <  4:179876/149‑1 (MQ=255)
                                                    cgcgCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCACTGATCGCTGCCtgattga                                                                                           <  7:311799/149‑7 (MQ=39)
                                                     gcgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCtt                                                                                          <  5:223401/149‑1 (MQ=255)
                                                      cgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTg                                                                                         <  7:212293/149‑1 (MQ=255)
                                                         cTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCACCGATCAGTGCCTTGACCTTGtt                                                                                       >  7:293197/1‑148 (MQ=255)
                                                            cAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCAt                                                                                    <  8:77149/148‑1 (MQ=255)
                                                            cAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                                   <  5:149844/149‑1 (MQ=255)
                                                               gCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGc                                                                                 >  3:15302/1‑148 (MQ=255)
                                                               gCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGc                                                                                 <  4:15302/148‑1 (MQ=255)
                                                                cAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAt                                                                               <  5:208467/149‑1 (MQ=255)
                                                                 aaCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAccga                                                                                                                             >  6:42662/1‑101 (MQ=255)
                                                                 aaCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAccga                                                                                                                             <  5:42662/103‑3 (MQ=255)
                                                                  aCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttc                                                                             <  6:200055/149‑1 (MQ=255)
                                                                    cGCATCCGCTGTAGGAGCGGTCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATATGGCCAGAGGCTTTTCTTTACAACCCAGACTC‑AATCAGGGAGCGATCAGTGCCTTGACCTTGTTCATCGCAttctt                                                                           <  5:283081/149‑1 (MQ=255)
                                                                    cGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttctt                                                                           <  6:317266/149‑1 (MQ=255)
                                                                     gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGc                                                                                                               >  6:350853/1‑112 (MQ=255)
                                                                      cATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGtt                                                                                       >  4:58070/1‑135 (MQ=255)
                                                                      cATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGtt                                                                                       <  3:58070/135‑1 (MQ=255)
                                                                      cATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATtcttct                                                                         >  4:286464/1‑149 (MQ=255)
                                                                              gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTgg                                                                 >  8:364906/1‑149 (MQ=255)
                                                                              gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTgg                                                                 >  3:132807/1‑149 (MQ=255)
                                                                               tAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAATAAGCCGCTCACCCCATTGGTCCAGCTGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGGCTTGACCTTGTTCATCGCATTCTTCTCCAACTgg                                                                 <  6:358453/148‑1 (MQ=255)
                                                                               tAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                <  7:364906/149‑1 (MQ=255)
                                                                                aGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                >  1:204118/1‑148 (MQ=255)
                                                                                 ggAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCg                                                               <  5:408946/148‑1 (MQ=255)
                                                                                       ggCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGc                                                        <  2:263520/149‑1 (MQ=255)
                                                                                              ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCa                                                                                                       >  7:218976/1‑95 (MQ=255)
                                                                                              ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCa                                                                                                       <  8:218976/95‑1 (MQ=255)
                                                                                              ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGt                                                                                                     <  6:179137/97‑1 (MQ=255)
                                                                                              ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGt                                                                                                     >  5:179137/1‑97 (MQ=255)
                                                                                              ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGcc                                                  >  1:185240/1‑148 (MQ=255)
                                                                                              ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGcc                                                  >  6:211644/1‑148 (MQ=255)
                                                                                               cgcgAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGAcc                                                                                            <  1:284778/105‑1 (MQ=255)
                                                                                               cgcgAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGAcc                                                                                            >  2:284778/1‑105 (MQ=255)
                                                                                                   aaGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACg                                            <  3:65251/149‑1 (MQ=255)
                                                                                                     gagGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCt                                          <  1:51338/149‑1 (MQ=255)
                                                                                                      agGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTg                                         <  2:81753/149‑1 (MQ=255)
                                                                                                      agGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTg                                         <  5:98392/149‑1 (MQ=255)
                                                                                                        gCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTa                                       >  5:283280/1‑149 (MQ=255)
                                                                                                         ccGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTA‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTAc                                      <  2:185240/149‑1 (MQ=255)
                                                                                                              aaaGCCAAACAAAAACGCCGCTCACCCCATACGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTc                                 <  4:245750/149‑1 (MQ=255)
                                                                                                                 gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGc                              <  2:119219/149‑1 (MQ=255)
                                                                                                                 gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGc                              <  2:204118/149‑1 (MQ=255)
                                                                                                                     aacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGCATACGCTCAGCCGATACGCTGTACTGGTCTGCCAGt                          >  3:330623/1‑149 (MQ=255)
                                                                                                                     aacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGt                          <  1:249746/149‑1 (MQ=255)
                                                                                                                        aaaaacTCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCa                       <  6:167543/149‑1 (MQ=255)
                                                                                                                         aaaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCAt                      <  3:286464/149‑1 (MQ=255)
                                                                                                                          aaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg                     >  2:233209/1‑149 (MQ=255)
                                                                                                                            aaGCCGCTGGCCCTATGGGGTGAGCGGCTTTTCTTTACAACCCAGACTCTAATCAGGCAGCGATCAGTGCCTTGACCCTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg                     >  1:95381/1‑148 (MQ=255)
                                                                                                                               ccGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                                  >  4:232690/1‑83 (MQ=255)
                                                                                                                               ccGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                                  <  3:232690/83‑1 (MQ=255)
                                                                                                                                  cTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTg                                                                                               <  7:274938/67‑1 (MQ=255)
                                                                                                                                  cTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTg                                                                                               >  8:274938/1‑67 (MQ=255)
                                                                                                                                   gcACCCCATCGGGCCAGCGGCTTTTCTTTACACCCCAGACTC‑ACTCAGGCAGCGATCAGTGCCTTGACCTTGTGCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGc            <  4:330623/148‑1 (MQ=255)
                                                                                                                                     aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                                   >  8:12805/1‑76 (MQ=255)
                                                                                                                                     aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                                   <  7:12805/76‑1 (MQ=255)
                                                                                                                                     aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGcc           <  5:211644/148‑1 (MQ=255)
                                                                                                                                       cccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGCTCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGACGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCTTGGCCTtc        >  8:186365/1‑149 (MQ=255)
                                                                                                                                         cATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtctc      <  6:283280/149‑1 (MQ=255)
                                                                                                                                          aTAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTctct     >  6:423029/1‑149 (MQ=255)
                                                                                                                                              ggCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGTCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTcc  <  7:147645/148‑1 (MQ=255)
                                                                                                                                              ggCCAGCGGCTTTTCTTTAAAACACAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTCCCGAATACGCTCAGCCAATAACCAGTACTTGTCTGCAAGTTCATTCAGCGTGGCCTTCTCTTcc  >  8:147645/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTC‑AATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCC  >  NC_002947/5827822‑5828111

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