Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,865,253	0	T	G	42.2%	74.1 / 12.7	45	G51G (GGT→GGG)	ygbN	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (6/20);  new base G (19/0);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 7.46e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACAAACGAAC  >  NC_000913/2865109‑2865402                                                                                                                                                 |                                                                                                                                                       aaTTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGGGGGGaaa                                                                                                                                                     >  1:768583/1‑148 (MQ=255)     aCATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGtgatga                                                                                                                                              <  1:779146/151‑1 (MQ=255)      cATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGTgaa                                                                                                                                                     >  2:453569/1‑143 (MQ=255)          gTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGc                                                                                                                                          <  1:152/150‑1 (MQ=255)            tATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGAGGGGGAAGAGGGGATTGCCgg                                                                                                                                       >  1:746767/1‑151 (MQ=255)            tATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGGGGGGGGGAAGGGATGATTGcc                                                                                                                                         >  2:901481/1‑149 (MQ=255)             aTGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGTGGGGGAAAGGATGATTGGCggg                                                                                                                                      >  1:224514/1‑151 (MQ=255)                   tGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAGGGGGGTGAAAAGGTGATAGGCGGGAg‑ggg                                                                                                                                >  1:568150/1‑150 (MQ=255)                      gTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGAAGTGGGGGAAGGGATGATTGGCGGGGT‑ggg                                                                                                                                >  1:455186/1‑147 (MQ=255)                               cGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGGGGGGGAAGGGGTGATTGCCggggagggggggg                                                                                                                            >  2:326540/1‑133 (MQ=255)                                  aaTCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTcc                                                                                                                <  2:479705/151‑1 (MQ=255)                                    tCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGTGGGGGAAGAGATGATTGCCGGGGT‑gggggggggggtttggt                                                                                                                  >  2:807647/1‑136 (MQ=255)                                     cATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTc                                                                                                             <  2:411470/151‑1 (MQ=255)                                           gTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAGGGAAGATTGGCGGGGT‑GGGGGGGGGGGGTGGGTTCGTGACtatt                                                                                                       >  1:466896/1‑151 (MQ=255)                                            ttGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGGGCGGGGGAAGGGAAGATTGCCCGGGG‑GGGGGGTGGGCCTTGTTCTGTccgtc                                                                                                         >  2:331369/1‑144 (MQ=255)                                             tgctgcTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGGGGGGAAAAGGATGGTCGCCGGGGTGGGGGGGGGGGCTTGTTTCCgccccatttt                                                                                                      >  1:173181/1‑142 (MQ=255)                                              gctgctGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACtattat                                                                                                     <  2:958239/150‑1 (MQ=255)                                                 gctgGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACtattatta                                                                                                   <  2:193918/149‑1 (MQ=255)                                                    ggTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGAAGGGGGGAAAGGGATGATCGCCGGGGG‑GGGGGGGGGGGTTGGGTTCGTCCCTatttttttttg                                                                                               >  1:405617/1‑148 (MQ=255)                                                       catcaAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGTGGGGGAAGGGATGATAGCCGGGGT‑GGGGGGGGGGCTTTGGTTCGTCACTATTTTTTTTTGTCt                                                                                            >  1:829068/1‑150 (MQ=255)                                                       catcaAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCt                                                                                            <  2:2944/150‑1 (MQ=255)                                                                 aaGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAAAGAAGATCGCCGGGGA‑GGGGGGGGTGGTTGGGTCCGTCACTCTTTTTTTTTGGCTGGgcg                                                                                       >  1:525499/1‑145 (MQ=255)                                                                 aaGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGc                                                                                 <  1:264944/151‑1 (MQ=255)                                                                        aaCCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGGAAGGGATTATCGCCCGGGG‑GGGGGGGGTGGTTTGTTTCCTCACTATTTTTTTTTGTGTGGGGGCTTTGCTgggggg                                                                          >  1:585306/1‑149 (MQ=255)                                                                           cATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTAt                                                                       <  1:294598/151‑1 (MQ=255)                                                                              tCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGa                                                                     <  2:380425/150‑1 (MQ=255)                                                                                 ttgctttgctCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAGGGATGATCGCCGGGGT‑GGGCGGGGGGGTTTGTTTCGTCCCTATTATTTTTTGGCTGGGGGCTATGCTGGGGGGgatgat                                                                    >  1:510543/1‑148 (MQ=255)                                                                                   gctttgctCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGaa                                                                <  1:161075/150‑1 (MQ=255)                                                                                              cGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCtggtg                                                      <  1:620587/149‑1 (MQ=255)                                                                                              cGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTggtggt                                                    <  1:311563/151‑1 (MQ=255)                                                                                                      tGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAGAGGTGATCGCCGGGGT‑GGGCGGGGGGGTTGGTTCCGTCCCTATTTTTTTTGGTCTGGGGGCTATGCTGGGGGGGATGATCGAACACTCTGGTGGGGGagag                                              >  1:105168/1‑149 (MQ=255)                                                                                                       gTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGt                                             <  2:671349/149‑1 (MQ=255)                                                                                                                  ttGCGGCAGGTATACCGGCGGGCGGAGGGGGGGAAGGGGTGTTAGCCGGGGG‑GGGGGGGGTGGTTTGGTTCGTCCCTCTTTTTATTTGTTTGGGGGCTATGGTGGGGGGGGTGATCGAAAACCCTTGGGGGGGAGAGGCa                                           >  2:376990/1‑140 (MQ=255)                                                                                                                    gCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAattatt                              <  2:250593/151‑1 (MQ=255)                                                                                                                    gCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAattatt                              <  1:453622/151‑1 (MQ=255)                                                                                                                         aGGTATACCGGCGGGCGAAGGGGGGGAAGTGATGATCGCCGGGGG‑GGGCGGGGGGCTTTGGTCCGTCACTATTTTTTTTGGGCTGGGGGGTATGCTGGGGGGGATGATCGAACACTCTGGGGGGGCAGAGTCACTGGCTAATTAtttta                           >  2:317647/1‑147 (MQ=255)                                                                                                                               aCCGGCGGGGGGAGGGGGGAAAGGGATGGTCTCCGGGGT‑GGGCGGGGGGCTTTGTTTCGTCACTATTATTTTTTGTCTTGGGGGCATGCTGGGGGGTATGATCGAACACTCTGGTGGTGCAGAGTCACCTGCGTATTAttttaatc                        >  2:791678/1‑141 (MQ=255)                                                                                                                                  ggcgggcgAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGtt                 >  2:968516/1‑150 (MQ=255)                                                                                                                                   gcgggcgAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGtt                 >  1:130946/1‑149 (MQ=255)                                                                                                                                      ggcgAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTg            <  1:447598/151‑1 (MQ=255)                                                                                                                                       gcgAAGCGGCTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGt             >  2:319194/1‑149 (MQ=255)                                                                                                                                         gAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACa         <  1:758152/151‑1 (MQ=255)                                                                                                                                           aGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACaa        >  1:342541/1‑150 (MQ=255)                                                                                                                                              gggTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACAAACGa    <  2:905904/151‑1 (MQ=255)                                                                                                                                                gTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACAAACGAAc  <  2:260067/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                       AATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACAAACGAAC  >  NC_000913/2865109‑2865402

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