Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,688	0	T	A	20.0%	90.1 / 20.1	45	intergenic (+31/+33)	lpp/ldtE	murein lipoprotein/murein L,D‑transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base A (4/5);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGA  >  NC_000913/1757542‑1757832                                                                                                                                                   |                                                                                                                                                 ttGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATtg                                                                                                                                                <  4:153110/149‑1 (MQ=255)     ccAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTgcg                                                                                                                                            >  4:30871/1‑149 (MQ=255)      cAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTgcgc                                                                                                                                           >  2:144484/1‑149 (MQ=255)        gCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                         <  1:168704/149‑1 (MQ=255)              cAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttt                                                                                                                                   <  1:144484/149‑1 (MQ=255)                  gacgTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGt                                                                                                                               <  3:127987/149‑1 (MQ=255)                            cAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAtt                                                                                                                                       <  1:200934/131‑1 (MQ=255)                            cAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAtt                                                                                                                                       >  2:200934/1‑131 (MQ=255)                             aaTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTAc                                                                                                                    <  8:133906/149‑1 (MQ=255)                             aaTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTAc                                                                                                                    >  7:133906/1‑149 (MQ=255)                                   ttCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAtttttca                                                                                                                                  >  1:57525/1‑127 (MQ=255)                                   ttCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAtttttca                                                                                                                                  <  2:57525/129‑3 (MQ=255)                                      cGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTg                                                                                                                            >  2:140925/1‑132 (MQ=255)                                      cGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTg                                                                                                                            <  1:140925/132‑1 (MQ=255)                                           ttCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCAGTTTACCGCTACTGCGTCAc                                                                                                      <  8:68104/149‑1 (MQ=255)                                           ttCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAt                                                                                                                                        <  5:97840/115‑1 (MQ=255)                                           ttCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCAt                                                                                                                                        >  6:97840/1‑115 (MQ=255)                                               ggCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCAcgcg                                                                                                   >  4:163398/1‑148 (MQ=255)                                                    ctAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttttt                                                                                                                                 >  2:142058/1‑113 (MQ=255)                                                    ctAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttttt                                                                                                                                 <  1:142058/113‑1 (MQ=255)                                                       aaGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttt                                                                                                                                   >  2:136094/1‑108 (MQ=255)                                                       aaGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAtttttt                                                                                                                                   <  1:136094/108‑1 (MQ=255)                                                       aaGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATAtt                                                                                          >  5:123884/1‑149 (MQ=255)                                                          aTGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTccc                                                                                       <  6:123884/149‑1 (MQ=255)                                                           tGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCt                                                                                      <  3:30871/149‑1 (MQ=255)                                                             aCGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                         <  2:65744/96‑1 (MQ=255)                                                             aCGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                         >  1:65744/1‑96 (MQ=255)                                                                   cTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTgg                                                                              <  3:163398/149‑1 (MQ=255)                                                                           aaCCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTg                                                                               <  1:31122/140‑1 (MQ=255)                                                                           aaCCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTg                                                                               >  2:31122/1‑140 (MQ=255)                                                                              cAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTccc                                                                                       <  4:85124/129‑1 (MQ=255)                                                                              cAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTccc                                                                                       >  3:85124/1‑129 (MQ=255)                                                                                   tCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTa                                                                                                                     <  8:68273/94‑1 (MQ=255)                                                                                   tCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTa                                                                                                                     >  7:68273/1‑94 (MQ=255)                                                                                       gACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGAct                                                          >  6:12687/1‑149 (MQ=255)                                                                                             aTGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACt                                                     >  1:99459/1‑148 (MQ=255)                                                                                                    ctaAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGgc                                               <  5:12687/147‑1 (MQ=255)                                                                                                             gCACGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAAc                                                                                               <  5:173558/90‑1 (MQ=255)                                                                                                             gCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAAc                                                                                               >  6:173558/1‑90 (MQ=255)                                                                                                                                  gtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCgg                               <  2:162472/133‑1 (MQ=255)                                                                                                                                  gtgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCgg                               >  1:162472/1‑133 (MQ=255)                                                                                                                                   tgaaAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTg                                                                                    <  8:153925/79‑1 (MQ=255)                                                                                                                                   tgaaAAATGGCGCACATAGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTg                                                                                    >  7:153925/1‑79 (MQ=255)                                                                                                                                              gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGACACCg    >  8:66679/1‑148 (MQ=255)                                                                                                                                               cccaTGGTGCGCCATTTTTTTTGTCTGCCGTGTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCGGCGGGAGGTGCTCCACTGCTCACCGAAACCGGa  <  7:66679/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 TTGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGGA  >  NC_000913/1757542‑1757832

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