Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,543 T→C 6.9% V136A (GTC→GCC)  PP_4611 → ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,233,5430TC 6.9% 204.0 / 4.7 72V136A (GTC→GCC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base C (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

CTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGT  >  NC_002947/5233397‑5233686
                                                                                                                                                       |                                                                                                                                               
cTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGCCg                                                                                                                                               >  2:115902/1‑149 (MQ=255)
   ccTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGAc                                                                                                                                             <  3:235065/148‑1 (MQ=255)
       cAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                          <  5:284214/147‑1 (MQ=255)
          gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTggg                                                                                                                                     <  3:90195/149‑1 (MQ=255)
             cTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCa                                                                                                                                  <  8:138633/149‑1 (MQ=255)
                gCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                               <  1:115902/149‑4 (MQ=255)
                 cacaGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGc                                                                                                                              >  8:28746/1‑149 (MQ=255)
                  acaGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGc                                                                                                                              >  6:258695/1‑148 (MQ=255)
                    aGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTCCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGa                                                                                                                            >  7:184967/1‑148 (MQ=255)
                      gTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAACCGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                          <  1:239133/149‑1 (MQ=255)
                       tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGAc                                                                                                                                             >  4:316019/1‑128 (MQ=255)
                       tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGAc                                                                                                                                             <  3:316019/128‑1 (MQ=255)
                       tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                        >  7:238363/1‑149 (MQ=255)
                         ccGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgc                                                                                                                      <  8:238363/149‑1 (MQ=255)
                                 ggTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGCACTCGCTCGCGAAc                                                                                                              >  2:245394/1‑149 (MQ=255)
                                        ccGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTc                                                                                                                                                <  1:217946/108‑1 (MQ=255)
                                        ccGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTc                                                                                                                                                >  2:217946/1‑108 (MQ=255)
                                         cGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAtt                                                                                                      <  1:121263/149‑1 (MQ=255)
                                          gTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAttt                                                                                                     <  6:295322/149‑1 (MQ=255)
                                              gAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGCCTAAAATCTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAttt                                                                                                     <  1:111138/149‑1 (MQ=255)
                                              gAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCt                                                                                                 <  8:173789/149‑1 (MQ=255)
                                                 gAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGc                                                                                              <  7:79103/149‑1 (MQ=255)
                                                  aaCAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGc                                                                                              >  5:262833/1‑148 (MQ=255)
                                                   aCAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCa                                                                                            >  1:179251/1‑149 (MQ=255)
                                                       aaGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                        >  5:301045/1‑149 (MQ=255)
                                                           cGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAg                                                                                    <  5:258695/149‑1 (MQ=255)
                                                            gATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgc                                                                                   >  3:217665/1‑149 (MQ=255)
                                                             atcatcTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgcg                                                                                  >  6:177136/1‑149 (MQ=255)
                                                              tcatcTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgcgc                                                                                 <  7:28746/149‑1 (MQ=255)
                                                               catcTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgcgc                                                                                 >  8:285932/1‑148 (MQ=255)
                                                                   ttCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                          >  6:103569/1‑103 (MQ=255)
                                                                   ttCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                          <  5:103569/103‑1 (MQ=255)
                                                                   ttCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCt                                                                            >  1:244256/1‑149 (MQ=255)
                                                                   ttCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCt                                                                            >  2:13076/1‑149 (MQ=255)
                                                                     cGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTCGCCACGGTCAAGCGCTACCTg                                                                           >  6:59253/1‑148 (MQ=255)
                                                                       agaCATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGaa                                                                         >  7:124182/1‑148 (MQ=255)
                                                                       agaCATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                        <  2:244256/149‑1 (MQ=255)
                                                                        gaCATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                        >  5:239988/1‑148 (MQ=255)
                                                                           aTTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTg                                                                                                            >  2:127731/1‑109 (MQ=255)
                                                                           aTTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTg                                                                                                            <  1:127731/109‑1 (MQ=255)
                                                                           aTTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                    <  6:301045/149‑1 (MQ=255)
                                                                                      tCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcgc                                                         >  3:50694/1‑149 (MQ=255)
                                                                                        gACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctg                                                       >  5:300304/1‑149 (MQ=255)
                                                                                         cccA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                      <  6:300304/148‑1 (MQ=255)
                                                                                         aCCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                      <  8:184967/149‑1 (MQ=255)
                                                                                           cA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTa                                                    >  1:194068/1‑149 (MQ=255)
                                                                                                       ggACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCGGCCATGCGCTGCTACTTCgcgc                                            >  4:115030/1‑149 (MQ=255)
                                                                                                        gACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCgcgc                                            <  7:285932/148‑1 (MQ=255)
                                                                                                           ggCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTa                                                    >  6:138857/1‑137 (MQ=255)
                                                                                                           ggCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTa                                                    <  5:138857/137‑1 (MQ=255)
                                                                                                             ccTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGaa                                      <  2:194068/149‑1 (MQ=255)
                                                                                                               ttCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACg                                    >  2:215488/1‑149 (MQ=255)
                                                                                                                 ccACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGc                                   <  8:124182/148‑1 (MQ=255)
                                                                                                                  cACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCt                                 >  5:281403/1‑149 (MQ=255)
                                                                                                                   aCCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                    <  8:237853/65‑1 (MQ=255)
                                                                                                                   aCCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                    >  7:237853/1‑65 (MQ=255)
                                                                                                                   aCCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                                >  5:152516/1‑149 (MQ=255)
                                                                                                                      gTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGtt                             >  7:228339/1‑149 (MQ=255)
                                                                                                                            cAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGcc                                                               >  1:132610/1‑109 (MQ=255)
                                                                                                                               a‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCt                    >  2:303147/1‑149 (MQ=255)
                                                                                                                                      cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCGATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc              <  5:59253/149‑1 (MQ=255)
                                                                                                                                         ccTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg             <  4:50694/147‑1 (MQ=255)
                                                                                                                                          gccttgGCCCAGGCCGTCGACCTGGGCCAGGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa          >  3:177210/5‑149 (MQ=255)
                                                                                                                                          ctgctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCGGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa          <  3:115030/149‑1 (MQ=255)
                                                                                                                                          ctgctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa          <  4:217665/149‑1 (MQ=255)
                                                                                                                                             ctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACa        <  8:228339/148‑1 (MQ=255)
                                                                                                                                                gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGcc    <  6:239988/149‑1 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt  <  2:179251/149‑1 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt  >  5:55225/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt  >  4:216939/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt  >  6:188724/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGAATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt  >  2:311941/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
CTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCA‑‑‑‑CTTGCTGGACGGCCTTCCACCGTTGGTCAAA‑CGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGT  >  NC_002947/5233397‑5233686

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