Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,397,538	0	C	T	47.8%	13.5 / 11.3	24	K86K (AAG→AAA)	rng	ribonuclease G
Reads supporting (aligned to +/- strand):  ref base C (11/1);  new base T (0/11);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 9.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCCCTGGCATAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTAC  >  NC_000913/3397388‑3397684                                                                                                                                                       |                                                                                                                                                   aGCCCCTGGCATAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGc                                                                                                                                                    <  2:364567/151‑1 (MQ=255)     ccTGGCATAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCtttt                                                                                                                                                >  2:157837/1‑151 (MQ=255)                                          tatCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCatgat                                                                                                           >  2:191265/1‑151 (MQ=255)                                                          cgcccttttAGTGCAAAGCGAATCTTTCACCCCCTCCACAATCAGATCTTGCCCTTGACAAACCATTTCCGAGATGTCGCCCACCTTGAATTGTTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATga                                                                                              <  2:74803/142‑1 (MQ=255)                                                                       ccacgcgtttCTTCCCCCCCCCGCCCCATCAGATTTTCCCCTTGCCAACCCATTTCCGAGTGGTCGCGCCCGGTGAAGAGGTGGGTGTCTTCACCCGCCACACATTCGGTGTGCGGCATGAGGTCGGATGCAGGAAGAAACGCGGcggt                                                                                <  1:55752/140‑4 (MQ=255)                                                                         aaGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATcc                                                                            >  1:281780/1‑151 (MQ=255)                                                                             ggttCTTCCACCCCCTGCCCCATCAGATCTTGCCCCTGCCGAACCATTTCCGAGATGTCGCGCCCCGTGATTTGTTTTTGTTCTCCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGcc                                                                        <  2:179172/148‑1 (MQ=255)                                                                             ggATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGcc                                                                        >  1:348875/1‑151 (MQ=255)                                                                                tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCa                                                                      >  2:14599/1‑150 (MQ=255)                                                                                      ccccccTGCCCCATCAGATCTTCCCCCGGACGACCCAGTTCCGAGATGTCGCGCACCGTGATTTGTTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGccc                                                                       <  1:32453/139‑1 (MQ=255)                                                                                      ccccccTGCCCCACCAGTTTTCGCCCCTGACGATCCATTTCCGAGTTGTCGCGCCCCGTAATTTGTTTTTGTTCTCCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCt                                                                <  2:201410/146‑1 (MQ=255)                                                                                                ccATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCCCCGTGAATTGTTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAgcc                                                       <  2:183277/149‑1 (MQ=255)                                                                                                    cAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCCCCGTAAATTGTTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCt                                                   <  2:240275/149‑1 (MQ=255)                                                                                                         cTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATAcc                                            >  2:279291/1‑151 (MQ=255)                                                                                                             ccccTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGCTGTGCGGCATGATGTCGGCTGCATTAAGAAACGCCGCTTTTTCCAGCCCAATATCTCCAAAAGACGCCTGCCTACCCgg                                         >  1:63876/1‑150 (MQ=255)                                                                                                             ccccGGACAACCCGGTCCCGAGATGTCGCGCCCCGTAATTGGTTTTTTTTCTCCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCAt                                               <  1:239254/144‑1 (MQ=255)                                                                                                                  gACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTac                                   >  2:23689/1‑151 (MQ=255)                                                                                                                        ccAGTCCCGAGATGTCGCGCCCCGTGAATTGTTTTTTTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACg                                <  2:345808/148‑1 (MQ=255)                                                                                                                          aGTCCCGAGATGTCGCGCCCCGTGATTTGTTTTTTTTCTCCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACtt                            <  2:49602/150‑1 (MQ=255)                                                                                                                            ttCCGAGATGTCGCGCACCGTGAATTGTTTTTGTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTa                           <  2:280736/149‑1 (MQ=255)                                                                                                                            tcccGAGATGTCGCGCACCGTGAATTGTTTTTGTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTa                           <  2:358232/147‑1 (MQ=255)                                                                                                                                     gTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGACGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGAACACGACTTACACGACCCtt                 >  1:199586/1‑150 (MQ=255)                                                                                                                                           cACCGTGATTTGCTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAg             <  2:371402/148‑1 (MQ=255)                                                                                                                                                  aaTTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCct    >  1:270135/1‑150 (MQ=255)                                                                                                                                                   aTTGCTTTTGTTCTTCACCCGCCACACATTCGGTGAGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCctac  >  1:210603/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   AGCCCCTGGCATAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGATATTGCCTAC  >  NC_000913/3397388‑3397684

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