Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,199,444 Δ1 bp 9.8% coding (114/582 nt) PP_2806 ← TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,199,4440T. 9.8% 246.0 / 20.2 82coding (114/582 nt)PP_2806TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base . (5/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

ACATTTCCCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGC  >  NC_002947/3199296‑3199590
                                                                                                                                                    |                                                                                                                                                   
aCATTTCCCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATctg                                                                                                                                                    >  2:1268082/1‑150 (MQ=255)
       ccccACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGc                                                                                                                                             >  2:851463/1‑150 (MQ=255)
       ccccACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCGCTGCGGCa                                                                                                                                            >  1:1061863/1‑150 (MQ=255)
       ccccACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGAACTTGTGGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGc                                                                                                                                             >  1:523256/1‑150 (MQ=255)
         ccACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAg                                                                                                                                           >  2:2055470/1‑150 (MQ=255)
            ccGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGa                                                                                                                                       <  1:1533010/150‑1 (MQ=255)
             cGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGAt                                                                                                                                      <  1:2052613/150‑1 (MQ=255)
               aTTCCGGCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTc                                                                                                                                    <  1:704518/150‑1 (MQ=255)
               aTTCCGCCGTCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTc                                                                                                                                    <  1:851463/150‑1 (MQ=255)
               aTTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCGCTGCGGCAGCAGATTc                                                                                                                                    <  2:1061863/150‑1 (MQ=255)
                   cgccgccgTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCg                                                                                                                                <  1:2071839/150‑1 (MQ=255)
                     ccgccgTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCggg                                                                                                                              >  2:1543561/1‑150 (MQ=255)
                     ccgccgTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCggg                                                                                                                              >  2:1140249/1‑150 (MQ=255)
                         cgTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTtgt                                                                                                                          <  1:1342913/150‑1 (MQ=255)
                          gTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGCCTAGCTGAGATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTtgtg                                                                                                                         <  2:919043/150‑1 (MQ=255)
                                 cTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTgg                                                                                                                  <  1:1653676/150‑1 (MQ=255)
                                   gTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATctgctg                                                                                                                                                 <  2:823712/118‑1 (MQ=255)
                                   gTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATctgctg                                                                                                                                                 >  1:823712/1‑118 (MQ=255)
                                             gtTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCANGtt                                                                                                      <  2:1514474/150‑1 (MQ=255)
                                                gCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCcg                                                                                                   >  1:452158/1‑150 (MQ=255)
                                                 cAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCcgc                                                                                                  <  1:2037091/150‑1 (MQ=255)
                                                  aaCGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCTCTGTGTTGTGGCGTTGGCCACCCCATGTTCcgcg                                                                                                 >  2:189904/1‑150 (MQ=255)
                                                  aaCGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCcgcg                                                                                                 >  1:1800741/1‑150 (MQ=255)
                                                   aCGATTGTTGCGGCATGTTGATCAAGGTCTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTACGCGa                                                                                                >  1:322363/1‑150 (MQ=255)
                                                   aCGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATctgct                                                                                                                                                  <  2:987451/101‑1 (MQ=255)
                                                   aCGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATctgct                                                                                                                                                  >  1:987451/1‑101 (MQ=255)
                                                      aTTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCa                                                                                             >  1:1435236/1‑150 (MQ=255)
                                                          ttgCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTTCTTGGTGTGGAAGACATAGTGCGACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCtt                                                                                         >  1:1747898/1‑150 (MQ=255)
                                                                cATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGt                                                                                   >  2:1743285/1‑150 (MQ=255)
                                                                 aTGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAAGCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTg                                                                                  <  1:997769/150‑1 (MQ=255)
                                                                    ttGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCTCGATCACCTTCACGGTGGcc                                                                               >  2:1764372/1‑150 (MQ=255)
                                                                      gATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTc                                                                             <  2:2031615/150‑1 (MQ=255)
                                                                        tCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGa                                                                           <  1:2055470/150‑1 (MQ=255)
                                                                            ggACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATaaa                                                                       <  1:2317288/150‑1 (MQ=255)
                                                                             gACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATaaa                                                                       >  2:1591280/1‑149 (MQ=255)
                                                                             gACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATaaa                                                                       <  1:1591280/149‑1 (MQ=255)
                                                                                tCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcc                                                                   >  2:1678764/1‑150 (MQ=255)
                                                                                tCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcc                                                                   >  2:2184865/1‑150 (MQ=255)
                                                                                  gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCgg                                                                                                                               <  2:1027956/88‑1 (MQ=255)
                                                                                  gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCgg                                                                                                                               >  1:1027956/1‑88 (MQ=255)
                                                                                  gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCAt                                                                                                         >  1:723714/1‑110 (MQ=255)
                                                                                  gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCAt                                                                                                         <  2:723714/110‑1 (MQ=255)
                                                                                  gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcctg                                                                 >  1:234101/1‑150 (MQ=255)
                                                                                  gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcctg                                                                 >  1:370534/1‑150 (MQ=255)
                                                                                      aCCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCa                                                                                                                                            <  2:2129592/72‑1 (MQ=255)
                                                                                      aCCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCa                                                                                                                                            >  1:2129592/1‑72 (MQ=255)
                                                                                            tAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGAAACGATTCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGctct                                                       <  1:982634/150‑1 (MQ=255)
                                                                                                    cTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCgg                                                                                                                               <  2:2212312/70‑1 (MQ=255)
                                                                                                    cTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCgg                                                                                                                               >  1:2212312/1‑70 (MQ=255)
                                                                                                    cTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGc                                                  >  1:1147604/1‑147 (MQ=255)
                                                                                                    cTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGc                                                  <  2:1147604/147‑1 (MQ=255)
                                                                                                       gTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAAt                                            >  1:628415/1‑150 (MQ=255)
                                                                                                       gTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAAt                                            <  2:1435236/150‑1 (MQ=255)
                                                                                                         ccTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATcc                                          >  2:2106228/1‑150 (MQ=255)
                                                                                                           ttGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                        >  1:1249869/1‑150 (MQ=255)
                                                                                                           ttGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                        <  1:189904/150‑1 (MQ=255)
                                                                                                           ttGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                        <  2:234101/150‑1 (MQ=255)
                                                                                                            tGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAg                                       <  1:1268082/150‑1 (MQ=255)
                                                                                                                 tgGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCAtt                                  >  2:2372582/1‑150 (MQ=255)
                                                                                                                   gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCa                                                 >  1:2214086/1‑133 (MQ=255)
                                                                                                                   gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCa                                                 <  2:2214086/133‑1 (MQ=255)
                                                                                                                   gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                        >  1:900866/1‑142 (MQ=255)
                                                                                                                   gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                        <  2:900866/142‑1 (MQ=255)
                                                                                                                   gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGa                                >  2:2120241/1‑150 (MQ=255)
                                                                                                                   gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAACGATTCGCGGGGTTGTGGCGTTGGGTACTCCATGTTACGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAAGCCGAGCCATTGa                                >  1:1994032/1‑150 (MQ=255)
                                                                                                                   gAAGACATAGTGCAACGATGCCAGCGGCGAATATGGAAACGCCGCGGTTCGGCGGTTGGTGGCGTTGGCGGCCCCATGTTCTACTCTTACCTTCACTGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGa                                >  1:1749932/1‑150 (MQ=255)
                                                                                                                    aaGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAg                               >  2:1925066/1‑150 (MQ=255)
                                                                                                                     aGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGc                              <  2:1481077/150‑1 (MQ=255)
                                                                                                                     aGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGc                              <  1:2372582/150‑1 (MQ=255)
                                                                                                                      gACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCg                             <  1:1140249/150‑1 (MQ=255)
                                                                                                                         aTAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTtc                          >  2:2052622/1‑150 (MQ=255)
                                                                                                                            gTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTcc                       <  1:2052622/150‑1 (MQ=255)
                                                                                                                               cAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCtt                    <  2:1747898/150‑1 (MQ=255)
                                                                                                                               cAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCtt                    <  1:2106228/150‑1 (MQ=255)
                                                                                                                                   gATGCCAGCGGCGAATCTTCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCg                >  1:1079174/1‑150 (MQ=255)
                                                                                                                                   gATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCg                >  1:1629500/1‑150 (MQ=255)
                                                                                                                                     tGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGCTGTGGCGTTGGCCACCCCATGTTCCGCGACCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCCCCAACAACCCGACCCACAGAGCGTTCTCCCTTTTCGGt              >  2:514672/1‑150 (MQ=255)
                                                                                                                                       ccAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTGCGGGAt            <  2:628415/150‑1 (MQ=255)
                                                                                                                                       aacccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTgg                                                                                 <  1:2135437/77‑1 (MQ=21)
                                                                                                                                       aacccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTgg                                                                                 >  2:2135437/5‑81 (MQ=21)
                                                                                                                                        cAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTAtt           >  1:527392/1‑150 (MQ=255)
                                                                                                                                        cAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTAtt           <  1:1543561/150‑1 (MQ=255)
                                                                                                                                          ccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGt                                                                                   <  1:96261/75‑1 (MQ=25)
                                                                                                                                          ccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGt                                                                                   >  2:96261/2‑76 (MQ=25)
                                                                                                                                            ggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTgccg       >  1:269527/1‑150 (MQ=255)
                                                                                                                                             gcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTgccgc      >  1:284848/1‑150 (MQ=255)
                                                                                                                                                 aTCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCAGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGc  <  1:1743285/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACATTTCCCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGC  >  NC_002947/3199296‑3199590

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: