Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	604,753	A→C	13.6%	intergenic (‑90/+18)	ybdG ← / ← nfsB	mechanosensitive channel protein, miniconductance/dihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	604,753	0	A	C	13.6%	124.1 / 15.8	44	intergenic (‑90/+18)	ybdG/nfsB	mechanosensitive channel protein, miniconductance/dihydropteridine reductase, NAD(P)H‑dependent, oxygen‑insensitive
Reads supporting (aligned to +/- strand):  ref base A (22/16);  new base C (4/2);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

TGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAA  >  NC_000913/604606‑604899                                                                                                                                                    |                                                                                                                                                   tggtgTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGAt                                                                                                                                                   <  2:9992/149‑1 (MQ=255)        atcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGc                                                                                                                                            >  1:106516/1‑149 (MQ=255)        atcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGc                                                                                                                                            >  1:536613/1‑149 (MQ=255)           gatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAag                                                                                                                                         >  4:547881/1‑149 (MQ=255)                cAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGAAAgagaga                                                                                                                                    <  3:539547/149‑1 (MQ=255)                 aaTACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGaa                                                                                                                                   <  2:68903/149‑1 (MQ=255)                      ccGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTaca                                                                                                                              >  3:166706/1‑149 (MQ=255)                                ttCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTaa                                                                                                                    >  3:475288/1‑149 (MQ=255)                                 tCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAg                                                                                                                   <  4:166706/149‑1 (MQ=255)                                 tCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAg                                                                                                                   <  3:120001/149‑1 (MQ=255)                                        gggATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGt                                                                                                            >  2:31205/1‑149 (MQ=255)                                              ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGct                                                                                                                                           >  1:74487/1‑111 (MQ=255)                                              ttAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCATCTGCCCGGct                                                                                                                                           <  2:74487/112‑2 (MQ=255)                                                 aaTGCTGCATAGCGGCTCATGTAAACAGCAGCAGCCTATGATTACGGCTGAATCAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCACGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCa                                                                                                   >  2:432562/1‑149 (MQ=255)                                                 aaTCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCa                                                                                                   <  1:31205/149‑1 (MQ=255)                                                      tGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACga                                                                                              >  1:226596/1‑149 (MQ=255)                                                          tAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGAtt                                                                                          >  1:279504/1‑149 (MQ=255)                                                           aaCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGAttt                                                                                         >  1:125979/1‑149 (MQ=255)                                                            aCGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTc                                                                                        >  2:386249/1‑149 (MQ=255)                                                              gTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCgg                                                                                      >  2:502499/1‑149 (MQ=255)                                                                       ataCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAaga                                                                                                                                        <  1:313239/90‑1 (MQ=255)                                                                       ataCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAaga                                                                                                                                        >  2:313239/1‑90 (MQ=255)                                                                        taCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgtagcg                                                                            <  3:547881/149‑1 (MQ=255)                                                                         acaacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCagcgtagcgt                                                                           >  1:187911/1‑149 (MQ=255)                                                                           aacaGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                         >  4:96684/1‑149 (MQ=255)                                                                               gcagcCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAAt                                                                     <  1:502499/149‑1 (MQ=255)                                                                                           gaCGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCtgt                                                         <  3:136602/149‑1 (MQ=255)                                                                                                              tGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGaacaac                                      <  2:279504/149‑1 (MQ=255)                                                                                                                 aaaTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                   >  3:2496/1‑149 (MQ=255)                                                                                                                 aaaTCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccac                                   >  1:365236/1‑149 (MQ=255)                                                                                                                        gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTg                            >  3:157422/1‑149 (MQ=255)                                                                                                                        gagaATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTg                            >  3:300045/1‑149 (MQ=255)                                                                                                                             tCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCAGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTa                       <  4:441105/149‑1 (MQ=255)                                                                                                                                  gagGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCtt                  >  3:550532/1‑149 (MQ=255)                                                                                                                                   agGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCtg                                                          >  3:508282/1‑108 (MQ=255)                                                                                                                                   agGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCtg                                                          <  4:508282/108‑1 (MQ=255)                                                                                                                                   agGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCttt                 >  2:434281/1‑149 (MQ=255)                                                                                                                                    ctctcttgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTtc                >  3:27931/8‑149 (MQ=255)                                                                                                                                       tcttgCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTctct             >  4:579918/5‑149 (MQ=255)                                                                                                                                         tAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCttt           <  2:531349/149‑1 (MQ=255)                                                                                                                                           gCCGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCa         <  4:27931/149‑1 (MQ=255)                                                                                                                                            ccGGGCATCTGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAg        >  4:201260/1‑149 (MQ=255)                                                                                                                                              gggCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGAc      <  1:546170/149‑1 (MQ=255)                                                                                                                                                  aGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCaaa  <  2:365236/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAACGTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAA  >  NC_000913/604606‑604899

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