Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,309,012	0	A	G	20.0%	81.3 / 15.4	40	intergenic (‑88/+21)	nlpI/pnp	lipoprotein involved in osmotic sensitivity and filamentation/polynucleotide phosphorylase/polyadenylase
Reads supporting (aligned to +/- strand):  ref base A (16/16);  new base G (4/4);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAC  >  NC_000913/3308865‑3309149                                                                                                                                                    |                                                                                                                                          aCTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAg                                                                                                                                          >  6:137063/1‑149 (MQ=255)       tCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAgg                                                                                                                                    <  2:211434/149‑1 (MQ=255)          tGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGc                                                                                                                                  >  8:64438/1‑148 (MQ=255)                           tCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCg                                                                                                                >  7:120229/1‑149 (MQ=255)                           tCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCg                                                                                                                >  3:173920/1‑149 (MQ=255)                               aaCGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCt                                                                                                            <  8:120229/149‑1 (MQ=255)                                   aaaCACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGcttttaa                                                                                                                            <  2:106893/129‑7 (MQ=255)                                   aaaCACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGcttttaa                                                                                                                            >  1:106893/1‑123 (MQ=255)                                   aaaCACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTc                                                                                                        <  5:137063/149‑1 (MQ=255)                                         cAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCATCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCagcagc                                                                                                  >  7:120567/1‑149 (MQ=255)                                               aaaaaaGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTAc                                                                                                                   <  8:160240/126‑1 (MQ=255)                                               aaaaaaGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTAc                                                                                                                   >  7:160240/1‑126 (MQ=255)                                                  aaaGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTc                                                                                         <  4:173920/149‑1 (MQ=255)                                                        ttCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgc                                                                                   >  6:91012/1‑149 (MQ=255)                                                            tttCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgctg                                                                                 >  6:143070/1‑147 (MQ=255)                                                               cccACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTgcagc                                                                             <  4:105750/148‑1 (MQ=255)                                                                   cTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGtt                                                                        >  7:36659/1‑149 (MQ=255)                                                                   cTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGtt                                                                        >  2:24817/1‑149 (MQ=255)                                                                      ccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                      <  1:24817/148‑1 (MQ=255)                                                                      ccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                      <  7:64438/148‑1 (MQ=255)                                                                         aaGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCaa                                                                                                                         >  2:16587/1‑94 (MQ=255)                                                                         aaGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCaa                                                                                                                         <  1:16587/94‑1 (MQ=255)                                                                         aaGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTAc                                                                                                                   <  3:79134/100‑1 (MQ=255)                                                                         aaGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTAc                                                                                                                   >  4:79134/1‑100 (MQ=255)                                                                              cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCgggggg                                                                                                                                   >  7:115304/1‑79 (MQ=38)                                                                              cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCgggggg                                                                                                                                   <  8:115304/79‑1 (MQ=38)                                                                              cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCaaa                                                                                                                               <  2:161402/83‑1 (MQ=255)                                                                              cACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCaaa                                                                                                                               >  1:161402/1‑83 (MQ=255)                                                                               aCGGTTGAATGAACGTCCTGTTCCCGGTTGCTATCAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTca                                                                                                       <  1:87924/106‑1 (MQ=255)                                                                               aCGGTTGAATGAACGTCCTGTTCCCGGTTGCTATCAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTca                                                                                                       >  2:87924/1‑106 (MQ=255)                                                                                   ttGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGc                                                        <  5:143070/149‑1 (MQ=255)                                                                                        tgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGctt                                                                                                                                <  6:58032/72‑3 (MQ=37)                                                                                        tgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGctt                                                                                                                                >  5:58032/1‑70 (MQ=37)                                                                                            cGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAAt                                               <  8:36659/149‑1 (MQ=255)                                                                                                                 cAAGGCATCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTg                          <  8:120567/149‑1 (MQ=255)                                                                                                                  aaGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTgg                         >  4:132363/1‑149 (MQ=255)                                                                                                                                    aaaaGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACt       >  3:53352/1‑149 (MQ=255)                                                                                                                                     tttgCCCTCCGCTGCGGCGGGGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTcc             <  4:56361/139‑1 (MQ=255)                                                                                                                                     tttgCCCTCCGCTGCGGCGGGGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTcc             >  3:56361/4‑142 (MQ=255)                                                                                                                                         ccccccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAc  <  5:91012/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                          ACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAC  >  NC_000913/3308865‑3309149

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