Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,199,450 C→G 10.5% A36A (GCG→GCC 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,4500CG10.5% 257.0 / 18.8 86A36A (GCG→GCCPP_2806TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (39/38);  new base G (6/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGCCCCTT  >  NC_002947/3199303‑3199595
                                                                                                                                                   |                                                                                                                                                  
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)
                                      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)
                                               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:2184865/1‑150 (MQ=255)
                                                                         tCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcc                                                                        >  2:1678764/1‑150 (MQ=255)
                                                                           gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCgg                                                                                                                                    >  1:1027956/1‑88 (MQ=255)
                                                                           gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCgg                                                                                                                                    <  2:1027956/88‑1 (MQ=255)
                                                                           gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCAt                                                                                                              <  2:723714/110‑1 (MQ=255)
                                                                           gTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCAt                                                                                                              >  1:723714/1‑110 (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                                             <  2:234101/150‑1 (MQ=255)
                                                                                                    ttGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                             <  1:189904/150‑1 (MQ=255)
                                                                                                    ttGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                             >  1:1249869/1‑150 (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                                             <  2:900866/142‑1 (MQ=255)
                                                                                                            gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                             >  1:900866/1‑142 (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                                                                                        >  2:96261/2‑76 (MQ=25)
                                                                                                                                   ccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGt                                                                                        <  1:96261/75‑1 (MQ=25)
                                                                                                                                     ggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTgccg            >  1:269527/1‑150 (MQ=255)
                                                                                                                                      gcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTgccgc           >  1:284848/1‑150 (MQ=255)
                                                                                                                                         gaatcgctgcGGCAGCAGATTCGCCGGGTGAATTCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGccgccg         >  1:62698/6‑150 (MQ=255)
                                                                                                                                          aTCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCAGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGc       <  1:1743285/150‑1 (MQ=255)
                                                                                                                                               ctgcCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTc                                                          <  1:1817325/94‑1 (MQ=255)
                                                                                                                                               ctgcCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTc                                                          >  2:1817325/1‑94 (MQ=255)
                                                                                                                                               ctgcCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGCCCCtt  <  2:527392/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGCCCCTT  >  NC_002947/3199303‑3199595

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: