Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,597 G→A 6.3% G138D (GGC→GAC)  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,5970GA 6.3% 222.6 / 4.5 79G138D (GGC→GAC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base G (47/27);  new base A (3/2);  total (50/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTTGA  >  NC_002947/5233449‑5233745
                                                                                                                                                    |                                                                                                                                                    
cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                                     >  1:1899300/1‑150 (MQ=255)
cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                                     >  1:1598340/1‑150 (MQ=255)
 tGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGcc                                                                                                                                                    <  1:3408/150‑1 (MQ=255)
    aaGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                                 <  2:1913863/150‑1 (MQ=255)
     ggcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCGCCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGGCGACGGCCTggg                                                                                                                                                <  2:297334/149‑1 (MQ=255)
     aGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTggg                                                                                                                                                >  1:1141992/1‑150 (MQ=255)
         cTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCaa                                                                                                                                            <  2:1738942/150‑1 (MQ=255)
           cGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                          <  2:1396566/150‑4 (MQ=255)
            gCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCACGGc                                                                                                                                         >  2:290324/1‑150 (MQ=255)
            gCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACGTGCTGGCCAGCCACAAACTCACAGTTCGTCTGGCCTTTCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGc                                                                                                                                         >  1:889751/1‑150 (MQ=255)
                aGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                                     >  1:944428/1‑150 (MQ=255)
                aGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                                     >  2:1284755/1‑150 (MQ=255)
                   tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                           <  2:1864660/141‑1 (MQ=255)
                   tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                           >  1:1864660/1‑141 (MQ=255)
                      cGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                                   <  2:1793299/146‑1 (MQ=255)
                      cGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                                   >  1:1793299/1‑146 (MQ=255)
                       ggccgagGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGACAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctcg                                                                                                                              >  2:503306/5‑150 (MQ=255)
                              gTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTg                                                                                                                                                  >  2:1601290/1‑123 (MQ=255)
                              gTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTg                                                                                                                                                  <  1:1601290/123‑1 (MQ=255)
                               tCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                                      <  1:1298993/150‑1 (MQ=255)
                                   gcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAt                                                                                                                  >  2:1769260/1‑150 (MQ=255)
                                                 aGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCCCTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTGGCGAACTGGGCATTTCGCTGGCCACgg                                                                                                    <  1:990353/150‑1 (MQ=255)
                                                  gAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCCCGAACTGGGCATTTCGCTGGCCACGGt                                                                                                   >  2:54192/1‑150 (MQ=255)
                                                   aaGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTc                                                                                                  >  1:784541/1‑150 (MQ=255)
                                                        gATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTGCACCGTCGGCCAAACACGCCGTCCGGGTGGCCCGGGGGGACAGCCTGGGCCAAGCCGAAATCGCTCGCGAGGTGGGCGTTTCGCTGACCACGGTCAAgcg                                                                                             >  2:1883069/1‑150 (MQ=255)
                                                          tcatcTTCGAGACATTGCTAGCGCTCGAGCACTTGCTGGACGGCCTTCCGCCGTTGGACACAAGCGCCTTAGTGACGGCGGAGGTTATTTGCAGGGCCCGAGGCGAAATCGCTCTCGAACTGGGCATTTAGCTGGCGACGGTTGAGCGCt                                                                                           >  1:896897/1‑150 (MQ=255)
                                                             tcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCATGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAcc                                                                                        >  2:279861/1‑150 (MQ=255)
                                                             tcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                           <  2:637347/99‑1 (MQ=255)
                                                             tcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                           >  1:637347/1‑99 (MQ=255)
                                                               ttCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                                      <  2:1832576/118‑1 (MQ=255)
                                                               ttCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                                      >  1:1832576/1‑118 (MQ=255)
                                                                   agaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATAGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACa                                                                                  >  1:395279/1‑150 (MQ=255)
                                                                    gaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACaa                                                                                 >  2:976353/1‑150 (MQ=255)
                                                                     aCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCCCGCGAACTGGGCATTTCGCTGGCCACGATCAATCGCTACCTGAACAAg                                                                                <  1:65007/150‑1 (MQ=255)
                                                                      cATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                               <  2:667679/150‑1 (MQ=255)
                                                                      cATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                               <  2:651117/150‑1 (MQ=255)
                                                                            tAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCa                                                                         >  2:1230326/1‑150 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGc                                                                                                                    <  2:1899076/103‑1 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGc                                                                                                                    >  1:1899076/1‑103 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGc                                                                                                             >  1:1904744/1‑110 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGc                                                                                                             <  2:1904744/110‑1 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCCTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcg                                                                     >  2:942774/1‑150 (MQ=255)
                                                                                gatCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcg                                                                     >  1:749296/3‑150 (MQ=255)
                                                                                  tCGACAACATGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgct                                                                   >  1:1113011/1‑150 (MQ=255)
                                                                                    gACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                                 >  2:1067075/1‑150 (MQ=255)
                                                                                    gACCACGGGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                                 >  1:1418461/1‑150 (MQ=255)
                                                                                      ccACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAcc                                                                                        >  2:1069151/1‑125 (MQ=255)
                                                                                      ccACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAcc                                                                                        <  1:1069151/125‑1 (MQ=255)
                                                                                         cTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCCTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                            <  1:942774/150‑1 (MQ=255)
                                                                                          ttGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCGGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTc                                                           <  2:819064/150‑1 (MQ=255)
                                                                                             cTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcgcg                                                        <  1:1067075/150‑1 (MQ=255)
                                                                                              tGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGATACTTCgcgc                                                       >  1:391924/1‑150 (MQ=255)
                                                                                                 aCGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                                    <  2:1899300/150‑1 (MQ=255)
                                                                                                      cTTNCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGTTCATTTAGCTGGCCAAGGGCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACg                                               >  2:1204697/1‑150 (MQ=255)
                                                                                                       ttCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTg                                                                                                           <  2:1777360/89‑1 (MQ=255)
                                                                                                       ttCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTg                                                                                                           >  1:1777360/1‑89 (MQ=255)
                                                                                                        tCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGcc                                             >  1:580288/1‑150 (MQ=255)
                                                                                                                   tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAAc                                                                                                                         >  2:1130811/1‑63 (MQ=255)
                                                                                                                   tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAAc                                                                                                                         <  1:1130811/63‑1 (MQ=255)
                                                                                                                   tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcccc                                  >  1:1028164/1‑150 (MQ=255)
                                                                                                                   tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcccc                                  >  2:72783/1‑150 (MQ=255)
                                                                                                                     aaaCGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGc                                >  1:41826/1‑150 (MQ=255)
                                                                                                                       aCGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTAGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtg                              >  1:808968/1‑150 (MQ=255)
                                                                                                                         gcgcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggt                            >  1:335596/1‑150 (MQ=255)
                                                                                                                         ccgcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggt                            <  1:316871/149‑1 (MQ=255)
                                                                                                                             cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg                        >  2:640429/1‑150 (MQ=255)
                                                                                                                             cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg                        >  2:37566/1‑150 (MQ=255)
                                                                                                                              ttCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCCAACACTACCGGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCgg                       >  2:1097236/1‑150 (MQ=255)
                                                                                                                              ttCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCgg                       >  1:545234/1‑150 (MQ=255)
                                                                                                                                  tgctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTa                                                                                          <  1:410226/79‑1 (MQ=255)
                                                                                                                                  tgctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTa                                                                                          >  2:410226/1‑79 (MQ=255)
                                                                                                                                   cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                                    >  2:294571/3‑116 (MQ=255)
                                                                                                                                   cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                                    <  1:294571/114‑1 (MQ=255)
                                                                                                                                   cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAg                  >  2:1514435/3‑150 (MQ=255)
                                                                                                                                       gCCCAGGCCGTCGACCTGGGCCAAGGCGAAATGGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCg              >  2:1068435/1‑150 (MQ=255)
                                                                                                                                        cccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt             >  1:952027/1‑150 (MQ=255)
                                                                                                                                          cgggTTGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCTACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGc           <  2:306894/148‑1 (MQ=255)
                                                                                                                                             gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGAGACGGCGGGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCAtt        <  2:1113011/150‑1 (MQ=255)
                                                                                                                                                 aCGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGtt    >  1:315864/1‑150 (MQ=255)
                                                                                                                                                   ggCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAATGCAGCAATTCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTTGa  <  2:395279/150‑1 (MQ=255)
                                                                                                                                                   ggCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTTGa  <  2:1598340/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTTGA  >  NC_002947/5233449‑5233745

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