Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,865,246	0	A	G	42.2%	62.8 / 10.3	45	E49G (GAA→GGA)	ygbN	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (7/19);  new base G (19/0);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACA  >  NC_000913/2865098‑2865395                                                                                                                                                     |                                                                                                                                                       cccATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGaa                                                                                                                                                       <  2:206772/150‑1 (MQ=255)  ccATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGt                                                                                                                                                     <  1:929468/151‑1 (MQ=255)   cATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAggg                                                                                                                                                    >  1:713518/1‑149 (MQ=255)            aaTTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGGGGGGaaa                                                                                                                                              >  1:768583/1‑148 (MQ=255)                aCATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGggggggaaaag                                                                                                                                           >  1:739593/1‑140 (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)                                                                                                                                                     |                                                                                                                                                       CCCATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCGCCGGGAT‑GGGCGGTGTGCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTGACA  >  NC_000913/2865098‑2865395

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