Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	53.8%	47.8 / 63.4	39	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (16/2);  new base C (0/21);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 4.06e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-07

TTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATT  >  NC_000913/501431‑501723                                                                                                                                                  |                                                                                                                                                     ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                                                                                 >  1:967866/1‑150 (MQ=255)        cgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCgg                                                                                                                                         <  2:718423/151‑1 (MQ=255)                         agggCTTTCCCGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATCCCTACA‑TATTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGtt                                                                                                                        <  1:188243/149‑1 (MQ=38)                              gTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGc                                                                                                                   >  2:207957/1‑151 (MQ=255)                                 cACGCCGCATCTGGCTTTCGCCTCTCGTTGCTTGAGCCGCCGCTGCCGCGTTTTATCAT‑CCTTCATTTTTTTTCATATTTTACACCCGACACCCCCCTTTTCCCCGGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATg                                                                                                                <  1:438750/151‑1 (MQ=255)                                    gccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                              >  2:315053/1‑150 (MQ=255)                                             ggCTTTTGCCTCTCGATGCCTGATGCGCCGCTGGCGCGTTTTATCATGCCTACA‑TTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                                       <  2:866656/148‑1 (MQ=255)                                               cttttGGCTCTCGATGCCGGATGCGACGCTGGCGCGTTTTATCTTGCCTACA‑TATTTTTCATTTTTTACACCCGGCAACCCCCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                                       <  2:539744/144‑1 (MQ=255)                                               cATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                                  >  2:148604/1‑151 (MQ=255)                                                      ctctCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATTTTTTACATCCGGCACCCACCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTc                                                                                            <  2:179776/150‑1 (MQ=39)                                                        ctcGATGCCTGATGCGACGCTGGCGCTTCTTATCATCCCTACA‑TTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                              <  1:697899/146‑1 (MQ=255)                                                          cGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTACATATTTTACACCCGGCAACCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                          <  2:522340/148‑1 (MQ=255)                                                           ggggCCTGGTGCGACGCGGGCGCGTCTTATCATACTTATA‑TATTTTTAATTTTTTACCGCCGACAGCCGCCGTTTCCCCCGGCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                                                                                           <  1:335122/143‑1 (MQ=255)                                                           gATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                                      >  2:724144/1‑151 (MQ=255)                                                           gATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                                                                                      >  2:927201/1‑151 (MQ=255)                                                            aTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATc                                                                                     >  2:184475/1‑151 (MQ=255)                                                                             ggCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                         <  2:167165/146‑1 (MQ=255)                                                                               cgcgTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcaccac                                                                  >  2:920670/1‑151 (MQ=255)                                                                                 cgccTTATCATGCCTACA‑TATTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                                                 <  1:26644/147‑1 (MQ=255)                                                                                 cgTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                                >  1:727343/1‑151 (MQ=255)                                                                                  gTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                               >  1:795414/1‑151 (MQ=255)                                                                                   tCTTATCATGCCTACA‑TTTTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                               <  1:929833/150‑1 (MQ=255)                                                                                    cTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                             >  1:483546/1‑151 (MQ=255)                                                                                      tATCATGCCTACT‑TTTTTTTCATTTTTTACACCCGCCACCCCCCTTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg                                                            <  1:275384/150‑1 (MQ=255)                                                                                        tCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                                          >  2:111459/1‑150 (MQ=255)                                                                                             ccTACA‑TATTTTTCATATTTTACATCCGGCAACCCCCGTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                      <  2:891495/149‑1 (MQ=255)                                                                                                  a‑tatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                               >  1:387416/1‑151 (MQ=255)                                                                                                      tttttCATTTTTTACACCCGCCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACtt                                             <  1:777202/150‑1 (MQ=255)                                                                                                        tttCATTTTTTACCCCCGGCACCCCCCGTTTACCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCa                                           <  1:396340/150‑1 (MQ=255)                                                                                                        ttccATATTTTCCACCCGGCACCCCCCGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACGACCTGATGCGCACCACGGGCAGGGAGAGACGCCACTTCa                                           <  2:812875/147‑1 (MQ=255)                                                                                                           cATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                                       >  2:419005/1‑151 (MQ=255)                                                                                                            atatTTTACACCCGGCAACCCCCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                                         <  2:514446/148‑1 (MQ=255)                                                                                                               ttttACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                                          <  2:961004/144‑1 (MQ=255)                                                                                                                   aCATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                               >  2:947051/1‑151 (MQ=255)                                                                                                                       ccGCCAACCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                                 <  1:631971/145‑1 (MQ=255)                                                                                                                         ggCACCCACCGTTTCCCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGcgcg                          <  2:250343/150‑1 (MQ=255)                                                                                                                         ggCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg                         <  1:778788/151‑1 (MQ=255)                                                                                                                                     ttcccccGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtctc              <  2:283111/147‑1 (MQ=255)                                                                                                                                             tCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCgg     >  2:21860/1‑151 (MQ=255)                                                                                                                                              caccacCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCgg     >  1:526922/1‑150 (MQ=255)                                                                                                                                                ccaccTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAtt  <  2:487316/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     TTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACA‑TATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATT  >  NC_000913/501431‑501723

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