Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913244,1330GA54.9% 158.3 / inf 395intergenic (+12/‑194)gmhA/yafJphosphoheptose isomerase/putative amidotransferase
Reads supporting (aligned to +/- strand):  ref base G (96/81);  new base A (113/102);  total (211/184)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 403 total aligned reads displayed.

GGAGATTCACATTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGC  >  NC_000913/244058‑244208
                                                                           |                                                                           
  agaTTCACATTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGcc                                                                           >  1:1025473/1‑76 (MQ=255)
  agaTTCACATTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGCGGAGGCCg                                                                          >  2:781661/1‑76 (MQ=255)
    aTTCACATTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCgg                                                                         <  2:5114014/76‑1 (MQ=255)
     ttCACATTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGa                                                                        <  1:2941851/76‑1 (MQ=255)
      tCACATTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGAt                                                                       <  1:5587893/76‑1 (MQ=255)
        acaTTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATaa                                                                     <  1:5018855/76‑1 (MQ=255)
            tAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACg                                                                 >  1:1333406/1‑76 (MQ=255)
              taGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGCCCGGATAAGACGtt                                                               >  2:282415/2‑76 (MQ=255)
               aGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGctt                                                              >  2:1115755/1‑76 (MQ=255)
                gTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTa                                                             >  1:4998410/1‑76 (MQ=255)
                gTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGAAGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTa                                                             <  2:3257220/76‑1 (MQ=255)
                 tGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTAc                                                            <  2:1529706/76‑1 (MQ=255)
                   aTCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTAcgc                                                          <  2:4125839/76‑1 (MQ=255)
                    tCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTAcgcc                                                         >  1:4221615/1‑76 (MQ=255)
                      cATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTAcgccgc                                                       <  2:2966121/76‑1 (MQ=255)
                       atatCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCa                                                      <  1:4850500/76‑1 (MQ=255)
                         atCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATc                                                    <  1:4369358/76‑1 (MQ=255)
                            cTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCgg                                                 <  1:231681/76‑1 (MQ=255)
                              gATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCa                                               >  2:4404245/1‑76 (MQ=255)
                               aTCCAGTTGATTGAAAAAGAGATGGTTAAGCAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCAt                                              >  1:1216377/1‑76 (MQ=255)
                                tatgcgcgtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCAtt                                             >  1:5420391‑M1/40‑76 (MQ=255)
                                tatgcgcgtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCAtt                                             <  2:1867994‑M1/37‑1 (MQ=255)
                                   aGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCCGCATTtgg                                          >  2:1777678/1‑75 (MQ=255)
                                    gTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgtg                                         >  2:3178780/1‑76 (MQ=255)
                                    gTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTtgtg                                         >  2:3684467/1‑76 (MQ=255)
                                     gcgtcttgtcaggtctacggtttatgggcgaagtGTAGGCCGGATAGGCCGTTTACGCCGCATCGGGCATTTGTgc                                        >  1:2587534‑M1/35‑76 (MQ=255)
                                     gcgtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgc                                        <  1:1785966‑M1/42‑1 (MQ=255)
                                      cgtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                       <  1:3943617‑M1/43‑1 (MQ=255)
                                      cgtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcg                                       >  1:2134025‑M1/34‑76 (MQ=255)
                                       gtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgc                                      <  1:3717828‑M1/44‑1 (MQ=255)
                                       gtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgc                                      <  2:3577314‑M1/44‑1 (MQ=255)
                                       gtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgc                                      <  2:5496261‑M1/44‑1 (MQ=255)
                                       gtcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTgcgc                                      >  1:5077695‑M1/33‑76 (MQ=255)
                                        tcttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCCGCATTTGtgcgcg                                     <  2:1313756‑M1/45‑2 (MQ=255)
                                         gtgaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                    <  2:4489849/75‑1 (MQ=255)
                                         cttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                    >  1:1413066‑M1/31‑76 (MQ=255)
                                         cttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                    >  2:3030060‑M1/31‑76 (MQ=255)
                                          ttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                   >  2:1650640‑M1/30‑76 (MQ=255)
                                          ttatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctga                                   <  1:569386‑M1/47‑1 (MQ=255)
                                           tatcaggcctgcggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                  >  2:5347369‑M1/29‑76 (MQ=255)
                                           tatcaggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                  >  2:4290449‑M1/29‑76 (MQ=255)
                                           gaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                  >  2:736674/1‑76 (MQ=255)
                                           gaAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                  <  2:2020854/76‑1 (MQ=255)
                                            aaaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg                                 <  1:3880464/76‑1 (MQ=255)
                                               aggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                              >  1:5234861‑M1/25‑76 (MQ=255)
                                               aggcctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcct                              >  1:5180846‑M1/25‑76 (MQ=255)
                                               aggcatacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgacgcct                              <  2:5788765‑M1/52‑1 (MQ=255)
                                                agagATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctg                             >  2:639436/1‑76 (MQ=255)
                                                   gaTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                          >  1:3379482/1‑76 (MQ=255)
                                                   gaTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                          <  1:3308005/76‑1 (MQ=255)
                                                   ctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                          <  2:3315800‑M1/56‑1 (MQ=255)
                                                   ctacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                          <  2:4829662‑M1/56‑1 (MQ=255)
                                                    tgcggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatgc                         >  2:889972‑M1/20‑76 (MQ=255)
                                                    tacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatgc                         <  1:985542‑M1/57‑1 (MQ=255)
                                                    tacggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatgc                         <  1:1158480‑M1/57‑1 (MQ=255)
                                                    aTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatgc                         <  1:3162068/76‑1 (MQ=255)
                                                     tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                        >  1:4798218/1‑76 (MQ=255)
                                                     acggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCg                        >  1:4512151‑M1/19‑76 (MQ=255)
                                                      ggTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGCTGCCTGATGCGa                       >  2:4643746/1‑76 (MQ=255)
                                                      cggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGa                       >  1:1103972‑M1/18‑76 (MQ=255)
                                                      cggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGa                       >  1:3462139‑M1/18‑76 (MQ=255)
                                                      cggtatatggggtgagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGa                       >  2:2269752‑M1/18‑76 (MQ=255)
                                                       ggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGAc                      <  1:5338090‑M1/60‑1 (MQ=255)
                                                       ggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGAc                      >  1:4120589‑M1/17‑76 (MQ=255)
                                                       ggtttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGAc                      <  1:713787‑M1/60‑1 (MQ=255)
                                                       gTTAAGTAAGTCTGGCGTAGGCCGGATGAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGAc                      >  1:3660096/1‑76 (MQ=255)
                                                       gTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGAc                      <  2:4922183/76‑1 (MQ=255)
                                                       gTTAAGTAAGTCTCGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGAc                      <  2:2216867/76‑1 (MQ=255)
                                                         tttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGc                    <  1:1572878‑M1/62‑1 (MQ=255)
                                                         taagtaagTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGc                    >  1:5828093/1‑76 (MQ=255)
                                                          ttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                   >  2:5562188‑M1/14‑76 (MQ=255)
                                                          ttatgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                   >  1:1623506‑M1/14‑76 (MQ=255)
                                                          aagtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                   <  2:5009359/76‑1 (MQ=255)
                                                           agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTg                  <  1:3065784/76‑1 (MQ=255)
                                                           agtaagtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTg                  <  2:3795701/76‑1 (MQ=255)
                                                            atgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGa                 <  1:1522548‑M1/65‑1 (MQ=255)
                                                            atgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCCGATGCGACGCtgg                 <  2:2675567‑M1/65‑2 (MQ=255)
                                                             tgggcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAc                <  1:3237579‑M1/66‑1 (MQ=255)
                                                              aagtCTGGCGTAGGCCGGATAAGACGTTTCCGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGCCGCTGAcg               >  2:5825054/1‑76 (MQ=255)
                                                               ggcggagtGTGGGCCGGATGAGGCGTTTAGGCCGCATCCGGCATTGGTGCGCTGATGCCTGATGCGACGCTGAcgc              >  2:471736‑M1/9‑76 (MQ=255)
                                                               ggcgaggtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcgc              >  2:3512872‑M1/9‑76 (MQ=255)
                                                                gcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcgcg             <  2:2449415‑M1/69‑1 (MQ=255)
                                                                gcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcgcg             <  1:844123‑M1/69‑1 (MQ=255)
                                                                gcgaagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAcgcg             >  1:5497956‑M1/8‑76 (MQ=255)
                                                                 tCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGt            >  1:802119/1‑76 (MQ=255)
                                                                 tCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGt            <  1:1932850/76‑1 (MQ=255)
                                                                 gcTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGt            <  2:3093288/75‑1 (MQ=255)
                                                                   tGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCt          >  2:254248/1‑76 (MQ=255)
                                                                   aagtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCt          <  2:2134466‑M1/72‑1 (MQ=255)
                                                                    agtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCtt         <  2:3437568‑M1/73‑1 (MQ=255)
                                                                     gtGTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa        <  1:916156‑M1/74‑1 (MQ=255)
                                                                     gCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTa        >  1:4605857/1‑76 (MQ=255)
                                                                       gTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATc      >  2:4350837/1‑76 (MQ=255)
                                                                       gTAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATc      >  2:3847924/1‑76 (MQ=255)
                                                                       gTAGACCGGAGAAGACGTTTACGCCGCATCCGGCATTTGGGCGCTGATGCCTGATCCGACGCGGACGCGTCGTATc      >  2:2165758/1‑76 (MQ=18)
                                                                        tAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCa     >  2:1376583/1‑76 (MQ=255)
                                                                        tAGACCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCa     <  1:2764261/76‑1 (MQ=255)
                                                                         aGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCAt    <  1:2505638/76‑1 (MQ=255)
                                                                         aGGCCAGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTCTCac    >  2:2387578/1‑75 (MQ=255)
                                                                           gCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGc  >  2:5145772/1‑76 (MQ=255)
                                                                           |                                                                           
GGAGATTCACATTAAAGTGATCCATATCCTGATCCAGTTGATTGAAAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGC  >  NC_000913/244058‑244208

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