Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	7,925	+A	coding (35/1431 nt)	yaaJ ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	7,922	1	.	A	95.7%	143.0 / ‑0.7	46	coding (38/1431 nt)	yaaJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (19/25);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCC‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCT  >  NC_000913/7792‑8055                                                                                                                                     |                                                                                                                                       gCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCC‑a‑aaaggac                                                                                                                               <  2:964584/139‑6 (MQ=255) gCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCC‑a‑aaaggac                                                                                                                               <  2:55834/139‑6 (MQ=255) gCCTGGCTGTGGATGAAAGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCC‑a‑aaaggac                                                                                                                               <  2:419506/139‑6 (MQ=255)  ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACg                                                                                                                               <  2:471805/139‑1 (MQ=255)  ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACg                                                                                                                               <  1:792102/139‑1 (MQ=255)  ccTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACg                                                                                                                               >  2:382783/1‑139 (MQ=255)   cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGc                                                                                                                              <  1:382783/139‑1 (MQ=255)   cTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGc                                                                                                                              >  1:878952/1‑139 (MQ=255)                aaTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTA‑‑GAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGgaga                                                                                                               >  1:298827/1‑139 (MQ=255)                 aTGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGgag                                                                                                                >  2:528159/1‑139 (MQ=255)                    cTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGaaa                                                                                                             >  2:314806/1‑139 (MQ=255)                         tttAAGACTTTTGCCAAACAGGCGGATGTAGCGAAACTGCACAAATCCGGCGCGAAAAGTGAACCAACAACCTGCGCCGAAGGGCAGGTAAATCATTACCGATCGCCA‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATc                                                                                                        >  2:781464/1‑139 (MQ=255)                            aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                     <  1:720999/139‑1 (MQ=255)                            aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                     >  2:588285/1‑139 (MQ=255)                            aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                     >  1:493601/1‑139 (MQ=255)                            aaGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTgg                                                                                                     <  1:528159/139‑1 (MQ=255)                            aaGACTTTTGCCAAACTGGCGGATGTA‑‑GAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGgca                                                                                                   <  2:298827/139‑1 (MQ=255)                             aGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGgc                                                                                                    <  1:673969/139‑1 (MQ=255)                               aCTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcat                                                                                                  >  2:383305/1‑139 (MQ=255)                                   ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                                                              <  2:209766/139‑1 (MQ=255)                                   ttGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgcat                                                                                              >  1:209766/1‑139 (MQ=255)                                     gCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCatat                                                                                            <  1:743200/139‑1 (MQ=255)                                      ccAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                           <  2:80764/139‑1 (MQ=255)                                              gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                          >  1:349469/1‑132 (MQ=255)                                              gCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                          <  2:349469/132‑1 (MQ=255)                                                     tAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGt                                                                            <  1:690444/139‑1 (MQ=255)                                                        cGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTcgc                                                                         <  2:1047728/139‑1 (MQ=255)                                                          aaaCTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGa                                                                       <  2:205896/139‑1 (MQ=255)                                                               gCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACt                                                                  <  1:398074/139‑1 (MQ=255)                                                                cacaAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAGTCATTACCGATC‑CCC‑ATAAGGACGCTGTTGCTGAAGGAGAAAAAATCTGGCATGCATATGCCTCTTATTGCCGGTCGCGATGACtt                                                                 >  2:766365/1‑139 (MQ=255)                                                                                          aacaacCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGttt                                       <  2:958907/139‑1 (MQ=255)                                                                                            caacCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTATAAACGTTACCAATTGTTTaa                                     <  2:878952/139‑1 (MQ=255)                                                                                              acCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                   <  1:398710/139‑1 (MQ=255)                                                                                                  gcgcCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                          >  2:563595/1‑80 (MQ=255)                                                                                                  gcgcCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATcc                                                                                          <  1:563595/80‑1 (MQ=255)                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                           <  1:268028/78‑1 (MQ=255)                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATc                                                                                           >  2:268028/1‑78 (MQ=255)                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                   >  2:324355/1‑134 (MQ=255)                                                                                                   cgcCGAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaaga                                   <  1:324355/134‑1 (MQ=255)                                                                                                       gAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATAc                          >  1:112088/1‑139 (MQ=255)                                                                                                       gAAGAGCAGGTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATAc                          >  1:576099/1‑139 (MQ=255)                                                                                                               ggTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAg                  >  1:915657/1‑139 (MQ=255)                                                                                                                gTAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAgg                 <  1:88100/139‑1 (MQ=255)                                                                                                                 tAAATCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGGCTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGt                <  1:647773/139‑1 (MQ=255)                                                                                                                   aaTCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTAc              <  1:234821/139‑1 (MQ=255)                                                                                                                     tCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                               <  2:470579/56‑1 (MQ=255)                                                                                                                     tCATTACCGATC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGgcatgca                                                                                               >  1:470579/1‑56 (MQ=255)                                                                                                                             atcc‑cca‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCaa                                             <  2:392103/94‑1 (MQ=255)                                                                                                                             atcc‑cca‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                    <  2:129311/103‑1 (MQ=255)                                                                                                                             atcc‑cca‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTaag                                    >  1:129311/4‑106 (MQ=255)                                                                                                                             atcc‑cca‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc   >  2:219085/4‑139 (MQ=255)                                                                                                                             atcc‑cca‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc   >  2:207800/4‑139 (MQ=255)                                                                                                                              tcc‑cca‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCt  <  1:588285/137‑1 (MQ=255)                                                                                                                              aTC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCaa                                             >  1:392103/1‑97 (MQ=255)                                                                                                                              aTC‑CCCAA‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTc   >  1:478232/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                       GCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATC‑CCC‑A‑AAGGACGCTGTTAATGAAGGAGAAAAAATCTGGCATGCATATCCCTCTTATTGCCGGTCGCGATGACTTTCCTGTGTAAACGTTACCAATTGTTTAAGAAGTATATACGCTACGAGGTACTTGATAACTTCT  >  NC_000913/7792‑8055

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