Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	44,763	0	T	G	43.5%	72.6 / 11.2	46	V195G (GTT→GGT)	fixC	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (7/19);  new base G (20/0);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 2.49e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGATGGCGATGTGATTGAAGCGAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCTGGTTTGTGGTTTGCATCATCTGC  >  NC_000913/44614‑44912                                                                                                                                                      |                                                                                                                                                        aGCCGATGGCGATGTGATTGAAGCGAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGtt                                                                                                                                                        <  2:329430/151‑1 (MQ=255)                    gAAGCGAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGGGAAGGAACAGAGCGAGTTACTGCAGTAGGGCA‑TTGAAGACCGGTTTCa                                                                                                                                      >  2:422911/1‑149 (MQ=255)                     aaGCGAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGtt                                                                                                                                   <  1:215773/151‑1 (MQ=255)                      aGCGAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTg                                                                                                                                  <  1:671035/151‑1 (MQ=255)                        cGAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGGGAAAGAAATGAACGAGATTACGCAGTCGGTTA‑ttttaaaccggtt                                                                                                                                         >  1:374377/1‑131 (MQ=255)                         gAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAg                                                                                                                               <  1:25427/151‑1 (MQ=255)                         gAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGATTGGGTGAAGGGGATGATCGAGTTTACGAAGACGGGGA‑TTTGAGAACGGTTTTAGTTggc                                                                                                                                >  2:720629/1‑148 (MQ=255)                             aCGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTa                                                                                                                           <  2:824734/151‑1 (MQ=255)                                    tCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGGAGGAACTGATCGAGTTACCGAAAGCGGGTT‑TTTAAGACCGTTTTTAGTTGCAGGGGAAACA‑ggg                                                                                                                   >  1:827477/1‑151 (MQ=255)                                      cTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTTGGGGGGAGGGAATGATCGAGGTTCCGCAGGGGGGGA‑TTTTAGACCGGTTTCTGTTTTAAGGGAATCA‑ggggg                                                                                                                 >  2:500816/1‑151 (MQ=255)                                      cTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGGAGGAACTGATCGAGGTAACGCAGTCGGGTA‑TTTAAGAACGTTTTCAGTTTCAGGGGAATAA‑ggggg                                                                                                                 >  2:410058/1‑151 (MQ=255)                                        tGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGggc                                                                                                                <  2:919756/150‑1 (MQ=255)                                           tGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGAAGGAAAAGAACGAGGTTACGGAGTCGGGTT‑TTGAAAAACGGTTTTAGTTGCAGGGGAAATA‑gggggggg                                                                                                              >  1:343811/1‑149 (MQ=255)                                            gATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGGAGGAACTGAACGAGGTACCGAAGTCGGGTT‑TTTAAAAACGTTTTTAGTTTCAGGGGAATCA‑gggggggg                                                                                                              >  1:671218/1‑148 (MQ=255)                                                     aaCTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGTGGAGGGAATGGTCGGGGTACCGAAGGCGGGTA‑TTTAAGAACGGTTTTTGTTTCAGGGGAATCA‑GGGGGGGGGTTGCCTGttt                                                                                                   >  1:33332/1‑150 (MQ=255)                                                     aaCTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTg                                                                                                  <  1:416748/151‑1 (MQ=255)                                                        tCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTTGGGGGGAGGAACAGAACGAGGTACCGCAGTCGGGTT‑TTTAAAACCGCTTTTTGTTTCAGGGGAAACA‑GGGGGGGGCTTGCCTGttt                                                                                                   >  2:316518/1‑147 (MQ=255)                                                         ccATCCTTGCCGAAAAATTGTGGATGGCAAAACGCGTCAAACCGACGGATGTGGCTGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCg                                                                                                <  2:295411/149‑1 (MQ=255)                                                         ccATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCggg                                                                                              <  2:293027/151‑1 (MQ=255)                                                             ccTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGTGGAGGGAATGGGCCAGTTTCCGGAGTGGGGTT‑TTTGAGAACGGTTTTAGTTGTAGGGGAATTA‑AGGGGGGGGTTGTCTtttttgcgg                                                                                              >  1:159025/1‑142 (MQ=255)                                                                 gCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATcaccca                                                                                      <  1:737707/151‑1 (MQ=255)                                                                  ccGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGGAGGGAATGATCGAGGTACCGAAGGCGGGTT‑TTGAAGACCGGTTTTTGTTTCAGGGGAATCA‑GGGGGGGGCTTGCCTGTTTTCGGGGTcacccac                                                                                     >  1:314486/1‑151 (MQ=255)                                                                  ccGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGATGGGGTGAAGGGAATGATCGAGTTTCCGAAGTCGGGTA‑TTTAAGAACGGTTTTTGTTGCAGGGGAATCA‑GGGGGGGGCTTTTCTGTTTTGGGGGTcaccc                                                                                       >  1:810080/1‑149 (MQ=255)                                                                   cGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGGAGGAAAAGATCGAGGTTCCGGAGGCGGGTT‑TTTGAGAACGGTTTTTGTTTCAGGGGAATCAGGGGGGGGGCTTG‑CTGttttggggggcaaccccc                                                                                    >  2:60206/1‑136 (MQ=255)                                                                   cGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGAAGGAAATGATCGAGTTTCCGAAGGCGGGTA‑TTTGAGAACGTTTTTAGTTGCAGGGGAATCA‑GGGGGGGGGTTGCCTGGTTTCGGGGTcacccccc                                                                                    >  1:342070/1‑151 (MQ=255)                                                                      aaaaTTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGAt                                                                                 <  2:156735/151‑1 (MQ=255)                                                                      aaaaTTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGGGGAGGAACTGATCGGGTTACCGAAGGCGGGTA‑TTGAAAACCGGTTTTAGTTGCAGGGGAAACA‑GGGGGGGGCTTGCCTGTTTGCGGGGTCACCCACCCAt                                                                                 >  1:583459/1‑151 (MQ=255)                                                                       aaaTTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGGAGGAAAAGATCCAGGTTACGAAGGCGGGTA‑TTGTAGAACGGTTTTAGTTTGAGGGGAAACA‑TGGGGGGGGTTGTCTGTTTTCggggtaaaacac                                                                                     >  2:771909/1‑137 (MQ=255)                                                                           tGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGGGGGGAGGAAAAGATTGAGTTACCCGAGTCGGGTATTTTAAAACGGTTTTC‑GTTGGGGGGGGAACA‑AGGGGGGGGTTTCCTGTTTTCGGGGGcacccccc                                                                                    >  2:655081/1‑143 (MQ=255)                                                                                   gCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGGAGGGAATGATCGAGTTACCCAAGACGGGTA‑TTGAAGAACGGTTTCAGTTGCCGGGGAATCA‑GGGGGGGGGGTGCCTGTTTTCGGGGGCACCCACCCATGGCCTGATggggg                                                                    >  1:424834/1‑149 (MQ=255)                                                                                         cgcgTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCt                                                              <  2:827069/151‑1 (MQ=255)                                                                                             tCAAACCGACGGATGTGGCGGTTGGGGGGGAGGAACTGATCGAGGTTCCGAAGTCGGGTA‑TTTGAGACCGTTTTTAGTTGCAGGGGAAACA‑GGGGGGGGGTTGCCTGTTTTGGGGGTCACCCACCCATGGGCTGGTGGGCGGCGGCTTcc                                                           >  1:748112/1‑150 (MQ=255)                                                                                              cAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCtt                                                         <  1:786963/151‑1 (MQ=255)                                                                                                      cGGATGTGGCGGTTGGCGTGGAGGAACTGGTCGAGGTTCCCCAGTCGGGTT‑TTTAAAAACGGTTTTTGTTTCCGGGGGATCA‑GGGGGGGGCTTTCCTGTTTTCGGGATCACCCACCCATGGCCTGGTGGGGGGGGGCCTCCttttt                                                      >  2:191485/1‑144 (MQ=255)                                                                                                             ggCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAac                                          <  1:321661/151‑1 (MQ=255)                                                                                                                ggTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACAcc                                       <  1:245207/151‑1 (MQ=255)                                                                                                                ggTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACAcc                                       >  2:520778/1‑151 (MQ=255)                                                                                                                    ggCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCt                                     <  2:11990/149‑1 (MQ=255)                                                                                                                      cGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCg                                 >  1:620826/1‑151 (MQ=255)                                                                                                                         gAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCt                               <  2:118726/150‑1 (MQ=255)                                                                                                                                  gATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCtggt                      <  2:695289/150‑1 (MQ=255)                                                                                                                                   aTCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCtggttt                    >  1:236067/1‑151 (MQ=255)                                                                                                                                      gAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCtggtttgtg                 >  1:629071/1‑151 (MQ=255)                                                                                                                                          tACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCtggtttgtggtt              <  2:828378/150‑1 (MQ=255)                                                                                                                                             cGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCTGGTTTGTGGTTTGc           >  1:706622/1‑150 (MQ=255)                                                                                                                                              gAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCTGGTTTGTGGTTTGcat         <  1:40101/151‑1 (MQ=255)                                                                                                                                                   cGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCTGGTTTGTGGTTTGCATCATCt    >  2:154256/1‑151 (MQ=255)                                                                                                                                                      ttA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCTGGTTTGTGGTTTGCATCATCTGc  >  2:364694/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                        AGCCGATGGCGATGTGATTGAAGCGAAAACGGTGATCCTTGCTGATGGGGTGAACTCCATCCTTGCCGAAAAATTGGGGATGGCAAAACGCGTCAAACCGACGGATGTGGCGGTTGGCGTGAAGGAACTGATCGAGTTACCGAAGTCGGTTA‑TTGAAGACCGTTTTCAGTTGCAGGGTAATCA‑GGGGGCGGCTTGCCTGTTTGCGGGATCACCCACCGATGGCCTGATGGGCGGCGGCTTCCTTTATACCAATGAAAACACCCTGTCGCTGGGGCTGGTTTGTGGTTTGCATCATCTGC  >  NC_000913/44614‑44912

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