Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,355,731	0	A	G	75.9%	45.1 / 12.4	30	T227T (ACT→ACC)	dtpC	dipeptide/tripeptide:H(+) symporter DtpC
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base G (13/9);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGTACTTCAATATTGAATGCCTGACGATTCACAAAGCGATCGATAAACAAGCTGATGGTACTGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGATGTTCGGGGAATTTGATCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGAGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACTGGTAAGGCAAATTTGACACTAGTGAGCGCTTTTTTATCCATACTACGTGTGGACTGGAAATGACGATGACCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACATGCCATCCATACCACTGAGCAGCCAGGCCGCAGGCGATTGGGGCTGCGATAGAACCGATATTGCCC  >  NC_000913/4355494‑4355977                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       aGGTACTTTAATATTGAATGCCTGACGATTCACAAAGCGATCGATAAACAAGCTGATGGTACTGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGATGTTCGGGGAATTTGATCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGAGTAAAAAATACt                                                                                                                                                                                                                                            <  2:43950/250‑1 (MQ=255) aGGTACTTCAATATTGAATGCCTGACGATTCACAAAGCGATCGATAAACAAGCTGATGGTACTGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACt                                                                                                                                                                                                                                            <  1:99798/250‑1 (MQ=255)     aCTTCAATATTGAATGCCTGACGACTCACAAAGCGATCGCTACAAAAGCTGATGGTAATACCGCCCTGTTGAGCCAGTACCCCGAACAATGTACAGACAAACATAAACAATACACTTTGCCAAAGAACCCGGCACTGTTCGGGGAATTTGATCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCAGGACCGGTCGTTCTCAAGAAGCAAAGTAACAAATACTGGGGc                                                                                                                                                                                                                                       >  2:361088/1‑251 (MQ=255)                               cacaAAGCGATCGATAAACAAGCTGATGGTACTGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGATGTTCGGGGAATTTGATCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGAGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGc                                                                                                                                                                                                             <  1:174017/251‑1 (MQ=255)                                                 cAAGCTGATGTGACTGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATGTGCCAAAGAGCACGGCGATGGTCGGGTAATTTGATCAGCATGCGGGCAATGATGTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATAGCCGGACCAGTCGTTCTCCAGCAGCAGAGTAAAAAATACTGGGGCTAAACAGAGCAGCACCACTAACCATCTCCATACTGGTAAGGCaa                                                                                                                                                                                           <  2:154762/251‑1 (MQ=255)                                                    gCTGATGGTACTGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAAtt                                                                                                                                                                                        <  1:386921/251‑1 (MQ=255)                                                               tGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAAtttcacgcaagtg                                                                                                                                                                             >  1:268525/1‑241 (MQ=255)                                                                     ccTGTTGTGCCAGTACCCAGAACAATGTAACGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGAGCTCGGGGCATTTGCACATCATGCGGGCAATGATTTGTGCGGCAATGAGGCCAACGATCGCCAGCAAATATCCGGACCCGTCGTTCTCCAGAAGCAGGGGAAAAAATACTGGGGCTAAACCGAGCATCCCCACTAACCAGCTCCATACCGGTGAGGCACATTTAAAGCTAGTGAGCGCt                                                                                                                                                                       >  2:351027/1‑251 (MQ=255)                                                                      cTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCtt                                                                                                                                                                      >  1:160574/1‑251 (MQ=255)                                                                               cAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTGTTTTATCCa                                                                                                                                                             <  2:202205/251‑1 (MQ=255)                                                                                   aCCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCAtact                                                                                                                                                         <  1:368557/251‑1 (MQ=255)                                                                                     ccAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGATGTTCGGGGAATTTGATCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGAGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACTGGTAAGGCAAATTTGACACTAGTGAGCGCTTTTTTATCCAtactac                                                                                                                                                       <  1:262876/251‑1 (MQ=255)                                                                                          aCAATATCCCGACAAACATCAACAATACAATTTGCCAACGAGCACGGCGATGTTCGGGGAATTTGATCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGAGTAACAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACTGGTACGGCAAATTTGACACTAGTGAGCGCTTTTTTATCCATACTACGtgtg                                                                                                                                                  >  1:301011/1‑251 (MQ=255)                                                                                            aaTGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCCAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTgg                                                                                                                                                 >  1:142007/1‑250 (MQ=255)                                                                                            aaTGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGGTGTTCGGGGAATGTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAATGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGa                                                                                                                                                <  1:370788/251‑1 (MQ=255)                                                                                                                caatacaatTTGCCAAAGAGCACGGCGGTGTTCGGGGAATTTGACCATCATGAGGGCAATGATTTGTGCGGCAATGAGGCAAAAGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGc                                                                                                                            >  1:365146/1‑251 (MQ=255)                                                                                                                             cAAAGAGCACGGCGATGTTCGGGGAAGTTGATCATCAGGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGAGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACTGGTAAGGCAAATTTGACACTAGTGAGCGCTTTTTTATCCATACTACGTGTGGACTGGAAATGACGATGACCGCTTAAGAAAATcaa                                                                                                               <  2:85918/251‑1 (MQ=255)                                                                                                                                           gtgtgcGGGGAACTTGAACATCATGCGTGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCTGTCTCCAGCAGCAGTGTAAAAAAGACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAAATTCACGCTAGTGAGCGCTTTTTTAACCATACTACGCCTGGACTGGAAATGACGATGCCCGCTTAATAAAATCAGCAAACCGATAAACa                                                                                                 <  2:365146/246‑1 (MQ=255)                                                                                                                                                      aTTTGACCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCAccc                                                                                      >  1:56505/1‑251 (MQ=255)                                                                                                                                                                      ggCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCaa                                                                      <  1:236953/251‑1 (MQ=255)                                                                                                                                                                       gCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCTCAAGGGCAAAGCCAAc                                                                     <  2:160574/251‑1 (MQ=255)                                                                                                                                                                        cAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAAc                                                                     >  1:338074/1‑250 (MQ=255)                                                                                                                                                                        cAATGATTTGGGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACa                                                                    <  1:351027/251‑1 (MQ=255)                                                                                                                                                                                           gAGGCAAACGATCGCCAGCAAATCTCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACCTGCCATCCGTACCACTg                                                   >  2:109249/1‑249 (MQ=255)                                                                                                                                                                                           gAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACATGCCATCCATACCACtacg                                                 >  2:243033/1‑248 (MQ=255)                                                                                                                                                                                                         ccAGCAAATATCCGGACCAGTCGATATCCTGCAGCAGGGTACACACTACTGGGGCTCAACAGAGCATCACCACTAACCAGCTCCATACAGGTAAGGCAAATTTAACGCTAGTGCGCGCTTTTTTATCCACACTACGCGTGGCCTGGTAATGCCGATGCCCGCTTCAGAAAATAAACAAACCGATAAACACGCCGCACCCCGCAAGGGAAAAGCCAACAAGCCATCCATACCACTGAGCAGCCAGGCCGc                                     >  2:21568/1‑249 (MQ=255)                                                                                                                                                                                                                      ggACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACATGCCATCCATACCACTGAGCAGCCAGGcc                                       <  1:37032/234‑1 (MQ=255)                                                                                                                                                                                                                      ggACCAGTCGTTCTCCAGCAGCAGGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACATGCCATCCATACCACTGAGCAGCCAGGcc                                       >  2:37032/1‑234 (MQ=255)                                                                                                                                                                                                                                         tgcagGGTAAAAAATACTGGGGCTAAACAGAGCATCACCAATAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAAAATGCCGCCACCCGCAAGGGCAAAGCAAACCTGACATCCATACCACTGAGCAGCCAGGCCGCAGGCGCTTGGGGCGGCGATAGCACCGatatatcc   >  2:77434/2‑247 (MQ=255)                                                                                                                                                                                                                                          gcagGGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACCGGTAAGGCAAATTTCACGCTAGTGAGCGCTTTTTTATCCATACTACGCGTGGACTGGAAATGACGATGCCCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACATGCCATCCATACCACTGAGCAGCCAGGCCGCAGGCGATTGGGGCGGCGATAGAACCGATATTGccc  >  1:273331/1‑251 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       AGGTACTTCAATATTGAATGCCTGACGATTCACAAAGCGATCGATAAACAAGCTGATGGTACTGCCGCCCTGTTGTGCCAGTACCCAGAACAATGTCCCGACAAACATCAACAATACAATTTGCCAAAGAGCACGGCGATGTTCGGGGAATTTGATCATCATGCGGGCAATGATTTGTGCGGCAATGAGGCAAACGATCGCCAGCAAATATCCGGACCAGTCGTTCTCCAGCAGCAGAGTAAAAAATACTGGGGCTAAACAGAGCATCACCACTAACCAGCTCCATACTGGTAAGGCAAATTTGACACTAGTGAGCGCTTTTTTATCCATACTACGTGTGGACTGGAAATGACGATGACCGCTTAAGAAAATCAACAAACCGATAAACATGCCGCCACCCGCAAGGGCAAAGCCAACATGCCATCCATACCACTGAGCAGCCAGGCCGCAGGCGATTGGGGCTGCGATAGAACCGATATTGCCC  >  NC_000913/4355494‑4355977

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