Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,199,444 Δ1 bp 8.9% 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. 8.9% 233.5 / 18.0 79coding (114/582 nt)PP_2806TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (41/30);  new base . (4/3);  total (46/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 = 1.81e-01

TTTCCCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGCC  >  NC_002947/3199299‑3199591
                                                                                                                                                 |                                                                                                                                                    
tttCCCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATctgctg                                                                                                                                                  >  2:2396587/1‑150 (MQ=255)
       cACCGATTCCGCCGCAGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGc                                                                                                                                           <  2:2207453/150‑1 (MQ=255)
        aCCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCg                                                                                                                                         >  2:1830963/1‑150 (MQ=255)
         ccGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGa                                                                                                                                        <  1:1435408/150‑1 (MQ=255)
           gATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTCGGGGTGGAAGACATCGTGCCGCGTCGTCAACGGCGAATGTGCTGCCGCAGCAAtt                                                                                                                                      >  2:443134/1‑150 (MQ=255)
              tCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGc                                                                                                                                   <  2:53830/150‑1 (MQ=255)
                 gccgccgTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGTGGCGAATCTGCTGCCGCAGCGATTCGCCgg                                                                                                                                <  1:676754/150‑1 (MQ=255)
                       gTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTtgtg                                                                                                                          <  2:948043/150‑1 (MQ=255)
                            gcCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGtt                                                                                                                     >  2:185269/1‑150 (MQ=255)
                            gcCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCGCTGCGGCAGCAGATTc                                                                                                                                     <  2:2048052/134‑1 (MQ=255)
                            gcCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCGCTGCGGCAGCAGATTc                                                                                                                                     >  1:2048052/1‑134 (MQ=255)
                                   ggCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGtt                                                                                                                             <  1:161089/135‑1 (MQ=255)
                                   ggCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGtt                                                                                                                             >  2:161089/1‑135 (MQ=255)
                                        gtgtTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATg                                                                                                         >  1:1183887/1‑150 (MQ=255)
                                             gCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAAACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCcg                                                                                                    >  2:2079424/1‑150 (MQ=255)
                                               aaCGATTGTTGCGGCATGTTGATCAAGGACTCTTAAACCGCATAGAACAGCTCGCACTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCcgcg                                                                                                  >  2:811512/1‑150 (MQ=255)
                                                aCGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGa                                                                                                 >  2:2113783/1‑150 (MQ=255)
                                                   aTTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCa                                                                                              <  1:2113783/150‑1 (MQ=255)
                                                      gttgCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCt                                                                                           <  2:1183887/150‑1 (MQ=255)
                                                       ttgCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCtt                                                                                          >  2:211965/1‑150 (MQ=255)
                                                         gCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCa                                                                                        >  2:842158/1‑150 (MQ=255)
                                                                gTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGc                                                                                 <  1:1830963/150‑1 (MQ=255)
                                                                     tCAAGGACTCGTAGACCCCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCAccc                                                                                                             <  2:1203774/117‑1 (MQ=255)
                                                                     tCAAGGACTCGTAGACCCCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCAccc                                                                                                             >  1:1203774/1‑117 (MQ=255)
                                                                     caacggACTCGCAGACCTCATAGAACATCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGa                                                                            >  2:890910/5‑150 (MQ=255)
                                                                       aaGGACTCGTAGACCGCATAGAACAGCTCGTCCTCGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCCGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGCTa                                                                          >  1:861260/1‑150 (MQ=255)
                                                                        aGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATaa                                                                         >  2:1864010/1‑150 (MQ=255)
                                                                         ggACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATaaa                                                                        >  1:904479/1‑150 (MQ=255)
                                                                             tCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGTGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcc                                                                    <  1:1449890/150‑1 (MQ=255)
                                                                             tCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGAGGCCCGCGGCCAATCTGGTGTCGCAGCGATGCGCCGGGTTGTGGCGTTGGCCGCCCCCGGTTCCGCGACCGCCTTCACGGTGGCCTCGATAAAGTcc                                                                    >  2:2077066/1‑150 (MQ=255)
                                                                                  gACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcctgcct                                                               <  1:842158/150‑1 (MQ=255)
                                                                                   aCCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTtgtg                                                                                                                          <  1:364235/90‑1 (MQ=255)
                                                                                   aCCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTtgtg                                                                                                                          >  2:364235/1‑90 (MQ=255)
                                                                                     cGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcctgcctccc                                                            >  2:358227/1‑147 (MQ=255)
                                                                                     cGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcct                                                                   <  1:1050075/143‑1 (MQ=255)
                                                                                     cGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTcct                                                                   >  2:1050075/1‑143 (MQ=255)
                                                                                        aTAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCTCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGctc                                                         >  1:1076733/1‑150 (MQ=255)
                                                                                        aTAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGctc                                                         <  1:211965/150‑1 (MQ=255)
                                                                                          aGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCtt                                                       <  1:1692917/150‑1 (MQ=255)
                                                                                            aatggCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCa                                                     >  2:1569276/5‑150 (MQ=255)
                                                                                            aaCAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCa                                                     >  1:476865/1‑150 (MQ=255)
                                                                                             aCAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCGGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAg                                                    <  2:889970/150‑1 (MQ=255)
                                                                                             aCAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAg                                                    <  1:890910/150‑1 (MQ=255)
                                                                                             aCAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGACGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAg                                                    >  1:326850/1‑150 (MQ=255)
                                                                                               gtttCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCtt                                                       <  1:2238514/142‑1 (MQ=255)
                                                                                               gtttCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCGCtt                                                       >  2:2238514/4‑145 (MQ=255)
                                                                                                gCTCTTCCTTAGTGGGGAATACATAGTGCACCGATGCCAGCTGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCAc                                                 <  2:1076733/150‑1 (MQ=255)
                                                                                                  tCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCa                                               >  2:376912/1‑150 (MQ=255)
                                                                                                    gTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTc                                                                                                                                     <  1:1793757/62‑1 (MQ=255)
                                                                                                    gTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTc                                                                                                                                     >  2:1793757/1‑62 (MQ=255)
                                                                                                    gTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAAt                                             >  2:2285846/1‑150 (MQ=255)
                                                                                                    gTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAAt                                             >  2:483706/1‑150 (MQ=255)
                                                                                                     tCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGa                                                                                                 >  1:1110998/1‑97 (MQ=255)
                                                                                                     tCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGa                                                                                                 <  2:1110998/97‑1 (MQ=255)
                                                                                                        ttGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGTGAATGGTATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                         <  2:64611/150‑1 (MQ=255)
                                                                                                        ttGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGa                                         <  1:1864010/150‑1 (MQ=255)
                                                                                                             gttgAAGACAAAGCGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCAt                                    <  2:476865/147‑1 (MQ=255)
                                                                                                              tgGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTgg                                                                                                                   >  2:1308224/1‑70 (MQ=255)
                                                                                                              tgGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTgg                                                                                                                   <  1:1308224/70‑1 (MQ=255)
                                                                                                              tgGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCAtt                                   <  1:1569276/150‑1 (MQ=255)
                                                                                                              tgGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCAtt                                   <  1:1917275/150‑1 (MQ=255)
                                                                                                                gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGa                                 >  1:2271360/1‑150 (MQ=255)
                                                                                                                gAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGa                                 >  2:2251304/1‑150 (MQ=255)
                                                                                                                  aGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCTCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGc                               >  2:1392032/1‑150 (MQ=255)
                                                                                                                      aTAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTtc                           >  2:2137372/1‑150 (MQ=255)
                                                                                                                         gTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGGCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTCCATGCGCACTTCAGCACCAATCCGAGCCATTGAGCGTTCTcc                        >  2:446289/1‑150 (MQ=255)
                                                                                                                         gTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTcc                        >  1:834357/1‑150 (MQ=255)
                                                                                                                             aaCGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCttt                    >  1:317735/1‑150 (MQ=255)
                                                                                                                             aaCGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCttt                    >  1:884502/1‑150 (MQ=255)
                                                                                                                              aCGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCtttt                   >  1:659177/1‑150 (MQ=255)
                                                                                                                                 aTGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCg                 <  1:1898486/149‑1 (MQ=255)
                                                                                                                                 aTGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCg                 >  2:1898486/1‑149 (MQ=255)
                                                                                                                                   caacccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGctct                                                        <  2:125975/103‑1 (MQ=37)
                                                                                                                                   caacccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGctct                                                        >  1:125975/6‑108 (MQ=37)
                                                                                                                                     cAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCAGCTGGGCCTCGATAAAGTCCTGCCTGAGATCTTGAGCAGCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTAtt            >  1:2340762/1‑150 (MQ=255)
                                                                                                                                       ccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCtt                                                       >  1:688315/2‑105 (MQ=38)
                                                                                                                                       ccggcgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCtt                                                       <  2:688315/104‑1 (MQ=38)
                                                                                                                                        cggagAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTgcc         <  1:446289/150‑1 (MQ=255)
                                                                                                                                           cgAATCGCTGCGGCAGCAGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTgccgcc      >  2:376437/1‑150 (MQ=255)
                                                                                                                                               tCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGcc  >  1:821364/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTTCCCCCACCGATTCCGCCGCCGTGGCGCCTGTCGGCACGTGTTGCAACGATTGTTGCGGCATGTTGATCAAGGACTCGTAGACCGCATAGAACAGCTCGTCCTTGGTGTGGAAGACATAGTGCAACGATGCCAGCGGCGAATCTGCTGCCGCAGCGATTCGCCGGGTTGTGGCGTTGGCCACCCCATGTTCCGCGATCACCTTCACGGTGGCCTCGATAAAGTCCTGCCTGCGCTCTTCAGCACCAATCCGAGCCATTGAGCGTTCTCCCTTTTCGGTATTGCCGCCGGCC  >  NC_002947/3199299‑3199591

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