Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,361,337	0	A	C	61.4%	30.4 / 51.6	45	G26G (GGT→GGG)	rhmD	L‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (11/6);  new base C (0/27);  total (11/34)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-07

GGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATA  >  NC_000913/2361189‑2361480                                                                                                                                                      |                                                                                                                                                ggCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCAc                                                                                                                                                <  2:175226/150‑1 (MQ=255)          cTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAgc                                                                                                                                      >  1:756005/1‑151 (MQ=255)            gCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAgcg                                                                                                                                     <  2:288761/150‑1 (MQ=255)                 gACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCtttt                                                                                                                               <  1:504741/151‑1 (MQ=255)                     tCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                                                                                            <  1:16800/150‑1 (MQ=255)                                         accGTTAATCCCAAAGGCCGGGCGTGACTGCTCGAAATCGCGGTATTTCCCCATCGGGGTGCAAATATGATCGTCTACCCGGTGATT‑CGCCCCCTGGCCGTAATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAc                                                                                                       <  1:390553/149‑1 (MQ=255)                                            gTTAAGGCAAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGAAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAcc                                                                                                      <  2:342762/149‑1 (MQ=255)                                             ttAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAgg                                                                                                   >  2:271952/1‑151 (MQ=255)                                               aaTGCCAAAGGACTGGCGTGACTGCTCGTAATCGGGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CCCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgc                                                                                                  <  2:79005/150‑1 (MQ=255)                                                aggCCAAATGACGGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCGGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcgc                                                                                                <  2:422446/149‑1 (MQ=255)                                                 tGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcgc                                                                                                >  2:421594/1‑150 (MQ=255)                                                   caaaaTGCCTGGCGTAACTGCTCGTAATCGCGGTTTTTACTCATGGGGGTGGAAATATGATCGCCTTTCCAGTGATC‑CCCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGa                                                                                              <  1:549959/148‑1 (MQ=255)                                                     aaaTGCCGGGCGTGACTGCCCGTAATCGGGTTATTTACTCACGGGGGTGGCAATATGATCGTCTTCCCAGTGATT‑CCCCCCCTGGTCGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGa                                                                                              <  1:770107/148‑1 (MQ=255)                                                     aaaTGACTGGCGTGACTGCTCGTAATCGCGGTATTTCCTCATCGGGGTGGCAATATGACCTTCTATCCAGTGATC‑CCCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc                                                                                            <  1:234757/150‑1 (MQ=255)                                                        tGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGt                                                                                        >  2:504459/1‑151 (MQ=255)                                                             ggCGTGACTGCTCTTAACCGCGTTTTTTCCTCATGGGGGTGGCAATATGTTCGTCTATCCGGTGATT‑CGCCCCCTGGCCGTAATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGttt                                                                                      <  1:543269/148‑1 (MQ=255)                                                             ggCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTaa                                                                                    <  2:591955/150‑1 (MQ=255)                                                                    ctgttctgaaTCGCGTTATTTTCTCATGGGGGTGGCAATATGATCGTCTACCCAGTTATT‑CGCCCCCTGGTCGTAATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCCCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg                                                                             <  2:663850/142‑1 (MQ=255)                                                                     tGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGt                                                                            >  1:575568/1‑150 (MQ=255)                                                                      gCTCGTAATCGCGTTATTTACCCATCGGGGTGGCAATATGATCGTCTACCCGGTGATC‑CGCCCCTTGGTCGTGATAACCACCGCCGCCAGCGCCTTTTTCTCCTGTCCCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGt                                                                            <  2:392535/149‑1 (MQ=255)                                                                         cgaaaccgcgGTATTTCCTCATCGGGGGGGCAATATGTTCTTTTACCCAGCGATC‑CGCCCCCTGGTGGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg                                                                              <  1:228236/138‑1 (MQ=255)                                                                         cGTAATCGGGGTATTTACTCATCGGGGTGGCAATATGATCGTCTACCCAGTGATT‑CGCCCCCTGGTCGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa                                                                           <  1:289971/147‑1 (MQ=255)                                                                          gTAATCGCGGTATTTACTCATGGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGAAAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc                                                                      <  2:279838/151‑1 (MQ=255)                                                                          gTAATCGCGGTATTTACTCATCGGGGGGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc                                                                      >  2:209030/1‑151 (MQ=255)                                                                          gTAACCGCGTTATTTACTCACCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAgg                                                                         <  1:880636/148‑1 (MQ=255)                                                                            aaTCGCGGTTTTTCCCCACGGGGGTGGAAATATGTTCGTCTACCCAGTGATC‑CCCCCCCTGTTGGTAATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa                                                                           <  1:970695/144‑1 (MQ=255)                                                                            aaTCGCGGTTTTTACTCATCGGGGTGGAAATATGATCGTCTCTCCAGTGATC‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCa                                                                     <  1:96746/150‑1 (MQ=255)                                                                              tCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatg                                                                   >  1:241836/1‑150 (MQ=255)                                                                                  ggTTTTTACCCATGGGGGGGGCAATATGATCGTCTTCCCAGTTATC‑CCCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGTTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat                                                                 <  1:1014649/148‑1 (MQ=255)                                                                                  ggTATTTACTCATCGGGGTGGAAATAGGATCTTCTATCCAGTGATC‑CGCCCCCTGGTCGTGAAAACCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGtt                                                              <  2:26629/151‑1 (MQ=255)                                                                                    tATTTACTCATCGGGGTGGAAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat                                                                 <  1:896252/146‑1 (MQ=255)                                                                                      tttccccATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTtct                                                            <  2:891740/143‑1 (MQ=255)                                                                                      tttACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTtct                                                            <  2:651015/149‑1 (MQ=255)                                                                                        tccTCATGGGGGTGGAAATATGATCGTCTACCCAGTGTTT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCa                                                         <  2:351956/148‑1 (MQ=255)                                                                                        tACTCATGGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAACCCCCGCCGCCAGCGCCTTTTTCTGTTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCa                                                         <  2:121068/150‑1 (MQ=255)                                                                                                   ggggaaaTATGATCGTTTTCCCAG‑GATTGCGCCCCCTGGCCTTGAAAACCCCCGCCGCCAGCGCTTTTTTCTACTGTCGAACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAgg                                              <  2:157232/145‑1 (MQ=255)                                                                                                          tatGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTa                                        >  1:400704/1‑149 (MQ=255)                                                                                                            tGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAg                                      <  2:799490/149‑1 (MQ=255)                                                                                                                 gTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAgg                               >  2:221748/1‑151 (MQ=255)                                                                                                                       ccAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACttt                         >  1:394991/1‑151 (MQ=255)                                                                                                                            gATT‑CGCCCCCTGGTCGTGATAATCCCCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCa                     <  2:420160/150‑1 (MQ=255)                                                                                                                              tt‑CGCCCCCTGGTCGTGATAACCCCCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGTGACTTTTc                       <  1:457561/146‑1 (MQ=255)                                                                                                                                        tGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACg         >  2:927888/1‑151 (MQ=255)                                                                                                                                             gTGATAATCACCGCCGCCAGCGCTTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGaaa    <  2:124662/151‑1 (MQ=255)                                                                                                                                                aTAACCGCCGCCGCAAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAata  <  1:660166/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                GGCCGTTTTCTGCTTCGACTTCAACAACCAACGTGCCAAGAACGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATT‑CGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGATGTTCTCCATTGTTATGAGGCTTGTAAGTCAAAGGGACTTTTCCATCCCAACAGACGTGAAATA  >  NC_000913/2361189‑2361480

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