Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,206,452	0	C	A	37.0%	24.2 / 10.0	28	intergenic (‑196/‑11)	ttdR/ttdA	transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (3/14);  new base A (10/0);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 3.39e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAAC  >  NC_000913/3206303‑3206596                                                                                                                                                      |                                                                                                                                                  aTTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACt                                                                                                                                                  <  1:161788/151‑1 (MQ=255)              gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCatga                                                                                                                                    <  2:40897/151‑1 (MQ=255)                        aaTATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATTGAG‑GTGTCATGGTGAGCGaaa                                                                                                                           >  1:645497/1‑150 (MQ=255)                                     aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAg                                                                                                             <  1:174125/151‑1 (MQ=255)                                        ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAATGGGG‑GTGGCATGATGAGGGAAAGGAAAAAACaa                                                                                                                >  2:537798/1‑145 (MQ=255)                                            cgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGGG‑GTTGCATGAAGAGCGAAAAGAATAAACAACAAGCaagg                                                                                                       >  2:617127/1‑147 (MQ=255)                                             gTTGATACACGAACTGTGAATAACATATGCCGGGTGGCAAAGAGTATAAAGTCAGAGTTTTAGTGGTAGTTTGTGGCGGTGAGGTTTTCCAAACAAATACCAAAACTGAAG‑TTGCCATGATGCGCGAACTTAATAAGCAACAGGCAgcg                                                                                                       <  1:80175/149‑3 (MQ=255)                                               tGATTCACGAACTGCTAATTACCTATTCACGGTGGCAAACAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAATTCTTCACGCCCAATATAAAAAATTGGG‑TTGCCATACTGAGCGAAACCAACAAGCAAAAgggcgt                                                                                                        >  1:584747/1‑142 (MQ=255)                                                gATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATa                                                                                                   <  2:241759/150‑1 (MQ=255)                                                    cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGtt                                                                                               <  1:505197/150‑1 (MQ=255)                                                       gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACa                                                                                           <  1:265113/151‑1 (MQ=255)                                                         aCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGAGAAATTTTTCAGACAAATTCAAAAACTGGAG‑TTTGCATGATGAGCGTAAGTAAAAAAGAACAAGCAGTGAATa                                                                                                   >  1:127533/1‑141 (MQ=255)                                                                 ttACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATTGAG‑GTGCCATGATGAGCGAAAATAATAAAAAACAAGCAGGGAATAAGTTTACAGAGAtt                                                                                     >  2:327985/1‑147 (MQ=255)                                                                   aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAG‑GTGCCCTGATGAGGGAAAGTAATAAACAACAGGCAGGGAATAAGTTGACAGAGATTGTCg                                                                                 >  2:478628/1‑149 (MQ=255)                                                                     atatTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAAAAAAAACCGGAGATTCCCAT‑ATGAGCGGAAGCAATAAACAAAAAGCAGTGAAAAAAATGACAAAGGTTGgcgc                                                                                >  2:606909/1‑148 (MQ=255)                                                                        ttCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAG‑TTGCCCTGATGAGCGAAAATAATAAACAACAGGCAGTGAAAAAATTGACAGAGATTGTCGCTAACtt                                                                          >  2:440403/1‑151 (MQ=255)                                                                           aCGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAG‑GTGCCATGATTAGCGAAAAAAATAAACAAAAAGCAGGGAAAAAGGTGAAAGAGAATGTCGCTAAAttt                                                                         >  2:568226/1‑149 (MQ=255)                                                                                 gCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCCGACAAAAACAAAAAATGGAG‑TTGCCCTGATGAGCGAAAATAAAAAACAACAAGCAGGGAATAAGTTGACAGAGATTGTTGCCAACTTTACCGCCAt                                                                 >  2:300497/1‑151 (MQ=255)                                                                                  cAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAG‑GTGCCATGATGAGGGAAAGGAATAAGCAAAAGGCAGTGAATAAGTTGAAAGAGATTGTCGCTAACTTTAACGcc                                                                   >  2:473381/1‑148 (MQ=255)                                                                                                  acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGc                                                <  2:553370/151‑1 (MQ=255)                                                                                                   caTTTTTAGTGGTGGTTTGTGGCGGTGAATTTTTCAAGACAAAAACAAAAACTGGAG‑TTCCCATGATGAGCGAAAGTAATAAGCAACGGGCAGTAAATAAGTGGACAGAGATTGTAGTTACCTTTACCGCCATGATTTCTACCCGaa                                                   <  1:454338/147‑1 (MQ=255)                                                                                                            tGGTAGTTTGGGGCGGTGAATTTTTCCAGACAAATACAAAACCTGGAG‑TTGCCATGAGGAGCAAAATAAAAAAGCAACAGGCAGTGAAAAAGTTGAAAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatga                                         <  1:263754/148‑1 (MQ=255)                                                                                                                     gtgGCGGTGAATTTTCCCAGAAAAATAAAAAAACGGGAG‑TTGCCATGATGAGCAAAAGTAAAAAGAAACAGGAATTGATTAGGTTGACAGGGATTGTCGCAAACTTAACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGa                                <  2:627174/148‑1 (MQ=255)                                                                                                                        ggggggAATTTTTCCAAACAATGAAAAAACGTGAAG‑TTCCAGTGATGAGCGAAATAAATAGGCTACAGCCAGTAAATAAGTGAAAGGAGATTGTCGCTAACTTGACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGAt                               <  1:608119/141‑1 (MQ=255)                                                                                                                          ggTGAATTTTTCAAGAAAAATACAAAAGCGGGAG‑TTGCCAGGAGAAGAGAAAGAAATAAGCAACAGGAAGGCAATAAGTTGACAGAAATGGTCGTTACCTTTGCCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTa                         <  2:61433/150‑1 (MQ=255)                                                                                                                            tGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa                       >  2:492304/1‑150 (MQ=255)                                                                                                                                     ccAGACAAATAAAAAAACGGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACGGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAg             <  2:314769/151‑1 (MQ=255)                                                                                                                                                 aaaaaCTGGAG‑TTGCAATAATGAGGGAAAGTAATAAGCAACAGTCAGTAAATAAGTTGACAGAGATTGTCGCGAATTTAACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAAc  <  1:351790/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAG‑TTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAAAGGATGCCGAAAC  >  NC_000913/3206303‑3206596

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