Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,223,668 Δ1 bp 9.1% coding (309/681 nt) PP_2824 ← TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,223,6680G. 9.1% 201.9 / 14.9 66coding (309/681 nt)PP_2824TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base . (3/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TAGCGCTGCTTGATGACACTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCATGGTGGCACCGGCATA  >  NC_002947/3223525‑3223817
                                                                                                                                                 |                                                                                                                                                     
taGCGCTGCTTGATGACACTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCg                                                                                                                                                 >  1:1275254/1‑150 (MQ=255)
 aGCGCTGCTTGATGACACTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTCGATCAGGTCGa                                                                                                                                                >  1:454500/1‑150 (MQ=255)
  gcgcTGCTTGATGACACTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTCGATCAGGTCGac                                                                                                                                               <  2:454500/150‑2 (MQ=255)
  gcgcTGCTTGATGACACTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGAt                                                                                                                                               >  2:760126/1‑150 (MQ=255)
     ctgctTGATGACACTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTc                                                                                                                                            >  2:56227/1‑150 (MQ=255)
                 acTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTc                                                                                                                                >  2:1254233/1‑150 (MQ=255)
                   tGGCGAACGCGGGCTGGTGCCGGGCATAGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCCGATCGAGGTCGATGTCCAGCAGTGCGTCGa                                                                                                                              <  2:1032434/150‑1 (MQ=255)
                        aaCGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTg                                                                                                                         >  1:1195240/1‑150 (MQ=255)
                         aCGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGt                                                                                                                        <  1:589834/150‑1 (MQ=255)
                          cgcgGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTc                                                                                                                       <  1:56227/150‑1 (MQ=255)
                             gggCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGGCTTTGTCGaa                                                                                                                    >  1:1154533/1‑150 (MQ=255)
                             gggCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGaa                                                                                                                    <  2:69203/150‑1 (MQ=255)
                                   gTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGt                                                                                                              <  1:318290/150‑1 (MQ=255)
                                   gTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTCGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGt                                                                                                              >  2:349164/1‑150 (MQ=255)
                                       cGGGCCTGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAAGCTGATCCATTGATGCCACTGCTCAGTAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTCTGTCGAACAGGGTATgg                                                                                                          >  1:99414/1‑149 (MQ=255)
                                       cGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAg                                                                                                                                                     <  2:104457/107‑1 (MQ=255)
                                       cGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAg                                                                                                                                                     >  1:104457/1‑107 (MQ=255)
                                       cGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTgc                                                                                                          <  2:1195240/150‑1 (MQ=255)
                                            aTGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTCGATCAGGTCGac                                                                                                                                               >  2:341038/1‑107 (MQ=255)
                                            aTGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTCGATCAGGTCGac                                                                                                                                               <  1:341038/108‑2 (MQ=255)
                                                    gCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTCGATCAGGTCGaca                                                                                                                                              <  1:955558/101‑3 (MQ=255)
                                                    gCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTCGATCAGGTCGaca                                                                                                                                              >  2:955558/1‑99 (MQ=255)
                                                    gCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACa                                                                                                                  <  2:1276017/129‑1 (MQ=255)
                                                    gCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACa                                                                                                                  >  1:1276017/1‑129 (MQ=255)
                                                        ttGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCAt                                                                                         >  1:878807/1‑150 (MQ=255)
                                                            aGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATttct                                                                                     >  2:1279773/1‑150 (MQ=255)
                                                            aGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATttct                                                                                     >  1:894856/1‑150 (MQ=255)
                                                              ccATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATttcttc                                                                                   <  2:388395/150‑1 (MQ=255)
                                                                  gCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCttg                                                                               >  2:702778/1‑150 (MQ=255)
                                                                    aaGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTc                                                                                                                                >  2:71210/1‑99 (MQ=255)
                                                                    aaGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTc                                                                                                                                <  1:71210/99‑1 (MQ=255)
                                                                     aGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGAAACGGCACTGACCATTTCTTCCttgttg                                                                            <  2:516708/150‑1 (MQ=255)
                                                                      gCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTgg                                                                           <  1:650415/150‑1 (MQ=255)
                                                                          gTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGaa                                                                       <  2:878807/150‑1 (MQ=255)
                                                                          gTTGCATGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGCCTCGAAACGGCAAGTCCGGTTTTCGATCTCCAGCAGTTTGTCGACTTTGTCGAATAGGGTCTTCGCAACTGCACTGACCATTTCTTCCTTGTTGGCGaa                                                                       <  1:1254233/150‑1 (MQ=255)
                                                                           ttGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAag                                                                      <  1:265294/150‑1 (MQ=255)
                                                                              ccTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAagtag                                                                   >  2:388526/1‑150 (MQ=255)
                                                                               cTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAagtagt                                                                  <  1:858736/150‑1 (MQ=255)
                                                                                 ggCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGtagtag                                                                >  2:23973/1‑150 (MQ=255)
                                                                                  gCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAgg                                                               <  1:388526/150‑1 (MQ=255)
                                                                                  gCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAgg                                                               >  2:1119495/1‑150 (MQ=255)
                                                                                    tCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTc                                                             >  1:758651/1‑150 (MQ=255)
                                                                                       tCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCAcc                                                          >  2:1032248/1‑150 (MQ=255)
                                                                                               gATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCtgt                                                  <  1:1032248/150‑1 (MQ=255)
                                                                                                 tCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTgtgg                                                <  1:1119495/150‑1 (MQ=255)
                                                                                                 tCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTgtgg                                                >  1:1221220/1‑150 (MQ=255)
                                                                                                               tGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGtgc                                  >  2:1025237/1‑150 (MQ=255)
                                                                                                                 cTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCtgtg                                                 >  1:502673/1‑133 (MQ=255)
                                                                                                                 cTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCtgtg                                                 <  2:502673/133‑1 (MQ=255)
                                                                                                                  tCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGtgctgc                               <  1:760126/150‑1 (MQ=255)
                                                                                                                     aTAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGcc                            >  1:361474/1‑150 (MQ=255)
                                                                                                                      tAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCgcgc                                                      <  2:274132/123‑1 (MQ=255)
                                                                                                                      tAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCgcgc                                                      >  1:274085/1‑123 (MQ=255)
                                                                                                                      tAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCgcgc                                                      <  2:274085/123‑1 (MQ=255)
                                                                                                                      tAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCgcgc                                                      >  1:274132/1‑123 (MQ=255)
                                                                                                                            aCTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTgc                     <  1:1279773/150‑1 (MQ=255)
                                                                                                                            aCTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTgc                     <  1:599986/150‑1 (MQ=255)
                                                                                                                              ttGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGaa                                                                       <  2:910527/98‑1 (MQ=255)
                                                                                                                              ttGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGaa                                                                       >  1:910527/1‑98 (MQ=255)
                                                                                                                                  tGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCAtgg               <  1:702778/150‑1 (MQ=255)
                                                                                                                                    tCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCAtggtg             <  2:1221220/150‑1 (MQ=255)
                                                                                                                                    tCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCAtggtg             >  2:349435/1‑150 (MQ=255)
                                                                                                                                    tCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCCCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCAtggtg             <  1:349435/150‑1 (MQ=255)
                                                                                                                                       aCCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACGTTGTGGAACAGGGTCTGCGCAAGGGCACTGACGTTTTCTTACTTGTCGCCGAATTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCAACCACGCACTTGGTGGCa          >  1:732255/1‑150 (MQ=255)
                                                                                                                                         cTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCATGGTGGCAcc        <  2:361474/150‑1 (MQ=255)
                                                                                                                                            gATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCTCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCATGGTGGCACCgg      >  2:1126209/1‑150 (MQ=255)
                                                                                                                                                 gAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGAGCAACGGCACTGACTATTTCTTCCTTGTTGGAGCAGTAGCAGGTCACTTCTCTGATTTTGCATCGTGCGTTCTGCGCCACCTTGCGCATGGTGGCACCGGCATa  <  1:23973/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TAGCGCTGCTTGATGACACTGGCGAACGCGGGCTGGTGCCGGGCATGGGTGAGCAGTTGCAGCCATGCAAGCCAGTTGCCTGGCTCGTCAAGGCTGATCCATTGGTGCCACTGCTCAATAAGCGACTTGATGTCGACCTGATCGAGGTCGATGTCCAGCAGTGCGTCGACTTTGTCGAACAGGGTCTGCGCAACGGCACTGACCATTTCTTCCTTGTTGGCGAAGTAGTAGGTCACCGCGCCTGTGGTGCACCCGGCGTGCTGCGCCACCTTGCGCATGGTGGCACCGGCATA  >  NC_002947/3223525‑3223817

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: