Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,946,790	T→G	17.0%	G180G (GGT→GGG)	yebB →	DUF830 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,946,790	0	T	G	17.0%	135.4 / 27.1	66	G180G (GGT→GGG)	yebB	DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base T (22/32);  new base G (4/7);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

GCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGC  >  NC_000913/1946698‑1946888                                                                                             |                                                                                                   gCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgtgt                                                                                             >  1:609204/1‑100 (MQ=255) gCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgt                                                                                               >  1:1444862/1‑98 (MQ=255)  cAAAACTCACTTTCTGGAAATTCTGGTTCTTAGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGtt                                                                                            >  2:2312652/1‑100 (MQ=255)  cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGtt                                                                                            >  1:254749/1‑100 (MQ=255)  cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGtt                                                                                            >  1:1005767/1‑100 (MQ=255)     aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAt                                                                                         <  1:1143486/100‑1 (MQ=255)      aCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGAtt                                                                                        <  2:234969/100‑1 (MQ=255)       gtCACTTTCTGGAGATTTTGGTTCGTGGGGTCTATTCTGTGGGAGCGTAAAACGGGCACGTCATCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTc                                                                                       <  2:2263040/99‑1 (MQ=255)       cTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAtt                                                                                        <  1:1916205/99‑1 (MQ=255)           cTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGc                                                                                   <  2:860706/100‑1 (MQ=255)           cTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGc                                                                                   <  1:1273067/100‑1 (MQ=255)             ttCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCg                                                                                  >  2:1710780/1‑99 (MQ=255)                tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCggagg                                                                              <  2:1603559/100‑1 (MQ=255)                 ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCggaggg                                                                             >  2:1518235/1‑100 (MQ=255)                  gAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCggaggg                                                                             <  1:584205/99‑1 (MQ=255)                    aaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGt                                                                          >  1:1250471/1‑100 (MQ=255)                     aTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTg                                                                         <  2:1875531/100‑1 (MQ=255)                     aTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTg                                                                         <  1:2190105/100‑1 (MQ=255)                      ttCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCTTCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTg                                                                         <  2:744282/99‑1 (MQ=255)                      ttCGGGTTCCTGGTTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGGCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTgg                                                                        <  2:2311388/100‑1 (MQ=255)                        cTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGaa                                                                      >  2:1032381/1‑100 (MQ=255)                        cTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGaa                                                                      <  2:527943/100‑1 (MQ=255)                          ggTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGaaag                                                                    >  1:527946/1‑99 (MQ=255)                           gTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACg                                                                   <  1:565586/100‑1 (MQ=255)                             tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcn                                                                 <  1:1283427/100‑2 (MQ=255)                             tCCTGGGTTCTATTCCGTCGGAGCTTAAAACCGTTCCGCCCCCCAGGTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACg                                                                   <  2:1769060/98‑1 (MQ=255)                                tGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCt                                                                <  1:436766/98‑1 (MQ=255)                                tGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCa                                                              <  1:1404419/100‑1 (MQ=255)                                 gggTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAg                                                             >  1:2361668/1‑100 (MQ=255)                                   gTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGcc                                                           <  2:2176617/100‑1 (MQ=255)                                   gTTCTATTCCGTGGGAGCGTAAAACCGTAACGCCAGCAAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGcc                                                           >  2:1452128/1‑100 (MQ=255)                                        aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                      >  1:1320492/1‑100 (MQ=255)                                           ccGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGa                                                   >  2:591507/1‑100 (MQ=255)                                                 gAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAgg                                             <  2:1005756/100‑1 (MQ=255)                                                  aGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGa                                               >  1:263869/1‑97 (MQ=255)                                                  aGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCACGCGGTGGGAGTTGAAACGCCTCAGCCTGAACTGACCGAGGc                                            <  2:1289360/100‑1 (MQ=255)                                                  aGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATACGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGc                                            >  2:1168838/1‑100 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                          <  1:1986461/100‑1 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCAAGTTTGTGGCATCATCAGGGCTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGACTCAGCCGGAACTGACCGACGc                                            >  2:1677080/1‑98 (MQ=255)                                                      tAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                          <  1:1921535/98‑1 (MQ=255)                                                        aaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtata                                      <  1:907706/100‑1 (MQ=255)                                                          aCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAAc                                    <  1:574755/100‑1 (MQ=255)                                                           cgggtaggccaggcgGTTTGTGCGATGATGCGGGTTTGAGGTTGATTCAGGCGGAGGGAGAGGAAACGCCTCAGCCGGAACTGACAGAGGCGGTATAAc                                    <  2:2060550/92‑1 (MQ=255)                                                            cGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACtt                                  >  2:1873022/1‑100 (MQ=255)                                                                aCGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTTACCGAGGCGGTATAACTTAACg                              <  2:846273/100‑1 (MQ=255)                                                                 cGCCTGCCAGTTTGTGGCATCATTCGGGTTTGGTGTTGATTCATGCGGAGGGAGGGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAAc                                    <  2:2070580/93‑1 (MQ=255)                                                                 cGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGc                             >  1:383636/1‑100 (MQ=255)                                                                 cGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGc                             >  2:129486/1‑100 (MQ=255)                                                                 cGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGc                             >  2:1823187/1‑100 (MQ=255)                                                                  gtcagccagTTTTTTGCATCATTCGGGTCTTGTTCTGATTCATGCGGAGGGAGTGTAAACGCCTCCGTCGGTATTGATCGAGGTGGTATAACTTACCGc                             <  2:1435041/97‑1 (MQ=255)                                                                  gccagccaGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCa                            >  1:1562041/1‑100 (MQ=255)                                                                    cagccagTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGt                          >  1:376529/1‑100 (MQ=255)                                                                     agccagTTTGTGGCATCATCCGGGTTTGGTGTGGATTCACGCTGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTc                         <  2:1683646/100‑1 (MQ=255)                                                                       ccagTTTGTGGCATCATCCGGGTTGGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGc                       <  2:1339060/100‑1 (MQ=255)                                                                        cagTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                      >  2:2075641/1‑100 (MQ=255)                                                                        cagTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                      <  1:2257678/100‑1 (MQ=255)                                                                              gtgGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg                <  2:1569281/100‑1 (MQ=255)                                                                                  catcatCGGGGTTTGGTGTTGATTCATGCGGATGGAGTGGAAACGCCTCCGGCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAg            <  2:1427782/100‑1 (MQ=255)                                                                                  catcatCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCa             >  2:530034/1‑99 (MQ=255)                                                                                  catcatCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAg            >  2:1796546/1‑100 (MQ=255)                                                                                   atcatcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg           <  1:237011/100‑1 (MQ=255)                                                                                    tcatcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGt          <  1:877600/100‑1 (MQ=255)                                                                                    tcatcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGt          <  2:132191/100‑1 (MQ=255)                                                                                     catcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGtt         <  1:473854/100‑1 (MQ=255)                                                                                      atcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGtt         <  2:2388785/99‑1 (MQ=255)                                                                                           ggTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGcc    >  1:22970/1‑99 (MQ=255)                                                                                             tttGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcgc  <  2:952475/99‑1 (MQ=255)                                                                                             |                                                                                                   GCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCGC  >  NC_000913/1946698‑1946888

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