Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,414,219 Δ1 bp 5.6% coding (256/561 nt) PP_1236 ← glycine cleavage system transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,414,2190T. 5.6% 277.1 / 10.2 90coding (256/561 nt)PP_1236glycine cleavage system transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (47/35);  new base . (2/3);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

CAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCA  >  NC_002947/1414073‑1414362
                                                                                                                                                   |                                                                                                                                               
cAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCg                                                                                                                                              >  2:2038796/1‑150 (MQ=255)
 aggtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGc                                                                                                                                               >  1:750880/1‑148 (MQ=255)
 aggtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGc                                                                                                                                               <  2:750880/148‑1 (MQ=255)
 aggtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTCGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCgg                                                                                                                                             >  1:1013554/1‑150 (MQ=255)
   gtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGAc                                                                                                                                           >  2:1361977/1‑150 (MQ=255)
     aggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCa                                                                                                                                         <  1:1465974/150‑1 (MQ=255)
       gtgtCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCagg                                                                                                                                       <  2:1406358/150‑3 (MQ=255)
         gtCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCaggaa                                                                                                                                     <  2:1603207/150‑5 (MQ=255)
            gCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCGAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAg                                                                                                                                  >  1:2187855/1‑150 (MQ=255)
             cAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAgg                                                                                                                                 >  2:1180221/1‑150 (MQ=255)
              aGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAggg                                                                                                                                >  1:1799690/1‑150 (MQ=255)
               ggTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGt                                                                                                                               <  1:1291261/150‑1 (MQ=255)
                  cATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAg                                                                                                                            >  1:324946/1‑150 (MQ=255)
                     ggCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCTGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGcc                                                                                                                         >  2:2024080/1‑150 (MQ=255)
                     ggCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCaggaa                                                                                                                                     <  2:1282054/138‑5 (MQ=255)
                     ggCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCaggaa                                                                                                                                     >  1:1282054/1‑134 (MQ=255)
                     ggCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGcc                                                                                                                         >  2:2166770/1‑150 (MQ=255)
                        ttCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTg                                                                                                                      >  2:1457255/1‑150 (MQ=255)
                            aGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGcttc                                                                                                                  <  2:1995218/150‑1 (MQ=255)
                               tCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTg                                                                                                               >  1:1734494/1‑150 (MQ=255)
                                 gACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGcc                                                                                                             >  2:1930158/1‑150 (MQ=255)
                                    gcgGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCaa                                                                                                          <  1:1457255/150‑1 (MQ=255)
                                      ggtggtCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGc                                                                                                        <  2:1517887/150‑1 (MQ=255)
                                       gtggtCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGcc                                                                                                       >  2:2040040/1‑150 (MQ=255)
                                         ggtCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTg                                                                                                     <  2:404114/150‑1 (MQ=255)
                                           tCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGc                                                                                                   >  2:1122714/1‑150 (MQ=255)
                                             aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAg                                                                                                 >  1:2069904/1‑150 (MQ=255)
                                             aaGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGACTGCGGGCGCACTTCCTGGTCCGCGCAGCGGCCCACCCTGAGCGCCAGTGCGTGCTTCTTGCCCAAGCCTGGCAg                                                                                                 >  1:1321360/1‑150 (MQ=255)
                                              aGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACACAGGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAgg                                                                                                >  2:563489/1‑150 (MQ=255)
                                               ggAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTg                                                                                                               <  1:885587/134‑1 (MQ=255)
                                               ggAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTg                                                                                                               >  2:885587/1‑134 (MQ=255)
                                               ggAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAggg                                                                                               >  2:1964193/1‑150 (MQ=255)
                                               ggAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAggg                                                                                               >  1:1483947/1‑150 (MQ=255)
                                                taagaaCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCTGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGt                                                                                              <  1:78618/149‑1 (MQ=255)
                                                gaagaaCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGt                                                                                              >  2:992721/1‑150 (MQ=255)
                                                gaagaaCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGt                                                                                              >  1:269172/1‑150 (MQ=255)
                                                   gaaCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGa                                                                                           <  2:324946/150‑1 (MQ=255)
                                                      cTGGCACAGTTTGTTTATGATGTCCAGGGGGGAGGCAGCGCTAACATAGGCCACCCACCGAAAGCCCTGCGCGCGCACTTCCTCGTCCACGCCGCGCCCCTCGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTc                                                                                        <  2:723273/150‑1 (MQ=255)
                                                        ggCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCga                                                                                      <  2:306421/150‑1 (MQ=255)
                                                         gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCg                                                                                                                                    >  1:1932958/1‑103 (MQ=255)
                                                         gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCg                                                                                                                                    <  2:1932958/103‑1 (MQ=255)
                                                         gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGa                                                                                                                                   <  1:1224933/104‑1 (MQ=255)
                                                         gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGa                                                                                                                                   >  2:1224933/1‑104 (MQ=255)
                                                         gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCaa                                                                                                          <  2:1828940/129‑1 (MQ=255)
                                                         gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCaa                                                                                                          >  1:1828940/1‑129 (MQ=255)
                                                         gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCgag                                                                                     >  2:2273667/1‑150 (MQ=255)
                                                          cacaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCgaga                                                                                    <  2:1013554/150‑1 (MQ=255)
                                                           acaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCa                                                                                                                                         >  2:1830231/1‑96 (MQ=255)
                                                           acaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCa                                                                                                                                         <  1:1830231/96‑1 (MQ=255)
                                                           acaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCt                                                                                                                                                 <  2:635688/88‑1 (MQ=255)
                                                           acaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGAc                                                                                   >  1:1673954/1‑150 (MQ=255)
                                                           acaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACGTCCTGGTCCGCGCt                                                                                                                                                 >  1:635688/1‑88 (MQ=255)
                                                                  gTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTc                                                                                                                              >  2:2156691/1‑100 (MQ=255)
                                                                  gTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAAGGCCTGCGTGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTc                                                                                                                              <  1:2156691/100‑1 (MQ=255)
                                                                   ttGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           >  2:123651/1‑150 (MQ=255)
                                                                   ttGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           >  1:1143507/1‑150 (MQ=255)
                                                                           gTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCa                                                                   <  1:1361977/150‑1 (MQ=255)
                                                                              cGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGcttc                                                                                                                  >  1:1584942/1‑100 (MQ=255)
                                                                              cGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGcttc                                                                                                                  <  2:1584942/100‑1 (MQ=255)
                                                                              cGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCt                                                                <  1:1180221/150‑1 (MQ=255)
                                                                                   ggTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTgccgc                                                           >  1:1509031/1‑150 (MQ=255)
                                                                                   ggTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTACGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTgccgc                                                           >  2:1147502/1‑150 (MQ=255)
                                                                                        gCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCAcc                                                      <  1:2057410/150‑1 (MQ=255)
                                                                                             gcTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcag                                                 >  2:160466/1‑150 (MQ=255)
                                                                                              cTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcagc                                                <  2:1799690/150‑1 (MQ=255)
                                                                                               tCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTgcagca                                               >  2:1147672/1‑150 (MQ=255)
                                                                                                     aGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTATTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAggg                                         <  1:164765/150‑1 (MQ=255)
                                                                                                       gCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCa                                                                                                                                         <  2:96887/52‑1 (MQ=37)
                                                                                                       gCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCa                                                                                                                                         >  1:96887/1‑52 (MQ=37)
                                                                                                        ccACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGgcgc                                      >  1:200346/1‑150 (MQ=255)
                                                                                                         cacaTACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAggg                                         <  2:2104345/146‑1 (MQ=255)
                                                                                                         cacaTACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAggg                                         >  1:2104345/1‑146 (MQ=255)
                                                                                                          acaTACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTg                                    >  1:1117486/1‑150 (MQ=255)
                                                                                                              caggAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGtct                                >  1:385161/3‑150 (MQ=255)
                                                                                                               cGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGtctc                               <  1:2040040/150‑1 (MQ=255)
                                                                                                                 gAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGc                             >  2:1719248/1‑150 (MQ=255)
                                                                                                                  aaGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           <  1:1381725/103‑1 (MQ=255)
                                                                                                                  aaGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           >  2:1381725/1‑103 (MQ=255)
                                                                                                                       ccTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggc                       <  2:165409/150‑1 (MQ=255)
                                                                                                                        cTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcg                      >  2:1012792/1‑150 (MQ=255)
                                                                                                                         tGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGTCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcgg                     <  2:2069904/150‑1 (MQ=255)
                                                                                                                         tGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGCCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTggcgg                     <  2:1321360/150‑1 (MQ=255)
                                                                                                                            gggCGCACTTCCTGGTCCGCGCTGCGGACCCCGTCGAGGGTCAGGGCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCATCTGCCGCCCCCCTGCAGCACCAGGGCGCTGGCCTCGCCGTGGCGGCTc                  >  1:456651/1‑150 (MQ=255)
                                                                                                                             ggCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCa                 <  1:1122714/150‑1 (MQ=255)
                                                                                                                                gcACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGAGCGCTGGTCTCGCCGTGGCGGCTCAgtc              >  2:15983/1‑148 (MQ=255)
                                                                                                                                  aCTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGgcgc            >  2:1368915/1‑150 (MQ=255)
                                                                                                                                     tCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTTTCTCCTTTGCGGCTCAGGCGCGAg         <  1:2166770/150‑1 (MQ=255)
                                                                                                                                     tCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAg         <  2:1143507/150‑1 (MQ=255)
                                                                                                                                          gTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGac    <  2:2187855/150‑1 (MQ=255)
                                                                                                                                            ccGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGacca  <  2:1117486/150‑1 (MQ=255)
                                                                                                                                              gcgcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           <  2:1503351/75‑1 (MQ=255)
                                                                                                                                              gcgcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           >  1:1503351/1‑75 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCA  >  NC_002947/1414073‑1414362

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