Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	145,758	0	G	T	22.0%	117.8 / 27.9	50	K226N (AAG→AAT)	yadE	putative polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (23/16);  new base T (8/3);  total (31/19)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACAATATGGCATGAAG  >  NC_000913/145661‑145857                                                                                                  |                                                                                                    caTCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGAt                                                                                                   >  2:601413/1‑100 (MQ=255)   tCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCa                                                                                                 >  1:752514/1‑100 (MQ=255)     ggTACGCGCTTTCAATAACCAGATGGCCTGGCTGCATGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCAAt                                                                                               >  1:1200740/1‑100 (MQ=255)         cgcgCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTcc                                                                                           >  2:1117308/1‑100 (MQ=255)            gcTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCAATCTCCCTg                                                                                        >  2:1051552/1‑100 (MQ=255)               ttCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCAATCTCCCTgcgc                                                                                     >  2:375196/1‑100 (MQ=255)               ttCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTgcgc                                                                                     >  2:264010/1‑100 (MQ=255)                tCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGcgcg                                                                                    <  2:68221/100‑1 (MQ=255)                tCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGcgcg                                                                                    <  2:422521/100‑1 (MQ=255)                   aTAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGc                                                                                 >  2:1092255/1‑100 (MQ=255)                      aCCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggt                                                                              >  2:567206/1‑100 (MQ=255)                         aGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggtggt                                                                           >  2:949765/1‑100 (MQ=255)                         aGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCggtggt                                                                           <  2:928982/100‑1 (MQ=255)                            tGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCAATCTCCCTGCGCGAGCGGTGGTGAt                                                                        >  1:1184560/1‑100 (MQ=255)                             ggCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGAt                                                                        <  2:109684/99‑1 (MQ=255)                                  ggCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCtt                                                                  >  1:1239353/1‑100 (MQ=255)                                   gCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCttt                                                                 <  2:484608/100‑1 (MQ=255)                                      gCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTtgat                                                              <  1:1044500/100‑1 (MQ=255)                                       cGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTtgatg                                                             <  1:477515/100‑1 (MQ=255)                                         tGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTgatgat                                                           >  2:218602/1‑100 (MQ=255)                                           aCAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATgg                                                         <  2:518084/100‑1 (MQ=255)                                            cAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGc                                                        >  1:56085/1‑100 (MQ=255)                                               gggATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTc                                                     >  1:1154064/1‑100 (MQ=255)                                                ggATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAATATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCa                                                    <  2:752509/100‑1 (MQ=255)                                                ggATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCa                                                    <  2:739856/100‑1 (MQ=255)                                                ggATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCa                                                    >  1:1143052/1‑100 (MQ=255)                                                ggATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCa                                                    >  1:467592/1‑100 (MQ=255)                                                     cgcgACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCg                                               >  2:407393/1‑100 (MQ=255)                                                         acacTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGa                                           <  1:1279917/100‑1 (MQ=255)                                                          cacTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAg                                          >  1:813240/1‑100 (MQ=255)                                                            cTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAgtc                                        <  1:1179760/100‑3 (MQ=255)                                                            cTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGcc                                        >  1:371286/1‑100 (MQ=255)                                                            cTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGcc                                        <  1:1081808/100‑1 (MQ=255)                                                                gCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTa                                    >  2:300299/1‑100 (MQ=255)                                                                  atgggtcAGCTGGAAGGCTACGTGAAGAATAAGATCAAGCTCCCTGCGCGAGCGGGGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTAt                                   <  2:1134984/93‑1 (MQ=255)                                                                      tGCAGCTGGAAGGCTACGTGAAGAATAATATCAATCTCCCTGCGCGAGCGGTGGTGATTAACTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTa                              >  2:245414/1‑100 (MQ=255)                                                                         aGCTGCAAGGCTACGTGAAGAATAAGATCAAGCTCCCTGCGCGAGCGGTGGTGATTAACTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATcc                           >  2:509803/1‑100 (MQ=255)                                                                          gCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCt                          >  1:403680/1‑100 (MQ=255)                                                                             ggAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCtgtg                       <  1:1092265/100‑1 (MQ=255)                                                                               aaGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTgtgt                      <  1:407397/99‑1 (MQ=255)                                                                                aGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTg                    >  1:519702/1‑100 (MQ=255)                                                                                  gCTACGTGAAGAATAATATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGaa                  <  2:577109/100‑1 (MQ=255)                                                                                  gCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCTTGTGTTGaa                  >  2:867968/1‑100 (MQ=255)                                                                                  gCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGaa                  >  2:1182822/1‑100 (MQ=255)                                                                                   cTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGaaa                 <  1:1117318/100‑1 (MQ=255)                                                                                    tACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAAc                <  1:169851/100‑1 (MQ=255)                                                                                     aCGTGAAGAATAATATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACa               >  1:1229209/1‑100 (MQ=255)                                                                                     aCGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACa               >  1:938441/1‑100 (MQ=255)                                                                                           agaaTAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACAATATgg         >  1:673395/1‑100 (MQ=255)                                                                                              aTAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACAATATGGCAt      >  1:706503/1‑100 (MQ=255)                                                                                                  gATCAATCTCCCTGCGCGAGCGGTGGTGATTACGTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACAATATGGCATGAAg  >  2:855377/1‑100 (MQ=255)                                                                                                  gATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACAATATGGCATGAAg  <  2:371283/100‑1 (MQ=255)                                                                                                  |                                                                                                    CATCGGTACGCGCTTTCAATAACCAGATGGCCTGGCTGCGTGACAGGGGATACGCGACACTGAGCATGGTGCAGCTGGAAGGCTACGTGAAGAATAAGATCAATCTCCCTGCGCGAGCGGTGGTGATTACCTTTGATGATGGCCTCAAGTCGGTGAGCCGCTATGCGTATCCTGTGTTGAAACAATATGGCATGAAG  >  NC_000913/145661‑145857

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