Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,151,589	0	A	T	40.0%	69.7 / 29.6	50	L135F (TTA→TTT)	yegJ	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (4/26);  new base T (19/1);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 5.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCC  >  NC_000913/2151446‑2151735                                                                                                                                                 |                                                                                                                                                     tttATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTAACATTGATTCtttttt                                                                                                                                                    >  2:366568/1‑146 (MQ=255)   tATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑c                                                                                                                                             <  1:63567/150‑1 (MQ=255)     ttGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTAACATTTATTCtttttt                                                                                                                                                    >  1:998248/1‑142 (MQ=255)       ggCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTAt                                                                                                                                         <  1:189589/150‑1 (MQ=255)        gCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATAACCTTTATTCTTTTTTTGTCC‑CCCATcc                                                                                                                                       >  1:528716/1‑151 (MQ=255)         cGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTACCCTTTATCCttttttttttacccctcc                                                                                                                                        >  1:540706/1‑138 (MQ=255)         cGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTACCATTTATTCTTTTTTTGtgccc                                                                                                                                             >  2:386395/1‑141 (MQ=255)         cGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTAACATTTTTTCTTTAttttttacccat                                                                                                                                          >  2:541248/1‑138 (MQ=255)               cAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAaagaa                                                                                                                                <  1:683750/151‑1 (MQ=255)                  cAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTAACATTTTTCCTTTTTTTGTTA‑CCCATCTAaaaaaaa                                                                                                                               >  2:938460/1‑148 (MQ=255)                     tGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGGTATAACATTTATTCTTTTTTTGGTA‑CCCATccccaaaaaaaa                                                                                                                             >  2:442231/1‑138 (MQ=255)                        cccTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCcgac                                                                                                                       <  1:456215/151‑1 (MQ=255)                           tGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATTTCCATGTTCCTTTTTTTGCTG‑CCTCTCCAAAAAAACAAGCcgac                                                                                                                       >  1:233067/1‑148 (MQ=255)                            ggTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTACCATTTATTCTTTTTTTGTTA‑CCCATTCCAAAAAAAAAACcgacga                                                                                                                     >  1:537796/1‑149 (MQ=255)                               aTTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTa                                                                                                                 <  1:9712/150‑1 (MQ=255)                                 tGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTAACATTTATCCTTTTTTTGTTACCCTTTCCAAAA‑AAAAAACCGACGAGCCaa                                                                                                                >  1:705049/1‑149 (MQ=255)                                     aCAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCG‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCaa                                                                                                          <  2:142005/151‑1 (MQ=255)                                             aGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGtt                                                                                                  <  2:392844/151‑1 (MQ=255)                                                  atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAGATTAGTTCGCTa                                                                                             <  2:8656/151‑1 (MQ=255)                                                  atGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTa                                                                                             <  2:646568/151‑1 (MQ=255)                                                        aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTAACAATGATTCTTTTTTTGTTA‑CCTATCCCAAAAAAGAAAGCGACGAGGTAAAAACAAAAATAGTTTTCTaa                                                                                            >  2:15167/1‑146 (MQ=255)                                                        aGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGc                                                                                       <  1:283493/151‑1 (MQ=255)                                                              ggAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTtat                                                                                 <  1:45164/151‑1 (MQ=255)                                                                  aaTCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTAtattt                                                                             <  2:583894/151‑1 (MQ=255)                                                                         aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTACCATTTATCCTTTTTTTGTCA‑CCTATCCAAAAAAAGAAAGCGACGAGCTAAAAGGAAAAATAGTTTGCTAATCAGCCGCTTTAtttttttttt                                                                      >  1:877570/1‑151 (MQ=255)                                                                         aaagaaGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGtt                                                                       <  1:225539/150‑1 (MQ=255)                                                                          aagaagATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATAACATTGGTCCTTTTTTTGTTA‑ACTTTCCCAAAAAAGAAACCGACGAGCTAAAAGGAAAAATTGTTCGCTAAACAGCCCTTATAtttttttttt                                                                      >  1:994114/1‑150 (MQ=255)                                                                             aagaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTAACATTTATTCTTTTTTTGTCA‑CCCATCCAAAAAAAAAAACCGACGAGCTAAGAGCAAAAATAGTTTGCTAATCAACCGTTTTATTTTTTTTTTtaa                                                                   >  1:60375/1‑150 (MQ=255)                                                                               gaTATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTg                                                                <  2:569118/151‑1 (MQ=255)                                                                                 tatTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTACCCTTTATTCTTTTTTTGTTA‑CCTATCCAAAAAAAAAAGACGACGAGCTAAGAGCAAAAATAGTTTGCTAATCAACCGGTATATTTTTTGTTTtaa                                                                   >  1:865678/1‑146 (MQ=255)                                                                                   tttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTAACCAATTATGCTTTTTTTGTGA‑CCCATCCAAAAAAAAAAGGCGAACAGCCAAAAACAAAAATAGTTCGCCAAACAACCGGTATATTTTTTTTTTtaaatagggg                                                            >  1:111311/1‑144 (MQ=255)                                                                                       aGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGGTTATACCTTTCTTTTTTTTTTGG‑GA‑CCCATCCAAAAAAAAAAACCGACCGGGTAAAAGCAAAAATATTTTGCTTATCAACCGGTATTTTTTTTTTTTTAAATTGC‑ggg                                                          >  2:939410/1‑147 (MQ=255)                                                                                       aGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTACCATTTATCCTTTTTTTGTCA‑CCCATCCAAAAAAAAAAACCGCCGAGCTAAAAGGAAAAATAGTTTGCTAATCAGCCGGTATATTTTTTGTTTtaaaatgggg                                                            >  1:559543/1‑140 (MQ=255)                                                                                       aGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGa                                                       <  1:291969/151‑1 (MQ=255)                                                                                        gACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTACCATTTATCCTTTTTTTTTGA‑CCCATCCAAAAAAAGAAAGCGACGAAGCAAAAGAAAAAATAGTTTGCTAATCAGCCGTTATATTTATTTTTTTAAATTGC‑gg                                                           >  2:760911/1‑146 (MQ=255)                                                                                        gACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATAACATTTATCCTTTTTTTGGCA‑CCTATCCAAAAAAAAAAGGCGCAGAGGCAAGAGGAAAAATAATTTGCTAAACAGCCGTTATTTTTTTTTTTTAAAATTGC‑gg                                                           >  1:632831/1‑146 (MQ=255)                                                                                         aCTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTTTATCATTTTTTCTTTTTTTGGTA‑CCTATCCCAAAAAAAAAGGCGACGAGGTAAAAGCAAAAATAATTTGCTAATCAACCGGTATTTTTTTTTTTTTAAATGGC‑ggg                                                          >  2:528721/1‑145 (MQ=255)                                                                                            ggATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAAttgt                                                  <  2:545111/151‑1 (MQ=255)                                                                                                gTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATtgttgt                                               <  1:934580/150‑1 (MQ=255)                                                                                                  ataCGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCt                                             <  1:739296/150‑1 (MQ=255)                                                                                                  ataCGCGCGGGGATAAAATTTACGGC‑GATTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCtt                                            <  1:986599/151‑1 (MQ=255)                                                                                                      gcgcgGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAAc                                        <  1:953378/151‑1 (MQ=255)                                                                                                        gcgGGGATAAAATTTACGGC‑GGTTATACCATTGATACTTTTTTTGTCA‑CCTATCCCAAAAAAGAAGCCGACGAGGCAAAAGGAAAAATAGTTTTCTCATCAACCGTTATATTTATTTTTTTAAATAGCGGGTTGAAttttgtttttt                                           >  2:874481/1‑138 (MQ=255)                                                                                                         cgGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGc                                     <  1:514756/151‑1 (MQ=255)                                                                                                          ggggATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGc                                     <  2:784179/150‑1 (MQ=255)                                                                                                              aaaaaaTTTGCGGC‑GGTTATACCTTTGTTCTTTTTTTGGTCC‑CCTACCAAAAAAAAGAAGCCGAGGAGCTAAGAACAAATTTATTTCGTTAACCAGCCTTTATTTTTATTTTTTAAACGGGC‑GTTTAAATTGTTGTCTTCAACCGCCATga                                <  2:568349/149‑1 (MQ=255)                                                                                                                                  taCCATTGATCTTTTATTGGTCA‑CCTATCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTTTTTTTTAAACGGGC‑GGTTAAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCcac              <  1:437106/150‑1 (MQ=255)                                                                                                                                      aTTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgcgc         >  1:281105/1‑151 (MQ=255)                                                                                                                                       ttactCCTTTATTGATAA‑CCTATCAAAAAAAAGAAGCCGACGAGTAAAGAGCAAAATTAGTTCTCTAATCAGCCTTTATATTTATTTTTTAAACGGGC‑GTTTAAATTGTTTTCTACAACCGCCATGAGAAACTCACTCTGCCCACAg            <  1:3329/143‑1 (MQ=255)                                                                                                                                       ttGATCTTTTATTGGTCC‑CCTACCAAAAAAAAGAACCCGACGAGAAAAGAGAAAAATTATTGCGCTACACAGCCGTTATATTTATTTTTAAAACTGGG‑GGTTGAATGGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgc           <  2:885299/148‑1 (MQ=255)                                                                                                                                       ttGATCCTTTATGGGTCC‑CCTACCAAAAAGAAGAAGCCGACGAGCAAAGAGAAAAATTAGTTCGCTAATCAGCCTTTATATTTATTTTTAAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAgcgc         <  1:184615/150‑1 (MQ=255)                                                                                                                                          aTCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAg      >  1:143244/1‑150 (MQ=255)                                                                                                                                               ttATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGcc  >  1:193654/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     TTTATTGGCGTAATCAGCAATGACCCTGGTATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGC‑GGTTATACCATTGATCCTTTATTGGTCA‑CCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGC‑GGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCC  >  NC_000913/2151446‑2151735

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