Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,614,055 A→T 6.5% intergenic (+37/+45) 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,0550AT 6.5% 248.6 / 2.5 82intergenic (+37/+45)PP_4933/hldEhypothetical protein/bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/43);  new base T (3/3);  total (36/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 = 1.00e+00

GCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATCAGCA  >  NC_002947/5613910‑5614199
                                                                                                                                                 |                                                                                                                                                
gcagcaGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAATGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAAc                                                                                                                                               <  3:35338/149‑1 (MQ=255)
    cagGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCAGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGc                                                                                                                                           <  6:333086/149‑1 (MQ=255)
       gggTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGctc                                                                                                                                         >  1:252574/1‑146 (MQ=255)
        ggTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCca                                                                                                                                       <  5:71313/149‑1 (MQ=255)
         gTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGCTCcag                                                                                                                                      <  8:153706/149‑1 (MQ=255)
             gtggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGCTCcagcaag                                                                                                                                  <  2:252574/149‑3 (MQ=255)
             gtggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                                  <  3:174052/149‑1 (MQ=255)
             gtggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCACGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                                  <  2:357104/149‑1 (MQ=255)
             gtcgAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGTTTACCCGCTCcagcaag                                                                                                                                  <  5:5621/146‑3 (MQ=255)
              tggAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCagag                                                                                                                                 >  2:14369/1‑149 (MQ=255)
                gAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAggg                                                                                                                               >  7:13139/1‑149 (MQ=255)
                gAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAgg                                                                                                                                <  5:21781/148‑1 (MQ=255)
                  aagcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                                  >  5:300503/1‑144 (MQ=255)
                  aagcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                                  <  6:300503/144‑1 (MQ=255)
                  aagcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgcc                                                                                                                             >  8:15806/1‑149 (MQ=255)
                   agcaagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgccg                                                                                                                            >  4:315759/1‑149 (MQ=255)
                     caagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                                  >  1:316795/1‑141 (MQ=255)
                     caagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                                  <  2:316795/141‑1 (MQ=255)
                      aagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTCCTGGAGCGGGTAAACCCGCCCCAGCAGAGGgccgcc                                                                                                                          >  6:1539/1‑148 (MQ=255)
                        gcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGaa                                                                                                                       <  5:1539/149‑1 (MQ=255)
                         cccGCACGTGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAg                                                                                                                      <  4:179178/149‑1 (MQ=255)
                         cccGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCGGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAg                                                                                                                      <  4:173362/149‑1 (MQ=255)
                          ccGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGc                                                                                                                     <  3:303933/149‑1 (MQ=255)
                              aCGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGAcc                                                                                                                 >  2:258749/1‑149 (MQ=255)
                                  gCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGAcccc                                                                                                               >  3:48106/1‑147 (MQ=255)
                                    gTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGCGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAgg                                                                                                           >  3:79344/1‑149 (MQ=255)
                                       gCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGt                                                                                                          >  5:207311/1‑147 (MQ=255)
                                          cTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGAt                                                                                                     <  5:376894/149‑1 (MQ=255)
                                           ttGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGAt                                                                                                     >  4:350415/1‑148 (MQ=255)
                                             gCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAg                                                                                                  >  5:261881/1‑149 (MQ=255)
                                             gCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAg                                                                                                  >  4:181264/1‑149 (MQ=255)
                                             gCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAg                                                                                                  >  4:223880/1‑149 (MQ=255)
                                              cAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGt                                                                                                 >  2:158813/1‑149 (MQ=255)
                                              cAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGt                                                                                                 >  4:160364/1‑149 (MQ=255)
                                                 cGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCt                                                                                              <  3:304213/149‑1 (MQ=255)
                                                    tttGGCTACAAGGAGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCACAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGc                                                                                           <  7:251501/149‑1 (MQ=255)
                                                          gACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAAttt                                                                                     >  7:319255/1‑149 (MQ=255)
                                                              aGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTtctc                                                                                 <  8:88092/149‑1 (MQ=255)
                                                               ggCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCg                                                                                <  4:307827/149‑1 (MQ=255)
                                                                    gAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAAt                                                                           <  2:313463/149‑1 (MQ=255)
                                                                          tgctgGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGt                                                                     >  6:292718/1‑149 (MQ=255)
                                                                                  cGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAg                                                              <  8:263361/148‑1 (MQ=255)
                                                                                       gcgGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTtc                                                        <  3:223880/149‑1 (MQ=255)
                                                                                       gcgGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTtc                                                        <  5:53617/149‑1 (MQ=255)
                                                                                           tCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGa                                                    >  3:17072/1‑149 (MQ=255)
                                                                                            ccTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTAGAc                                                   >  3:269846/1‑149 (MQ=255)
                                                                                             cTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAATCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGAc                                                   >  6:9463/1‑148 (MQ=255)
                                                                                             cTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGAc                                                   <  6:207311/148‑1 (MQ=255)
                                                                                              tCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGAc                                                   >  5:357081/1‑147 (MQ=255)
                                                                                                gCTGAACAAGTGCCATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCa                                                 <  7:15806/147‑1 (MQ=255)
                                                                                                    aaCAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCa                                           <  4:48106/149‑1 (MQ=255)
                                                                                                      cAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGc                                         <  1:158813/149‑1 (MQ=255)
                                                                                                             cATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTTTTCTATCTTTTCTCGACCAGCCCCATCACCTTCa                                  <  6:240004/149‑1 (MQ=255)
                                                                                                              aTCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTtc                                                        <  3:34723/126‑1 (MQ=255)
                                                                                                              aTCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTtc                                                        >  4:34723/1‑126 (MQ=255)
                                                                                                              aTCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCAc                                 <  6:357081/149‑1 (MQ=255)
                                                                                                               tCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACg                                <  8:319255/149‑1 (MQ=255)
                                                                                                               tCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACg                                <  3:160364/149‑1 (MQ=255)
                                                                                                                cccGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACgg                               >  4:351575/1‑149 (MQ=255)
                                                                                                                 ccGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGt                              >  5:138239/1‑149 (MQ=255)
                                                                                                                     aaGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTgccg                          <  3:350415/149‑1 (MQ=255)
                                                                                                                     aaGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCACCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGCGCTGTTCTCGACCAGCCCCCGCACCTTAACTGTgccg                          >  6:337651/1‑149 (MQ=255)
                                                                                                                        cTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCGTGAAGGTgctgggg                       <  4:17072/149‑6 (MQ=255)
                                                                                                                         tttCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTTGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCa                       >  5:29839/1‑148 (MQ=255)
                                                                                                                            cACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAg                    <  1:258749/148‑1 (MQ=255)
                                                                                                                             aCCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGc                   <  1:14369/148‑1 (MQ=255)
                                                                                                                                  tctgctGGGGCGGGTTTACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa             >  1:358034/3‑149 (MQ=255)
                                                                                                                                  cttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTg              <  3:315759/148‑1 (MQ=255)
                                                                                                                                  cttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa             <  5:292718/149‑1 (MQ=255)
                                                                                                                                  cttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa             >  7:246916/1‑149 (MQ=255)
                                                                                                                                  cttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa             <  3:295301/149‑1 (MQ=255)
                                                                                                                                  cttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa             >  1:80191/1‑149 (MQ=255)
                                                                                                                                      ctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGa          >  6:262295/1‑148 (MQ=255)
                                                                                                                                        ggggCGGGTTTACCCGCTCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGatat       >  5:285093/4‑149 (MQ=255)
                                                                                                                                         gAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATc      <  6:261881/149‑1 (MQ=255)
                                                                                                                                         gAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATc      <  3:181264/149‑1 (MQ=255)
                                                                                                                                         gAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATc      <  1:104655/149‑1 (MQ=255)
                                                                                                                                           gCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCa                                                 <  7:91118/104‑1 (MQ=255)
                                                                                                                                           gCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCa                                                 >  8:91118/1‑104 (MQ=255)
                                                                                                                                             gggTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTTGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATCAGca  <  6:29839/149‑1 (MQ=255)
                                                                                                                                             gggTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCTTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATCAGca  <  5:9463/149‑1 (MQ=255)
                                                                                                                                             gggTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATCAGca  >  2:214333/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATCAGCA  >  NC_002947/5613910‑5614199

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