Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,594 A→T 7.5% D137V (GAC→GTC)  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,5940AT 7.5% 231.6 / 7.8 79D137V (GAC→GTC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (45/28);  new base T (4/2);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

ACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTT  >  NC_002947/5233447‑5233743
                                                                                                                                                   |                                                                                                                                                     
cccTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACg                                                                                                                                                     <  2:432974/149‑1 (MQ=255)
  cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                                   >  1:1598340/1‑150 (MQ=255)
  cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                                   >  1:1899300/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                                                                                                                                   >  2:1284755/1‑150 (MQ=255)
                  aGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                                   >  1:944428/1‑150 (MQ=255)
                     tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                         >  1:1864660/1‑141 (MQ=255)
                     tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                         <  2:1864660/141‑1 (MQ=255)
                        cGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                                 >  1:1793299/1‑146 (MQ=255)
                        cGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                                 <  2:1793299/146‑1 (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                                                                                                                                         >  1:637347/1‑99 (MQ=255)
                                                               tcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                         <  2:637347/99‑1 (MQ=255)
                                                                 ttCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                                    >  1:1832576/1‑118 (MQ=255)
                                                                 ttCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTgg                                                                                                                    <  2:1832576/118‑1 (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                                                                                                                  >  1:1899076/1‑103 (MQ=255)
                                                                                  gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGc                                                                                                                  <  2:1899076/103‑1 (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                                                  <  1:294571/114‑1 (MQ=255)
                                                                                                                                     cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                                  >  2:294571/3‑116 (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)
                                                                                                                                                   |                                                                                                                                                     
ACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGTT  >  NC_002947/5233447‑5233743

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