Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,543 T→C 6.4% 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.4% 229.6 / 4.3 78V136A (GTC→GCC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (2/3);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCAT  >  NC_002947/5233397‑5233690
                                                                                                                                                       |                                                                                                                                                   
cTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGGTCAAACGCGCCTTCCTGCTGGCCCAGGTc                                                                                                                                                    <  1:226571/149‑1 (MQ=255)
       cAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCt                                                                                                                                            >  6:111286/1‑149 (MQ=255)
        aaGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTg                                                                                                                                           >  2:134036/1‑149 (MQ=255)
          gcgcTCGCACAGGTCCCTGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGTTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTggg                                                                                                                                         <  6:300386/149‑1 (MQ=255)
          gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                          <  1:134036/148‑1 (MQ=255)
          gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGACTGGGCCACGGTTGGTCAAACGCGCCTACCTGCTGGCCCAGGTCGCCGGCCTg                                                                                                                                           >  1:312549/1‑149 (MQ=255)
           cgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGc                                                                                                                                        <  4:10057/149‑1 (MQ=255)
              tCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAAGAGGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGc                                                                                                                                        >  1:307456/1‑148 (MQ=255)
              tCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGcc                                                                                                                                       >  6:321754/1‑147 (MQ=255)
                gCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAgg                                                                                                                                   >  3:240860/1‑149 (MQ=255)
                acacaGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                   <  4:139504/148‑4 (MQ=255)
                 cacaGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACTGCCTGGGCCAAGGc                                                                                                                                  >  1:361866/1‑149 (MQ=255)
                     ggTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                              <  5:158978/149‑1 (MQ=255)
                      gTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                    >  4:143663/1‑142 (MQ=255)
                      gTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                    <  3:143663/142‑1 (MQ=255)
                         ccGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgc                                                                                                                          >  5:215813/1‑149 (MQ=255)
                           ggAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                        <  6:298641/149‑1 (MQ=255)
                            gAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctcg                                                                                                                       >  4:154689/1‑149 (MQ=255)
                              agagGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTcgcg                                                                                                                     <  5:321754/149‑1 (MQ=255)
                                     gcgcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                               <  6:215813/148‑1 (MQ=255)
                                        cggTGCGAGGAACAGAAGGCTATCATATTCGAGACATTGGTAGCGCTGGACCAGTTGCTGGACGGCGT‑‑TCCACCGTTGGTAAAACGCGCCTTCCTGCTGGCCAAGGTCTACGGCGTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAt                                                                                                           <  2:361866/147‑1 (MQ=255)
                                        ccGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                   <  7:197943/125‑4 (MQ=255)
                                        ccGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                   >  8:197943/1‑122 (MQ=255)
                                              gAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCt                                                                                                     <  3:154689/149‑1 (MQ=255)
                                                  aaCAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGcc                                                                                                 <  8:350896/149‑1 (MQ=255)
                                                    cAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCAc                                                                                               <  2:307456/149‑1 (MQ=255)
                                                    cAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGc                                                                                                                                        <  5:386409/108‑1 (MQ=255)
                                                    cAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGc                                                                                                                                        >  6:386409/1‑108 (MQ=255)
                                                             atcatcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgcg                                                                                      >  4:71247/1‑149 (MQ=255)
                                                              tcatcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgcgc                                                                                     <  3:192186/149‑1 (MQ=255)
                                                                  cTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAcc                                                                                 >  4:94570/1‑149 (MQ=255)
                                                                    tCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTg                                                                                                                >  6:116999/1‑116 (MQ=255)
                                                                    tCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTg                                                                                                                <  5:116999/116‑1 (MQ=255)
                                                                     cGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCt                                                                                <  4:73454/147‑1 (MQ=255)
                                                                       agaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGaa                                                                             >  1:63084/1‑148 (MQ=255)
                                                                        gaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACgg                                                                                             <  5:2754/131‑1 (MQ=255)
                                                                        gaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACgg                                                                                             >  6:2754/1‑131 (MQ=255)
                                                                            ttGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGc                                                                       <  5:111286/149‑1 (MQ=255)
                                                                             tGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                        >  1:37651/1‑147 (MQ=255)
                                                                               cTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAg                                                                     <  2:312549/148‑1 (MQ=255)
                                                                               cTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAg                                                                     <  2:37651/148‑1 (MQ=255)
                                                                               cTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGc                                                                    >  4:189851/1‑149 (MQ=255)
                                                                                tAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                        <  4:292909/96‑1 (MQ=255)
                                                                                tAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctc                                                                                                                        >  3:292909/1‑96 (MQ=255)
                                                                                              aCTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGc                                                                                                                                  >  8:165694/1‑74 (MQ=255)
                                                                                              aCTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGc                                                                                                                                  <  7:165694/74‑1 (MQ=255)
                                                                                              aCTTGCTGGACGGCCT‑‑CCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTAc                                                       >  1:260057/1‑149 (MQ=255)
                                                                                                ttGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                     >  6:171097/1‑149 (MQ=255)
                                                                                                ttGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                     >  5:57247/1‑149 (MQ=255)
                                                                                                     ggACGGCCT‑‑TCCACAGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCgcgc                                                >  2:184917/1‑149 (MQ=255)
                                                                                                      gACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCt                                               <  5:338862/149‑1 (MQ=255)
                                                                                                         ggCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                             >  3:296421/1‑148 (MQ=255)
                                                                                                           ccT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGaa                                          <  3:94570/149‑1 (MQ=255)
                                                                                                             t‑‑tCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACg                                        >  4:108775/1‑149 (MQ=255)
                                                                                                                tCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGc                                       >  7:86665/1‑149 (MQ=255)
                                                                                                                 ccACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGcc                                      >  6:129733/1‑149 (MQ=255)
                                                                                                                  cACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCt                                     <  1:184917/149‑1 (MQ=255)
                                                                                                                   aCCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                                    <  3:71247/149‑1 (MQ=255)
                                                                                                                   aCCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                                    <  4:240860/149‑1 (MQ=255)
                                                                                                                      gTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGtt                                 <  4:296421/149‑1 (MQ=255)
                                                                                                                           gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTc                              <  3:108775/148‑1 (MQ=255)
                                                                                                                           gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTc                              <  2:260057/148‑1 (MQ=255)
                                                                                                                           gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTc                              <  5:129733/148‑1 (MQ=255)
                                                                                                                              aaaCGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCg                          >  1:351534/1‑149 (MQ=255)
                                                                                                                                 cgcgCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtg                       <  6:57247/149‑1 (MQ=255)
                                                                                                                                    gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACACGGCCGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggtg                    >  3:246010/1‑149 (MQ=255)
                                                                                                                                    gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGTCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggtg                    <  2:63084/149‑1 (MQ=255)
                                                                                                                                          ctgctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa              <  5:171097/149‑1 (MQ=255)
                                                                                                                                            gctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgc                                                               >  5:41268/1‑98 (MQ=255)
                                                                                                                                            gctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgc                                                               <  6:41268/98‑1 (MQ=255)
                                                                                                                                                gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGcc        >  1:95037/1‑149 (MQ=255)
                                                                                                                                                gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACCGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGcc        >  7:247786/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt      >  4:183026/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt      >  2:141789/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt      <  3:189851/149‑1 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt      >  4:91749/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGAGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt      >  5:72701/1‑149 (MQ=255)
                                                                                                                                                      gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCAt  >  8:209476/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CTACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCT‑‑TCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCAT  >  NC_002947/5233397‑5233690

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: