Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,546 A→T 9.9% 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,5460AT 9.9% 228.7 / 12.3 89D137V (GAC→GTC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base T (5/4);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

ACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATT  >  NC_002947/5233399‑5233691
                                                                                                                                                    |                                                                                                                                                 
aCCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGAc                                                                                                                                                  >  4:343647/1‑149 (MQ=255)
 ccTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACg                                                                                                                                                 <  6:36880/149‑1 (MQ=255)
   tGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCAAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                               >  5:243730/1‑149 (MQ=255)
    gCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGcc                                                                                                                                              >  6:123042/1‑149 (MQ=255)
       aGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                           >  2:388058/1‑149 (MQ=255)
       aGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                           <  5:337224/149‑1 (MQ=255)
        gcgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTggg                                                                                                                                          <  6:103058/149‑1 (MQ=255)
          gcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGc                                                                                                                                         >  5:120769/1‑148 (MQ=255)
          gcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCATGCTGGACCAGGTCGACGGCCTGGGc                                                                                                                                         >  6:379065/1‑148 (MQ=255)
           cTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGATCAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCa                                                                                                                                       <  1:138039/149‑1 (MQ=255)
            tCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTg                                                                                                                                            <  8:222920/143‑1 (MQ=255)
            tCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTg                                                                                                                                            >  7:222920/1‑143 (MQ=255)
               cacaGGGCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGc                                                                                                                                   <  3:307474/149‑1 (MQ=255)
                     tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                             >  1:239103/1‑149 (MQ=255)
                     tCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGACCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtc                                                                                                                             >  1:123404/1‑149 (MQ=255)
                      cccGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcg                                                                                                                            >  4:26598/1‑149 (MQ=255)
                           aaGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctcgc                                                                                                                       <  5:219157/149‑1 (MQ=255)
                                gTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCag                                                                                                                                      >  5:310371/1‑128 (MQ=255)
                                gTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCag                                                                                                                                      <  6:310371/129‑2 (MQ=255)
                                   gcgcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGGGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTggg                                                                                                               >  1:85762/1‑149 (MQ=255)
                                      ccGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAt                                                                                                            <  3:167929/149‑1 (MQ=255)
                                        gTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAtt                                                                                                           <  1:280777/149‑1 (MQ=255)
                                        gTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAttt                                                                                                          <  2:222476/149‑1 (MQ=255)
                                        gTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGGCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAttt                                                                                                          <  3:25122/149‑1 (MQ=255)
                                            gAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCt                                                                                                      >  4:305759/1‑149 (MQ=255)
                                              ggAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTTGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTgg                                                                                                    <  5:379065/149‑1 (MQ=255)
                                                   aGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCa                                                                                                             <  2:286717/135‑1 (MQ=255)
                                                   aGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACg                                                                                               <  2:239103/149‑1 (MQ=255)
                                                     aaGGCTACCTTCTTCCAGACATTGCTAGCGCCCGACCACTTGCTGGACGGCCTTCCCCCGATGGACAAACGCGCCTTCCTGCTGGCCCAGGTCGCCGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                             <  2:197338/149‑1 (MQ=255)
                                                          gATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgc                                                                                        >  2:248678/1‑149 (MQ=255)
                                                          gATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgc                                                                                        >  7:290648/1‑149 (MQ=255)
                                                              atcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCt                                                                                     <  2:97695/148‑1 (MQ=255)
                                                              atcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCt                                                                                     <  4:154883/148‑1 (MQ=255)
                                                                cTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAcc                                                                                  >  3:164605/1‑149 (MQ=255)
                                                                  tCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCt                                                                                 >  2:27798/1‑148 (MQ=255)
                                                                   cGATACATTGCTAGCGCTCTACCACTTGCTGGACGGCCTGCCACCTTTGGTCAAACGCGCCTTCCTGCTGGCCAAGGTCGACGCCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAc                                                                                   <  1:388058/145‑1 (MQ=255)
                                                                   cGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAc                                                                                   <  1:279383/145‑1 (MQ=255)
                                                                   cGAGACATTGCTAGCGCTCGACCACTTACTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTg                                                                                >  1:10574/1‑148 (MQ=255)
                                                                     agaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGaa                                                                              >  7:29517/1‑148 (MQ=255)
                                                                       aCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                             >  5:275900/1‑147 (MQ=255)
                                                                        cATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAg                                                                          <  2:123404/149‑1 (MQ=255)
                                                                         aTTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                         >  2:357717/1‑149 (MQ=255)
                                                                          ttGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGc                                                                        <  8:290648/149‑1 (MQ=255)
                                                                            gCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAg                                                                      >  5:359794/1‑149 (MQ=255)
                                                                              tAGCGCTCGACCACTTGCTGGGCGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGc                                                                     >  7:177978/1‑148 (MQ=255)
                                                                                   cTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTc                                                                                            >  6:147224/1‑120 (MQ=255)
                                                                                   cTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTc                                                                                            <  5:147224/120‑1 (MQ=255)
                                                                                      gACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcgc                                                              <  3:26598/147‑1 (MQ=255)
                                                                                      gACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctg                                                            >  8:13314/1‑149 (MQ=255)
                                                                                       aCCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                    <  1:123345/76‑4 (MQ=37)
                                                                                       aCCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGCCagca                                                                                                                                    >  2:123345/1‑73 (MQ=37)
                                                                                            ttGCTGGGCGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                      <  8:177978/149‑1 (MQ=255)
                                                                                            ttGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                      >  7:239788/1‑149 (MQ=255)
                                                                                            ttGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACtt                                                      >  6:150271/1‑149 (MQ=255)
                                                                                              gCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGACCTGGGCATTTCGCTGGCCACGGTCAAGCGCGACCTGAACAAGGCAGCCCTGCGCGGCTACTTcg                                                    >  6:227958/1‑149 (MQ=255)
                                                                                              gCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcg                                                    >  1:377819/1‑149 (MQ=255)
                                                                                              gCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTc                                                     <  2:172414/148‑1 (MQ=255)
                                                                                               cTGGACGGCCTTCCACCGTGGGTCAACCGCGCCTTCCTGCTGGCCCAGGTCGCCGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcg                                                    <  5:227958/148‑1 (MQ=255)
                                                                                                tGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcgcg                                                  >  2:13933/1‑149 (MQ=255)
                                                                                                  gACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcgcg                                                  <  2:10574/147‑1 (MQ=255)
                                                                                                  gACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCt                                                <  6:120769/149‑1 (MQ=255)
                                                                                                    cGGCCTTCCACCGTTGGTCAAACGCGGCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCTCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                              <  8:230250/149‑1 (MQ=255)
                                                                                                    cGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                              <  4:133071/149‑1 (MQ=255)
                                                                                                    cGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                              <  2:377819/149‑1 (MQ=255)
                                                                                                    cGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                              <  4:164605/149‑1 (MQ=255)
                                                                                                     ggCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                              >  5:6149/1‑148 (MQ=255)
                                                                                                       ccTTCCACCGTTGGTCTAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCt                                                                                                      <  4:320963/90‑1 (MQ=255)
                                                                                                        cTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAAc                                          >  7:192420/1‑149 (MQ=255)
                                                                                                             aCCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGAGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTc                                     >  1:309602/1‑149 (MQ=255)
                                                                                                                   ggTCAAACGCGCCATCCTGCTGGCCCAGGTCGACGGCATGGGCCAAGGCGAAATCGCTCGCGAACTGTGCATTTCGCTGGCCACGGTCAACCGGTACCAGAACAAGGCAGCCATGCGCTGCTACTTAGCGCTGAGAACGCCTCGTTTTc                               >  4:7251/1‑149 (MQ=255)
                                                                                                                    gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTc                               <  3:305759/148‑1 (MQ=255)
                                                                                                                            cgcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggt                      <  5:123042/149‑1 (MQ=255)
                                                                                                                                 tCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg                  >  5:274669/1‑148 (MQ=255)
                                                                                                                                  ccTGCTGGTCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                           >  8:308630/1‑106 (MQ=255)
                                                                                                                                  ccTGCTGGTCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                           <  7:308630/106‑1 (MQ=255)
                                                                                                                                    ccttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa               >  1:25521/4‑148 (MQ=255)
                                                                                                                                     cttgGCCCAGGCCGTTGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCCCGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa               >  2:50656/3‑147 (MQ=255)
                                                                                                                                     cttgGCCCAGGCCGTCGACCTGGGCCAGGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgct                                                             >  7:267978/3‑101 (MQ=38)
                                                                                                                                     cttgGCCCAGGCCGTCGACCTGGGCCAGGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgct                                                             <  8:267978/99‑1 (MQ=38)
                                                                                                                                     cttgGCCCAGGCCGTCGACCTGGACCAGGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgct                                                             <  6:106719/99‑1 (MQ=255)
                                                                                                                                     cttgGCCCAGGCCGTCGACCTGGACCAGGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgct                                                             >  5:106719/3‑101 (MQ=255)
                                                                                                                                         gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGcc         <  3:343647/149‑1 (MQ=255)
                                                                                                                                         gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGc          <  6:6149/148‑1 (MQ=255)
                                                                                                                                         gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGc          <  6:359794/148‑1 (MQ=255)
                                                                                                                                           ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt       <  8:192420/149‑1 (MQ=255)
                                                                                                                                           ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt       >  5:61141/1‑149 (MQ=255)
                                                                                                                                           ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt       >  5:282810/1‑149 (MQ=255)
                                                                                                                                               ggTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCa    <  6:243730/149‑1 (MQ=255)
                                                                                                                                                 tCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCAtt  <  6:275900/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACCTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATT  >  NC_002947/5233399‑5233691

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: