Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,270,4550GT40.8% 69.1 / 19.3 49intergenic (+853/+170)yhaC/garKpentapetide repeats‑related protein/glycerate kinase I
Reads supporting (aligned to +/- strand):  ref base G (6/23);  new base T (20/0);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.63e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGC  >  NC_000913/3270305‑3270599
                                                                                                                                                      |                                                                                                                                                   
gTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTg                                                                                                                                                     <  1:320780/151‑1 (MQ=255)
  tCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGc                                                                                                                                                    <  1:455553/150‑1 (MQ=255)
   cAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGctc                                                                                                                                                  <  1:885008/151‑1 (MQ=255)
         gggCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGctctctct                                                                                                                                             >  1:459277/1‑148 (MQ=255)
         gggCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGCGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGt                                                                                                                                             >  1:253483/1‑150 (MQ=255)
          ggCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTg                                                                                                                                           <  1:421057/151‑1 (MQ=255)
                aCCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTGCCGCACCCTTTCGCCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTgg                                                                                                                                     <  2:614343/151‑1 (MQ=255)
                aCCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGGTTTTTCTCTGTTGCAc                                                                                                                                        >  2:855745/1‑148 (MQ=255)
                 ccTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGt                                                                                                                                    <  2:724796/151‑1 (MQ=255)
                   ttATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCAc                                                                                                                                        >  2:768047/1‑145 (MQ=255)
                   ttATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTGTGccgtctttttga                                                                                                                                         >  1:636160/1‑133 (MQ=255)
                       gAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACCGGG‑CCT‑ggg                                                                                                                            >  1:545497/1‑151 (MQ=255)
                        aaCCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGT‑CGT‑gg                                                                                                                             >  2:68356/1‑149 (MQ=255)
                           cccTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTc                                                                                                                        <  2:934038/151‑1 (MQ=255)
                                tttGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGT‑CGT‑GGGGTTcccccc                                                                                                                   >  1:394581/1‑151 (MQ=255)
                                 ttGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCCCCCTTCCCATGCTCCGG‑GGGTTTcccccc                                                                                                                   >  1:421573/1‑151 (MQ=255)
                                   ggCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGT‑CGT‑GGGGTTccccccc                                                                                                                  >  2:517866/1‑149 (MQ=255)
                                   ggCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCa                                                                                                                 <  1:903941/150‑1 (MQ=255)
                                    gCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAgg                                                                                                               <  1:401204/151‑1 (MQ=255)
                                     ccTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGc                                                                                                              <  1:882841/151‑1 (MQ=255)
                                      cTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTTGT‑CGG‑GGGGTTCCCCCCCAgg                                                                                                               >  2:349857/1‑149 (MQ=255)
                                      cTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGCTTTCTTCCCGTTGCTCTCTG‑GGT‑CGGGTTCCCCCCCCGGc                                                                                                              >  1:188060/1‑150 (MQ=255)
                                        tGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTTTCTCTCTTTCACTGGT‑CGG‑GGGGTTTCCCCCCCGGCGtt                                                                                                           >  2:760574/1‑151 (MQ=255)
                                        tGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCg                                                                                                             <  1:216092/149‑1 (MQ=255)
                                         gCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACCTGT‑CGG‑GGGGTTCCCCCCCCgg                                                                                                               >  2:468176/1‑146 (MQ=255)
                                         gCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTa                                                                                                          <  2:642077/151‑1 (MQ=255)
                                            ccGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTAcc                                                                                                        <  1:108846/150‑1 (MQ=255)
                                             cGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTg                                                                                                      <  1:81462/151‑1 (MQ=255)
                                              gggTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTgg                                                                                                     <  2:306902/151‑1 (MQ=255)
                                                       ttACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTGCCCCTCTTGCACCGGG‑CCT‑GGGGTTCCCCCCCCAGGGTTAACCTGCAccccccc                                                                                            >  1:424328/1‑147 (MQ=255)
                                                               ccACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAg                                                                                    <  1:749356/151‑1 (MQ=255)
                                                                       tGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGT‑CGT‑GGGGTTTCCCCCCAGGCGTTTCCCGGCACCCCGCCCTATGGAGCCCGGGCt                                                                            >  1:811703/1‑151 (MQ=255)
                                                                       tGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGGTTTCCCCCCCGGCGGTACCTGGCACCCTGCCCTCTGGAGCCCGGGCt                                                                            >  2:314636/1‑151 (MQ=255)
                                                                                 cgcgcgGTGCGCTCTTACCGCCCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTcct‑cccct                                                                 <  1:936927/151‑1 (MQ=255)
                                                                                    gcgGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGG‑CGG‑GGGGTTCCCCCCCCGGCGCTTCCTGGCACCCCGCCCTATGGAGCCCGGGCTTTcct‑cccctcc                                                               >  2:952875/1‑150 (MQ=255)
                                                                                      ggTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGcc                                                            <  2:575262/151‑1 (MQ=255)
                                                                                               ttACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTTCACCGGG‑CGT‑GGGGTTTCCCCCCCGGCGTTACCTCGCACCCCCCCCCATGGGGGCCCGGCTTTCCT‑CCCCCCCCCCCGTCTcccc                                                   >  1:957577/1‑151 (MQ=255)
                                                                                                     cACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCGGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAaga                                             <  1:157767/151‑1 (MQ=255)
                                                                                                      aCCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAagag                                            <  1:540630/151‑1 (MQ=255)
                                                                                                        ccTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTGGT‑CGT‑GGGGTTTCCCCCCCGGCGGTTCCTCGCACCCCGCCCTCTGGAGCCCGGGCTTTTCTCCCCCCCCCCCCGCCCCCCCCGagggg                                            >  1:608838/1‑146 (MQ=255)
                                                                                                          tttCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACCGGG‑CCT‑GGGTTTTCCCCCCCGGCGGTTCCTGGCGCCCCCGCCCCTGGAGGCCGGCCTTTCCT‑CCCCCCCGCCCCTCCCCCCCCAAGGGgacg                                        >  1:802054/1‑151 (MQ=255)
                                                                                                                    aCCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTTCTCTCTTTTGCACTTGG‑CGT‑GGGTTTCCCCCCCCGGCGGTACCCGGCCCCCCGCCCCATGGGGGCCCGGCTTTTCT‑CCCCCCCCCCCCGCCCCCCCCCAGAGgacgacgacca                                 >  1:509040/1‑146 (MQ=255)
                                                                                                                       tGATCCCGCTTGCGCGGGCCATCGGCGGTTTGTTCTCTGTTGCACCTGT‑CGT‑GGGTTTCCCCCCCAAGGGTTTCCTGGCCCCCTGCCCCATAGAGCCCCGACCTTCCT‑CCCCTCCCCCCGCCTTCCCCCAAGAGgacgaagcacaaacg                             >  2:464210/1‑139 (MQ=255)
                                                                                                                           cccGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACt                       <  2:291415/151‑1 (MQ=255)
                                                                                                                              gCTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTGGG‑CGT‑GGGGTTCCCCCCCCGGCGGTTCCCGGCACCCCGCCCCCTTGGGCCCGGGCCTTTCT‑CCCCCCCCCCCCCCTCCCCCCCAGAGGACGACGACGAAgcgg                            >  1:546734/1‑143 (MQ=255)
                                                                                                                               cTTGCGCGGGCCATCGGCGGTTTTCTCTCTGTTGCACTTGG‑CGTGGGGTTTCCCCCCCCGG‑GTTACCTGGCACCCCGCCCCCTGGAGGCCCGGCTTTCCT‑CCCCCCCCCCCCCCTCCCCCCCAGAGGGCGACGACCAAgcggcgcc                        >  1:398563/1‑144 (MQ=255)
                                                                                                                                    gcgGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCa              <  2:41490/151‑1 (MQ=255)
                                                                                                                                             cggcggTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgc      <  1:695591/150‑1 (MQ=255)
                                                                                                                                               gcggTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcgcg   >  2:646521/1‑151 (MQ=255)
                                                                                                                                                cggTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGgcgcgc  <  2:291198/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GTTCAGTACGGGCCGTACCTTATGAACCCCTATTTGGCCTTGCTCCGGGTGGAGTTTACCGTGCCACGGACTGTTACCAGCCGCGCGGTGCGCTCTTACCGCACCCTTTCACCCTTACCTGATCCCGCTTGCGCGGGCCATCGGCGGTTTGCTCTCTGTTGCACTGGT‑CGT‑GGGTTTCCCCCCCAGGCGTTACCTGGCACCCTGCCCTATGGAGCCCGGACTTTCCT‑CCCCTCCGCCCGTCTCCCCCGAAGAGGACGACGACGAAGCGGCGACTGTCTGGTCAGCTTCGGCGCGC  >  NC_000913/3270305‑3270599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑