Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,624,517	0	G	T	38.7%	40.0 / 11.3	31	E492* (GAA→TAA)	ppk	polyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base G (5/14);  new base T (12/0);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 8.83e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCT‑ATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTG  >  NC_000913/2624373‑2624648                                                                                                                                                  |                                                                                                                                      aaaGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAAc                                                                                                                                 <  1:333716/150‑1 (MQ=255)  aaGAAAACGGTGAAGTTGAACTTTATTACCACATCGGGACCGGAATCTTACACGAAAAAACCGGGCTGCTGTAAACTCACTAAGCGTTGTTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GTGTATTTAACTTTATTG‑AAAcccc                                                                                                                               <  1:145251/151‑1 (MQ=255)  aaGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGATTTTTTAATTTTTTG‑AAAAacc                                                                                                                               >  1:295684/1‑151 (MQ=255)   aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTAAATTTTTTT‑AAAAacc                                                                                                                               >  2:226423/1‑150 (MQ=255)      aaaCGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTACTTTTTTT‑AAAAccc                                                                                                                               >  1:687115/1‑147 (MQ=255)       aaCGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATAcc                                                                                                                          <  1:512349/151‑1 (MQ=255)         cGGTGTAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTACTTTTTTT‑AAAAacc                                                                                                                               >  2:111437/1‑144 (MQ=255)         cGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTTACTTTATTT‑AAAATCCaccac                                                                                                                          >  1:140873/1‑145 (MQ=255)              aaGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGt                                                                                                                   <  2:391672/151‑1 (MQ=255)                  ggtgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTTACTTTATTG‑AAAACCCATAACGTcc                                                                                                                      >  2:503507/1‑144 (MQ=255)                  ggtgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTAACTTTTTTT‑AAAAAACCTACCCTCCgggg                                                                                                                  >  2:258412/1‑146 (MQ=255)                                    cGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTACTTTTTTT‑TAAAACCCCACCCCCCGGTGGCACTTTATTATTTTATGGt                                                                                              >  1:528813/1‑150 (MQ=255)                                      ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATc                                                                                           <  1:269435/151‑1 (MQ=255)                                      ggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGTTTTAACATTATTT‑TAAACCCATAACCTCCGGGCACATTTGATTTTTTAATGGTAt                                                                                            >  1:269628/1‑150 (MQ=255)                                           gggAGCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATcgccgc                                                                                      <  1:623545/151‑1 (MQ=255)                                              aaCTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATCCGTTGCTAACCTCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCaaa                                                                                   <  1:185419/151‑1 (MQ=255)                                                   tAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTc                                                                               <  2:604186/150‑1 (MQ=255)                                                     aCGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACt                                                                                <  1:553880/147‑1 (MQ=255)                                                       gAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgcc                                                                          <  2:274196/151‑1 (MQ=255)                                                            aaCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTATTT‑ACTGTTTTT‑TAAACGCATACCCTCCGGTGGAATTTGTTTTTTTTATGGTATTACCGGACAACCCcctcccc                                                                        >  1:275714/1‑146 (MQ=255)                                                                       tttATACTGACTATTCGTTGCGGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCT‑ATTGTAt                                                              <  1:267739/146‑1 (MQ=255)                                                                         tataCTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTACTTTTTTT‑AAAAACCCTACCGGCCGGGGGCATTTTTTTTTTTTATTGTATTGCCCCCAAAACCCCCCCCCCC‑ATTGTTTGaaa                                                          >  2:272263/1‑148 (MQ=255)                                                                                   aTTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTACTTTTTTT‑TAAACCCATAACCTCCGGGGACATTTGTTTTTTTTATGGTATCGCCGCCAAACTCCCCCCGCCC‑ATTGTTTGAAAAGGTGGGcccc                                               >  2:233506/1‑147 (MQ=255)                                                                                     tCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCT‑ATTGTATGAAATGGTGGACCGCgag                                            <  2:660568/150‑1 (MQ=255)                                                                                        ttGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTACCTTTTTGAAAAACCCAAACCTCCCGGGGAAATTTG‑TTATTTTATTGGATTCACGCAAAAACCCCCCCCCCTCTTTGTTTGAAAGGGGGGcccc                                                >  1:244399/1‑143 (MQ=255)                                                                                            tGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGGATTTTACTTTATTT‑AAAACCCCTACCGTCCCGGGGCATTTTGTTTTTTAATTGTATAGCCCCCAAAACCCCCCCCCCT‑ATTTTATGAAATGGGGGGCCGCGGGAACGcc                                      >  2:465894/1‑149 (MQ=255)                                                                                                        cgcgcATCACCAACGAAGTACGGC‑GGGGATTTTACTTTTTTT‑AAAAACCAAACCGGCCGGGGACATTTGATTATTTAATGGTATCGCCGCAAAAATCCCCCCCCCC‑ATTTTATTAAATGGGGGAACGCGAGATCGCCCACGCCCAGCaa                          >  2:125490/1‑149 (MQ=255)                                                                                                                           aCGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCT‑ATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCAc      <  2:140921/150‑1 (MQ=255)                                                                                                                              gc‑ggGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCT‑ATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTg  <  1:624809/151‑1 (MQ=255)                                                                                                                               c‑ggGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCT‑ATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTg  >  1:622429/1‑150 (MQ=255)                                                                                                                               c‑ggGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGAATATTTAATGGTATCGCCGAAAAACTCCCGCCGCCT‑ATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTg  >  1:532533/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                      AAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGC‑GGGTATTTAACTTTATTG‑AAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCT‑ATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTG  >  NC_000913/2624373‑2624648

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