Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,079,468	0	A	G	55.9%	20.5 / 40.5	34	D55G (GAT→GGT)	putP	proline:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (4/11);  new base G (19/0);  total (23/11)
Fisher's exact test for biased strand distribution p-value = 4.77e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-06

CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079320‑1079616                                                                                                                                                     |                                                                                                                                                       ccGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGata                                                                                                                                                      <  2:608826/151‑1 (MQ=255)    aTGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATg                                                                                                                                                    <  2:519662/150‑1 (MQ=255)                ttttGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGA‑GCggctgg                                                                                                                                          >  1:294861/1‑147 (MQ=255)                   tgtgTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTGGGGTATGA‑GCGGGTGGGTGtttt                                                                                                                                   >  2:354699/1‑149 (MQ=255)                    gtgtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGG‑GCGGGTGGCTGTTTAt                                                                                                                                  >  1:85100/1‑151 (MQ=255)                      gtCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGGCGGGGATGA‑GCGGGTGGCTGtttttgg                                                                                                                                >  1:567601/1‑146 (MQ=255)                       tCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGGTATTA‑GCGGGTGGCTGTTTATggg                                                                                                                               >  1:520947/1‑151 (MQ=255)                        cTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATgggg                                                                                                                              <  1:244780/151‑1 (MQ=255)                                   gCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGgcg                                                                                                                    >  2:216570/1‑150 (MQ=255)                                       gATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGA‑GGGGGGGGCTGTTTATTGGGGTGGCGGGGGCT‑ggt                                                                                                              >  2:447556/1‑149 (MQ=255)                                        atatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGA‑GGGGCTGGCGGTTTATTGGGGTGGCGGGCGCT‑Ggtt                                                                                                             >  1:432269/1‑151 (MQ=255)                                         tatTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGGGGGTATTA‑GGGGGTGGGTGTTTATGGGGGTGGCGGGCGCG‑Ggttt                                                                                                            >  2:403385/1‑151 (MQ=255)                                            tgattgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGGATGT‑GGGGGTGGCGGTTATTGGGGGTGCCGGGGGCT‑Gttttttt                                                                                                         >  1:421078/1‑149 (MQ=255)                                              attgGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTtttctt                                                                                                        <  2:294856/150‑1 (MQ=255)                                                 gggTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGTATGA‑GCGGGTGGGTGTTTATGGGGGTGCCGGGGGCT‑GTtttttttttc                                                                                                     >  1:222799/1‑148 (MQ=255)                                                                                    cgacTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACa                                                                 <  2:634546/151‑1 (MQ=255)                                                                                      acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGGGGCGGGGATTG‑GGGGGTGGCTGTTTATAGGGGGGGCGGGCGGT‑GGTTTTTTTTTCGGGGTTTCCGCCAGCTGGATGGCCactcggctgtcaat                                                               >  2:373868/1‑138 (MQ=255)                                                                                      acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACAt                                                                <  2:45179/150‑1 (MQ=255)                                                                                      acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGTCGGGTATGA‑GGGGGTGGGTGTTTATGGGGTTGGCGGGGGGC‑GTTTTTTTTTTCTGGGTTTTCGGAAGCTGGGTTGCCCCTTGCCTGGCAtt                                                               >  2:638926/1‑151 (MQ=255)                                                                                      acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATTA‑GCGGCTGGGTGTTTATGGGGGTGGCGGGGGCG‑GTTTTTTTTTTCGGGGTTTCCGGAAGCTGGGTGGCCATTTGCCTGTCAtt                                                               >  2:562459/1‑151 (MQ=255)                                                                                      acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGA‑GGGGCTGGCTGTTTATGGGGGTGCCGGGGGCT‑GTTTTTTTTTTCGGGGTTTTCGAAAGCTGGATCGCCCTTGGCCTGACAtt                                                               >  1:478964/1‑151 (MQ=255)                                                                                      acTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGA‑GCGGCCGGCTGGTAATGGGGGTGTCGGGGGGTGGTTTTTTTTTTCGGG‑TTTTCGAAAGCGGGATCTGCATTGGCCTGCCAtt                                                               >  1:10565/1‑151 (MQ=255)                                                                                         atatTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAgg                                                            <  2:296607/151‑1 (MQ=255)                                                                                           atTCTGTGTGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGGGGCGGGTATGAGGGGGCGGGTGGTTATTGGGTTTGCGGGGGCGG‑GTTTTTTTTTCGGGGGTTTCCG‑AAGCTTGGTCGCCCTTGGGCTGACATTTGgcg                                                          >  2:500274/1‑151 (MQ=255)                                                                                            ttCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGGGCGTCGGGTATGA‑GGGGCTGGCTGTTAATGGGGGTGGCGGGCGCG‑GTTTTTTTTTTCGGGGTTTTCGAAAGCTGGATCGCCATTGGCCTGGCATTAGgcgg                                                         >  1:517515/1‑150 (MQ=255)                                                                                                         tAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTg                                             <  2:423945/150‑1 (MQ=255)                                                                                                             cTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGc                                        <  2:95260/151‑1 (MQ=255)                                                                                                                     aTTCGTGACGGCATTATCGGCGGGTGCGGCGGGTATGA‑GCGGCGGGGTGTTTATGGGGGTGGCGGGCGCC‑GTTTTTTTTTCCGGGGTTTTCGAAAGGTGGGTCGCCATTGGCCTGACATTTGGCGCGGGGGTTTACTTGGAGCtggtggg                                 >  1:262895/1‑149 (MQ=255)                                                                                                                          tGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGGTTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCg                           <  2:535887/151‑1 (MQ=255)                                                                                                                             cGGCATTATCGGCGGGTGCGGCGGGTATGA‑GCGGCTGGCTGTTTATGGGGGTGCCGGGGGCT‑GTTTTTTTTTTCGGGGTTTTCGAAAGCTGGGTCGCCCTTTGCCTGGCATTAGGCGCGTGGGTTTAATTGGAGCTGGTGGGCgggggggt                        >  1:425220/1‑149 (MQ=255)                                                                                                                               gCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGGTGCCGGGCGCT‑GTTTTTCTTTTCGGGGTTTCCGAAAGCTGGATCGCCCTTGGCCTGGCATTAGGCGCGGGGGTTAACTGGAAGCTGGTGGCCGGGCGGGTg                       >  1:571987/1‑150 (MQ=255)                                                                                                                                               cGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATa      >  2:376544/1‑151 (MQ=255)                                                                                                                                               cGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATa      >  1:2392/1‑151 (MQ=255)                                                                                                                                                 tCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATac     <  1:344399/150‑1 (MQ=255)                                                                                                                                                    gATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATacaac  <  1:309209/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CCGATGTTGGTGACATTTTGTGTCTATATCTTTGGCATGATATTGATTGGGTTTATCGCCTGGCGATCAACGAAAAACTTTGACGACTATATTCTGGGCGGTCGTAGTCTTGGGCCATTCGTGACGGCATTATCGGCGGGTGCGTCGGATATGA‑GCGGCTGGCTGTTAATGGGGTTGCCGGGCGCT‑GTTTTTCTTTCCGGGATTTCCGAAAGCTGGATCGCCATTGGCCTGACATTAGGCGCGTGGATTAACTGGAAGCTGGTGGCCGGGCGGTTGCGTGTGCATACCGAATACAAC  >  NC_000913/1079320‑1079616

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