Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,093,139	0	C	A	28.9%	65.6 / 11.0	38	V100V (GTG→GTT)	rhaM	L‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base C (10/17);  new base A (11/0);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 3.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGC  >  NC_000913/4092932‑4093361                                                                                                                                                                                                                |                                                                                                                                                                                                                                    cAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTTAAAAACTTCTTTCAACTCg                                                                                                                                                                                                                      >  1:63827/1‑223 (MQ=255)                        aaTGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTGAAAAACTTCTTGCAGCTCGCCACT‑CACACGGCTGTTTTCCCCGTTAg                                                                                                                                                                                         >  2:447263/1‑228 (MQ=255)                                                                aTTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTTCAGCTCGCCACT‑CACCCGGCTGTGATCACGGTTAGCGGGGATAAAATCGGTCATAGA‑TTTCCACCAC‑AGTTGGcc                                                                                                                                               >  1:231651/1‑227 (MQ=255)                                                                aTTCATTATTTTTGTGATTTACTCACAAGGTATAAGGGAAATATAACTTACCGGGAGATGTAATGTGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCa                                                                                                                                               <  2:635160/228‑1 (MQ=255)                                                                                               aTCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTTAAAAACCTCTTTCAACCTGCTACT‑CACCCGTCTGTTATTCCGGTTATCGGGGCTAACATCGGGCTCATA‑TTTCTCCCCA‑CAATGGCAAAAATCA‑TGTCTGTCAAACGCCTTcccacgt                                                                                                               >  2:651894/1‑223 (MQ=255)                                                                                                   gggAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑Gttctt                                                                                                          <  1:447996/228‑1 (MQ=255)                                                                                                      aaaTATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCCCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAAAATCGGGCATATA‑TTTTCACCAACCGTTGGCACAAACCA‑‑TGC‑GGCAACCGCATTCCaagcgtcctcag                                                                                                        >  2:277399/1‑216 (MQ=255)                                                                                                        atatAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGc                                                                                                                                                                                                                            >  2:405843/1‑114 (MQ=255)                                                                                                        atatAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAAATTCTTGCAGGTCGCTAAT‑CACCGGGCTGTTATCCGGGTTAGCGGGGATAACATAGGTCACATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGGAAACGCATTCCAAA‑GTTTTTCTGa                                                                                                      >  1:19642/1‑227 (MQ=255)                                                                                                        atatAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACT‑CACCGGGCt                                                                                                                                                                                                       >  2:124830/1‑134 (MQ=255)                                                                                                         tataACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCa                                                                                                                                                                                                                           <  1:406505/114‑1 (MQ=255)                                                                                                         tataACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTg                                                                                                                                                                                                      <  1:125013/134‑1 (MQ=255)                                                                                                              cTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTTCAGCTCGCTACTCCCCCGGCCTGTATTCCGGG‑TAGCGGGCATCACATCGGTCATATA‑TTTTCACCCA‑CACTTGGAAACATTA‑ATGCTGGCAAACGCATTTCAGC‑CTTCTTCAGAt                                                                                                     >  2:679930/1‑222 (MQ=255)                                                                                                               ttACCGGGATATGTAATAGGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGAtt                                                                                                    <  1:681043/222‑1 (MQ=255)                                                                                                                tACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATGTTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTGCAGCTCGCTACT‑CAc                                                                                                                                                                                                             >  1:145501/1‑120 (MQ=255)                                                                                                                 aCCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATGTTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CAcc                                                                                                                                                                                                            <  2:145290/120‑1 (MQ=255)                                                                                                                 aCCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTCTTTCAGCTCCCTACT‑CACCCGGCTGTTTTCCCGGGTAGCGGGCATAACAACGGTCATATA‑TTTCCACCAC‑AGTTGGGAAACATCA‑ATGCTGGCAACCGCCATCCAGA‑GTGCTTCAGAATCAATatct                                                                                            >  2:648551/1‑228 (MQ=255)                                                                                                                     ggAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGGAAAAAACTTCTTGCAGCTCGCTACT‑CACCCGGCTGTTATCCGGGGTAGCgg                                                                                                                                                                                      >  2:634728/1‑138 (MQ=255)                                                                                                                      gagaTGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCggg                                                                                                                                                                                     <  1:635739/138‑1 (MQ=255)                                                                                                                                atGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTCAAACAATTCTTTTAGCTCGCCACT‑CACCGGGGAGTTAt                                                                                                                                                                                                  >  2:13266/1‑115 (MQ=255)                                                                                                                                 tGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATc                                                                                                                                                                                                 <  1:13286/115‑1 (MQ=255)                                                                                                                                  gTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAAAACTTTTTTCAGCCCGCTACT‑CACCGGGGTGTTTTCCCGGGTAGCGGGCATAACATCGGGCCTATATTTTCCACCAC‑GCTTGGCAAAC‑TCAGGGGCTGG‑AAACGCATTCCAACGGTATTTCAGTTTC‑ACATCTATCATGGCAAAcagcag                                                                           >  2:639581/1‑228 (MQ=255)                                                                                                                                   tttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAg                                                                                                                                                                                         >  1:389260/1‑121 (MQ=255)                                                                                                                                    ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTTTTGCAGCTCGCCACT‑CTACCGGCTGTTATTCCGGTTAACGGGCCTAACATCGGTCATATA‑TTTTCACCAA‑CCTTTGCAAACATCA‑GTGCCGGCAACCGACTTCCTGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAAcag                                                                              >  2:84964/1‑223 (MQ=255)                                                                                                                                    ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTAc                                                                                                                                                                                                                  >  2:363512/1‑96 (MQ=255)                                                                                                                                    ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGc                                                                                                                                                                                        <  2:388619/121‑1 (MQ=255)                                                                                                                                    ttAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCCGGCTGGTATTCGGGGTAGCGGGCATAACATCGGGCATATA‑TTTTCACCAA‑CGTTGGCAAACATCA‑ATGCTTGCAACCGCCTTCCATC‑GTTCTTCATATTCAATCTCTACCATGGccacaagc                                                                             >  2:327137/1‑217 (MQ=255)                                                                                                                                     tAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACt                                                                                                                                                                                                                 <  1:364116/96‑1 (MQ=255)                                                                                                                                                                              tctcTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGt                               <  2:63735/228‑1 (MQ=255)                                                                                                                                                                                 ctGGCACATTGGGCAATTCCGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTga                             <  2:231295/227‑1 (MQ=255)                                                                                                                                                                                            ggCAATTACGGCAGGTAAAACACTTCTTGCAGCCCCCTACT‑CACCCGGCTGTTTTCCGGGGTAGCGGGCATAAAAACGGTCATATA‑TTTCCCCCCA‑CGTTGGCAAACATCA‑ATGCTGGCAAACGCATTCCACC‑GTTCTTTCGATTTAAACTCTACCAAGGCAAACTACAGGATACGCGCTTCGTCGAGATAGATGGc                                                >  2:385959/1‑197 (MQ=255)                                                                                                                                                                                             gCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCg                                               <  1:386599/197‑1 (MQ=255)                                                                                                                                                                                               aaTTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCtt                                                                <  1:398274/178‑1 (MQ=255)                                                                                                                                                                                                 ttACGGCAGGGAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCCCTTTGTCGAGATAGATGGCGTAGTTATTCGGAACGTg                              >  1:403947/1‑210 (MQ=255)                                                                                                                                                                                                  tACGGCAGGTAAACCACTTCTTGCAGCTCGCTACT‑CACCGTGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTga                             <  2:403289/210‑1 (MQ=255)                                                                                                                                                                                                  tACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTAcgcg                                                                   >  1:4930/1‑172 (MQ=255)                                                                                                                                                                                                   aCGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACgcgc                                                                  <  2:4924/172‑1 (MQ=255)                                                                                                                                                                                                           tAAAACACTTCTTGCAGCTCGCTCCT‑CACCGGGTTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGc  <  1:239555/228‑1 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                                    CAGATTTGCGGAGTAATGCACATAATGGTTATTTAAATAAACCACATGAATCATTAATGGGTTATTCATTATTTTTGTGATTTACTCACAATGTATCAGGGAAATATAACTTACCGGGAGATGTAATATGTTTAATTTTTCAAATCGAAATTTAAAATATTGTGCCGGAGGCATCTCTGGCACATTGGGCAATTACGGCAGGTAAAACACTTCTTGCAGCTCGCTACT‑CACCGGGCTGTTATCCGGGTTAGCGGGCATAACATCGGTCATATA‑TTTCCACCAA‑CGTTGGCAAACATCA‑GTGCTGGCAACCGCATTCCAGC‑GTTCTTCAGATTCAATCTCTACCATGGCAAACAGCAGATTACGCGCTTTGTCGAGATAGATGGCGTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGC  >  NC_000913/4092932‑4093361

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