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

CCTGGCATAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGA  >  NC_000913/3397392‑3397674                                                                                                                                                   |                                                                                                                                         ccTGGCATAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCttt                                                                                                                                       <  1:374514/150‑1 (MQ=255)       aTAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGttctt                                                                                                                                >  2:588775/1‑151 (MQ=255)        tAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGttcttc                                                                                                                               >  2:166122/1‑151 (MQ=255)         aaaCACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCa                                                                                                                              <  1:126857/151‑1 (MQ=255)                     gCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGAGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCcacaca                                                                                                                   >  1:520285/1‑150 (MQ=255)                               aGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCg                                                                                                        >  2:187021/1‑151 (MQ=255)                                gCGTGATATCGGTGGTCAGGCGCGCCCCTTTAGTGCAAAGCGGATCTTTCCCCCCCTGCCCCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCCCCGTGAATTGTTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCgg                                                                                                       <  2:148029/151‑1 (MQ=255)                                     atatCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCatg                                                                                                   >  1:668692/1‑150 (MQ=255)                                      tatCGGTGGTCAGGCGCGCCCCTTTAGTGCCAAGCGGATCTTCCCCCCCCGGCCCCATCAGATCTTCCCCCTGACGAACCAGTGCCGAGATGTCGCGCACCGTGAATTGTTTTTGTCCTCCACCCGCCACACATTCGGTGTGCg                                                                                                        <  1:651209/144‑1 (MQ=255)                                        tCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGt                                                                                               >  1:298422/1‑151 (MQ=255)                                          ggtggtCAGGCGCGCCCCTTTAGTGCCAAGCGGATCTTCCCCCACCTCCACCATCAGATCTTCCCCCTGACGAACCATTTCCGAGATGTCGCGCCCCGTAAATTGTTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCg                                                                                             <  2:571231/151‑1 (MQ=255)                                           gtggtCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCgg                                                                                            >  2:214726/1‑151 (MQ=255)                                               tcgggCGCGCACCTTTAGTCCAAGGCGGATCTTCCCCCCCCTCCACCATCAGATCTTCCCCCTGCCGACCCATTTCCGAGATGTCGCGCACCGTGAATTGCTTTGTTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATg                                                                                         <  1:31061/147‑1 (MQ=255)                                                  ggCGCGCCCCTTTAGTGCCAAGGGGATCTTTCCCCCCCTGCCCCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATTTCGCGCCCCGTAAATTGTTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATg                                                                                     <  2:285325/151‑1 (MQ=255)                                                   gcgcgcACTTTTAGTGCCAAGGGGATCTTCCCCCCCCTGCCCCATAAGATCTGGCCCCTGACGACCCAGTCCCGAGATGTCGCGCGCCGTGAATTCTTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATga                                                                                    <  2:170299/151‑1 (MQ=255)                                                   gcgcgcACCTTTAGTGCAAAGCGGATCTTTCCCCCCCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGTTGTCGCGCCCCGTAAATTGTTTTTGTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATga                                                                                    <  2:185662/151‑1 (MQ=255)                                                          ccTTTAGTGCCAAGCGGATCTTTCACCCCCTGCCCCATCAGATCTTCCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGTTTTTGTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAc                                                                              <  2:653302/150‑1 (MQ=255)                                                            tttAGTGCAAGGCGGATCTTCCCCCCCCTGCCCCATCAGATCTTGCCCCTGACGACCCAGTTCCGAGATGTCGCGCACCGTGAATTGTTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAAcg                                                                             <  2:276273/149‑1 (MQ=255)                                                                     aaGCGGATTTTTCCCCCCCTGCACCATCAGATCTGGCCCCTGACGACCCATTTCCGAGATGTCGCGCACCGTGAATTGTTTTTTTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCgg                                                                          <  2:131634/143‑1 (MQ=255)                                                                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCa                                                            >  2:70674/1‑150 (MQ=255)                                                                            tCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCa                                                            >  2:36092/1‑150 (MQ=255)                                                                               ttCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAatat                                                        >  2:186125/1‑151 (MQ=255)                                                                                   cctccgcccccttCAGATATCCCCCCTGACGAACCGTTTCGGAGGTGTCGAGCACAGTGAATGGTTTTTGTTCTTCGCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATc                                                       <  2:211772/136‑1 (MQ=255)                                                                                          cACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAACCCAATATCTACAAAAgcc                                             >  1:182868/1‑151 (MQ=255)                                                                                              accAGTTCTTGCCCCTGACAAACCATTCCCGAGATTTCGGGCACGGTGAATTGCTTTTTTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAgccgc                                           <  1:524685/147‑1 (MQ=255)                                                                                                    ttttcccccTGCCGACCCAGTTCCGAGATGTCGCGCACCGTGAATTGTTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATa                                    <  1:261997/145‑1 (MQ=255)                                                                                                          cccTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGa                              >  2:204438/1‑150 (MQ=255)                                                                                                           ccTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGa                              >  2:203152/1‑149 (MQ=255)                                                                                                                 gAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACg                      >  2:506968/1‑151 (MQ=255)                                                                                                                  aaCCAGTTCCGAGATGTCGCGCCCCGTGAATTGTTTTTGTTCTTCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACg                      <  2:225274/150‑1 (MQ=255)                                                                                                                     cAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGACGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACtt                  >  2:631688/1‑151 (MQ=255)                                                                                                                          ccGAGATGTCGCGCCCCGTGAATTGTTTTTTTTCTCCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACg             <  2:567481/151‑1 (MQ=255)                                                                                                                              gaTGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGAccc         >  2:325805/1‑151 (MQ=255)                                                                                                                                   cgccccccGTGATTTGTTTTTTTTCTCCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGt     <  2:402072/144‑1 (MQ=255)                                                                                                                                      gcACCGTGATTTGTTTTTTTTCTCCCCCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGa  <  1:258993/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                         CCTGGCATAAACACCAGATAGCGAGAAGGGAGCGTGATATCGGTGGTCAGGCGCGCACCTTTAGTGCCAAGCGGATCTTTCACCACCTGCACCATCAGATCTTGCCCCTGACGAACCAGTTCCGAGATGTCGCGCACCGTGAATTGCTTTTGTTCTTCACCCGCCACACATTCGGTGTGCGGCATGATGTCGGATGCATGAAGAAACGCGGCTTTATCCAGCCCAATATCTACAAAAGCCGCCTGCATACCCGGAAGTACACGACTTACACGACCCTTGTAGA  >  NC_000913/3397392‑3397674

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: ‑