Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,414,216:1 +A 5.6% coding (259/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,2161.A 5.6% 370.3 / 15.2 90coding (259/561 nt)PP_1236glycine cleavage system transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (48/37);  new base A (2/3);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GGCGCCAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCA  >  NC_002947/1414068‑1414362
                                                                                                                                                     |                                                                                                                                                  
ggcgCCAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCcgcc                                                                                                                                                    >  1:1561211/1‑149 (MQ=255)
     cAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCg                                                                                                                                              >  2:2038796/1‑150 (MQ=255)
      aggtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGc                                                                                                                                               <  2:750880/148‑1 (MQ=255)
      aggtaggtGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGc                                                                                                                                               >  1:750880/1‑148 (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                                                                                                                                     >  1:1282054/1‑134 (MQ=255)
                          ggCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCaggaa                                                                                                                                     <  2:1282054/138‑5 (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                                                                                              >  1:269172/1‑150 (MQ=255)
                                                     gaagaaCTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGt                                                                                              >  2:992721/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                                                                                                                                    <  2:1932958/103‑1 (MQ=255)
                                                              gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCg                                                                                                                                    >  1:1932958/1‑103 (MQ=255)
                                                              gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGa                                                                                                                                   >  2:1224933/1‑104 (MQ=255)
                                                              gCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGa                                                                                                                                   <  1:1224933/104‑1 (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                                                                                                                                         <  1:1830231/96‑1 (MQ=255)
                                                                acaGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCAGCGCGGACCa                                                                                                                                         >  2:1830231/1‑96 (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                                                                           >  1:1143507/1‑150 (MQ=255)
                                                                        ttGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           >  2:123651/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                                                                           >  1:1503351/1‑75 (MQ=255)
                                                                                                                                                   gcgcTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAgg                                                                           <  2:1503351/75‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGCGCCAGGTAGGTGTCGCAGGTCATGGCTTCCAGCTCGACGCGGTGGTCAAGGAAGAACTGGCACAGCTCGTTGATGATGTCCGGGCGGTAGGCAGCGCTCACGTAGGCCACATACGGAAGGGCCTGCGGGCGCACTTCCTGGTCCGCGCTGCGGACCACGTCGAGGGTCAGGCCGTGCTTCTTGCCCAAGCCTGGCAGGGTGGATTCGAGACGTGCCAGGGCGTCCCAGCTGCCGCCCACCTGCAGCACCAGGGCGCTGGTCTCGCCGTGGCGGCTCAGGCGCGAGGTGACCA  >  NC_002947/1414068‑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: