Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,614,048 A→G 7.3% 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,0480AG 7.3% 239.8 / 6.6 83intergenic (+30/+52)PP_4933/hldEhypothetical protein/bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (31/45);  new base G (3/3);  total (35/48)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.51e-01

GGTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATC  >  NC_002947/5613904‑5614195
                                                                                                                                                |                                                                                                                                                   
ggTTAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                 <  7:332171/149‑1 (MQ=255)
ggTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCgg                                                                                                                                                 <  5:332036/149‑1 (MQ=255)
 gTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCggg                                                                                                                                                <  1:348492/149‑1 (MQ=255)
 gTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTAGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGACCggg                                                                                                                                                >  4:365920/1‑149 (MQ=255)
  tGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCTCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGt                                                                                                                                               <  6:266835/149‑1 (MQ=255)
  tGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGGGCGGGt                                                                                                                                               >  7:153706/1‑149 (MQ=255)
  tGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGt                                                                                                                                               <  3:114311/149‑1 (MQ=255)
      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                                                                                                                              <  2:316795/141‑1 (MQ=255)
                           caagcCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCaga                                                                                                                              >  1:316795/1‑141 (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                                                                                              >  4:181264/1‑149 (MQ=255)
                                                   gCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAg                                                                                              >  5:261881/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         >  1:80191/1‑149 (MQ=255)
                                                                                                                                        cttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa         >  7:246916/1‑149 (MQ=255)
                                                                                                                                        cttgctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGa         <  3:295301/149‑1 (MQ=255)
                                                                                                                                            ctGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGa      >  6:262295/1‑148 (MQ=255)
                                                                                                                                               gAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATc  <  3:181264/149‑1 (MQ=255)
                                                                                                                                               gAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATc  <  6:261881/149‑1 (MQ=255)
                                                                                                                                               gAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATc  <  1:104655/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGTGAAGCAGCAGGGGTTGGTGGAAAGCAAGCCCGCACGGGCGTCGCACTTGCAACGCTTTGGCGACAAGGCGCGAGCGCTGCTGGAGCGGGCGCGGTCCTCGCTGAACAAGTGACATCCCGAAAGCTTTCACCTGCTTGCTGGAGCGGGTAAACCCGCCCCAGCAGAGGGCCGCCGAAGCGACCCCAAGGTCAGATCAGTTCTTGCGAATTTTCTCGACAATCGCCGTGGTCGAGCTGTTCTCGACCAGCCCCAGCACCTTCACGGTGCCGCCATAGCCCTTGACGATATC  >  NC_002947/5613904‑5614195

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: