Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,614,099 A→G 7.8% intergenic (+30/+52) PP_4933 → / ← hldE hypothetical protein/bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,614,0990AG 7.8% 253.1 / 8.8 90intergenic (+30/+52)PP_4933/hldEhypothetical protein/bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (41/42);  new base G (3/4);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

TTGGGTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATAT  >  NC_002947/5613952‑5614245
                                                                                                                                                   |                                                                                                                                                  
ttGGGTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGc                                                                                                                                                  <  2:2332183/150‑1 (MQ=255)
  gggTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                <  2:1867817/150‑1 (MQ=255)
  gggTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                <  1:1119442/150‑1 (MQ=255)
   ggTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCggg                                                                                                                                               >  1:1534329/1‑150 (MQ=255)
   ggTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCggg                                                                                                                                               <  2:1475536/150‑1 (MQ=255)
      gAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCggg                                                                                                                                               >  1:1200545/1‑147 (MQ=255)
      gAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCggg                                                                                                                                               <  2:1200545/147‑1 (MQ=255)
      gAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTaa                                                                                                                                            >  1:1733649/1‑150 (MQ=255)
        agcagcAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGtttac                                                                                                                                          <  1:1172482/150‑5 (MQ=255)
            gcagGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGc                                                                                                                                      >  2:2154943/1‑150 (MQ=255)
            gcagGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGc                                                                                                                                      >  1:1934700/1‑150 (MQ=255)
            gcagGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGc                                                                                                                                      <  2:2339596/150‑1 (MQ=255)
              agGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGccc                                                                                                                                    <  1:2177982/150‑1 (MQ=255)
                gggTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGCTCca                                                                                                                                  <  1:2154943/150‑1 (MQ=255)
                gggTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCca                                                                                                                                  <  2:1733649/150‑1 (MQ=255)
                gggTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCca                                                                                                                                  <  1:657262/150‑1 (MQ=255)
                     tctgggaagcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCGCCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                             >  2:343206/7‑150 (MQ=255)
                       tggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGCTCcagcaagca                                                                                                                           <  2:103343/150‑5 (MQ=255)
                       tggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCggg                                                                                                                                               >  2:2230214/1‑130 (MQ=255)
                       tggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCggg                                                                                                                                               <  1:2230214/130‑1 (MQ=255)
                       tggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAgg                                                                                                                           <  1:1292028/150‑1 (MQ=255)
                       tggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAgg                                                                                                                           >  2:690115/1‑150 (MQ=255)
                         gAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGCTCcagcaagcag                                                                                                                          >  1:1890321/1‑144 (MQ=255)
                         gAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGCTCcagcaagcag                                                                                                                          <  2:1890321/149‑6 (MQ=255)
                         gAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgc                                                                                                                         >  2:109644/1‑150 (MQ=255)
                           aagcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgccg                                                                                                                       >  1:842639/1‑150 (MQ=255)
                           aagcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgccg                                                                                                                       >  2:1478095/1‑150 (MQ=255)
                            agcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCCCCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgccgc                                                                                                                      >  1:873153/1‑150 (MQ=255)
                             gcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgccgcc                                                                                                                     <  1:394235/150‑1 (MQ=255)
                                agcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGaa                                                                                                                  <  2:2419491/150‑1 (MQ=255)
                                agcCCGCACGGGCGTCGCACTTGCAACGCTGTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACACCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGaa                                                                                                                  <  1:1086352/150‑1 (MQ=255)
                                 gcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAg                                                                                                                 >  2:1577493/1‑150 (MQ=255)
                                  cccGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAACAAATGCCGACTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGc                                                                                                                >  1:1970133/1‑150 (MQ=255)
                                           gCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAACACCGAACCAGACGAGGGCAAGCGAAGCGACCCCAAg                                                                                                       >  2:1719359/1‑150 (MQ=255)
                                           gCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAg                                                                                                       <  1:1006959/150‑1 (MQ=255)
                                             gTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                >  1:2372876/1‑107 (MQ=255)
                                             gTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                <  2:2372876/107‑1 (MQ=255)
                                             gTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGt                                                                                                     >  1:1277408/1‑150 (MQ=255)
                                                  aCTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAggg                                                                                                                          >  2:137968/1‑124 (MQ=255)
                                                  aCTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAggg                                                                                                                          <  1:137968/124‑1 (MQ=255)
                                                  aCTTCCAACGCTTTGTCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGAt                                                                                                <  1:377505/150‑1 (MQ=255)
                                                    ttGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCa                                                                                              <  2:329111/150‑1 (MQ=255)
                                                     tGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAg                                                                                             <  1:408112/150‑1 (MQ=255)
                                                     tGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAg                                                                                             <  1:109644/150‑1 (MQ=255)
                                                       cAACGCTTTGGCGACAAGGCGCGAACGCTGCTGGAACGAGCGCAGACCTCGGTGAACAACGGATAACTCAAAAACAATCACCTGCTTGCTGGAACGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGtt                                                                                           >  2:2102588/1‑150 (MQ=255)
                                                        aaCGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                <  1:1472989/96‑1 (MQ=255)
                                                        aaCGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                >  2:1472989/1‑96 (MQ=255)
                                                        aaCGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGa                                                                                                 >  2:968947/1‑143 (MQ=255)
                                                        aaCGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGa                                                                                                 <  1:968947/143‑1 (MQ=255)
                                                        aaCGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTc                                                                                          >  1:1820002/1‑150 (MQ=255)
                                                         aCGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCt                                                                                         >  2:2059966/1‑150 (MQ=255)
                                                          cGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCtt                                                                                        <  2:1534329/150‑1 (MQ=255)
                                                             tttGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                <  2:266701/91‑1 (MQ=255)
                                                             tttGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                >  1:266701/1‑91 (MQ=255)
                                                                 gCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAAtt                                                                                 >  1:1034876/1‑150 (MQ=255)
                                                                                   tgctgGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTg                                                               >  1:1389309/1‑150 (MQ=255)
                                                                                     ctgGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTACTCGACAATCGCCGTGGt                                                             <  1:54939/150‑1 (MQ=255)
                                                                                        gAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTCTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGa                                                          <  1:1719359/150‑1 (MQ=255)
                                                                                         aGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAg                                                         <  2:1389309/150‑1 (MQ=255)
                                                                                          gCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGc                                                        <  1:513268/150‑1 (MQ=255)
                                                                                            gggCGCGGTCCTCGCTGAACAAGTGAAAAGACGATAACGTTCACCTGCTTGCTGGAGCCCGTCAACCCGCCCAAGCAGTCATGGGCCGAAATGACGGCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCTAGGTCTAGCTg                                                      >  1:2388751/1‑150 (MQ=255)
                                                                                             ggCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGt                                                     >  2:1180332/1‑150 (MQ=255)
                                                                                                gcgGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTtctc                                                                            <  2:1402052/124‑1 (MQ=255)
                                                                                                gcgGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTtctc                                                                            >  1:1402052/1‑124 (MQ=255)
                                                                                                gcgGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTtct                                                  <  1:458638/150‑1 (MQ=255)
                                                                                                 cgGTCCTCGCTGAACAAGTGACATCCCGAGAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCATTCCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTtctc                                                 <  1:1223296/150‑1 (MQ=255)
                                                                                                 cgGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTtctc                                                 <  1:2416015/150‑1 (MQ=255)
                                                                                                   gTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGa                                               <  2:67004/150‑1 (MQ=255)
                                                                                                   gTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGa                                               >  1:1781675/1‑150 (MQ=255)
                                                                                                    tCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGAc                                              <  1:2059966/150‑1 (MQ=255)
                                                                                                          cTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGt                                                                <  2:483217/126‑1 (MQ=255)
                                                                                                          cTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGt                                                                >  1:483217/1‑126 (MQ=255)
                                                                                                            gAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCa                                      <  2:1034876/150‑1 (MQ=255)
                                                                                                               cAAGTGAAATCCAGAACGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCa                                   >  2:1797747/1‑150 (MQ=255)
                                                                                                                aaGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCAc                                  >  1:872627/1‑150 (MQ=255)
                                                                                                                  gTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCt                                >  2:227783/1‑150 (MQ=255)
                                                                                                                  gTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCt                                >  2:164897/1‑150 (MQ=255)
                                                                                                                  gTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCt                                >  1:790936/1‑150 (MQ=255)
                                                                                                                           cGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTgc                       >  1:1194352/1‑150 (MQ=255)
                                                                                                                                 cTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGt                                                     <  2:295091/114‑1 (MQ=255)
                                                                                                                                 cTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGt                                                     >  1:295091/1‑114 (MQ=255)
                                                                                                                                 cTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCAt                 >  2:2065815/1‑150 (MQ=255)
                                                                                                                                      aCCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGccc            <  1:1180332/150‑1 (MQ=255)
                                                                                                                                        cTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGt                                                                <  2:207816/96‑1 (MQ=255)
                                                                                                                                        cTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGt                                                                >  1:207816/1‑96 (MQ=255)
                                                                                                                                        cTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCtt          >  1:1858968/1‑150 (MQ=255)
                                                                                                                                          gcttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa        <  1:2102588/150‑1 (MQ=255)
                                                                                                                                          ctctgctGGGGCGGGTTTACCCGCTCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa        >  2:1172823/4‑150 (MQ=255)
                                                                                                                                             tgctGGAGCGGGTAACCCCGCACCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCAAATTAACGGTGCCGCCATCGCCCTTGACgc     <  2:1114714/150‑2 (MQ=255)
                                                                                                                                                tGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGatat  >  2:815731/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGGGTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATAT  >  NC_002947/5613952‑5614245

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: