Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	604,759	G→A	intergenic (‑96/+12)	ybdG ← / ← nfsB	miniconductance mechanosensitive channel YbdG/NAD(P)H‑dependent nitroreductase NfsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	604,759	0	G	A	96.9%	101.6 / NA	32	intergenic (‑96/+12)	ybdG/nfsB	miniconductance mechanosensitive channel YbdG/NAD(P)H‑dependent nitroreductase NfsB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/11);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCACCCAATGCAAAATAATATGCACCACGACGGCGGTCAGAAAAATAATACCGAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTG  >  NC_000913/604534‑605001                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   ccaccCAATGCAAAATAATATGCACCACGACGGCGGTCAGAAAAATAATACCGAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGATAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGtt                                                                                                                                                                                                                            >  1:69773/1‑250 (MQ=255)      cAATGCAAAATAATATGCACCACGACGGCGGTCAGAAAAATAATACCGAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTg                                                                                                                                                                                                                      >  1:48732/1‑251 (MQ=255)        aTGCAAAATAATATGCACCACGACGGCGGTCAGAAAAATAATACCGAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGAt                                                                                                                                                                                                                    >  2:7464/1‑251 (MQ=255)                       caccacGACGGCGGTCAGAAAAATAATACCGAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACg                                                                                                                                                                                                     <  1:332687/251‑1 (MQ=255)                         ccacGACGGCGGTCAGAAAAATAATACCGAAAAGCATCACCATCGAGGTGGTGGGATCGAGCTCAATACCCGCGAAATCTTCAACCTGGGATATGAAATCCTGCATAAC‑GGCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGGTAGGCAAATCAGGAGGATGGGAGAGGAAATAGCCGGGCAGATGGCCGGCAAGAGGGGATTGCACGGCGGTTAAGGTGATGTTGTGCTGCAGACgag                                                                                                                                                                                                   <  2:193084/251‑1 (MQ=255)                                              taataCCGAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTa                                                                                                                                                                               >  2:309411/1‑250 (MQ=255)                                                     gAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTc                                                                                                                                                                       >  2:86476/1‑251 (MQ=255)                                                                 cATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGtgatg                                                                                                                                                           >  2:159334/1‑251 (MQ=255)                                                                       ggtggtGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATTCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTGAAAATCTTCAACGCTGTGATGACCTAc                                                                                                                                                     <  1:320542/251‑1 (MQ=255)                                                                                 gcgatcTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAacc                                                                                                                                           <  2:102528/250‑1 (MQ=255)                                                                                     tctcAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccacc                                                                                                                                        >  2:21302/1‑250 (MQ=255)                                                                                        cAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAaccacc                                                                                                                                        <  1:11908/247‑1 (MQ=255)                                                                                                 cTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGc                                                                                                                           >  2:118159/1‑251 (MQ=255)                                                                                                 cTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGc                                                                                                                           >  1:243147/1‑251 (MQ=255)                                                                                                   aaaTCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCt                                                                                                                         <  2:298508/251‑1 (MQ=255)                                                                                                    aaTCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCtt                                                                                                                        <  2:324936/251‑1 (MQ=255)                                                                                                          cAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATCCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCtt                                                                                                                  >  2:18875/1‑251 (MQ=255)                                                                                                             ccTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTc                                                                                                                >  2:308927/1‑250 (MQ=255)                                                                                                             ccTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTc                                                                                                                >  2:308907/1‑250 (MQ=255)                                                                                                                  gATAGTAAATCCTGCATAACTGT‑TCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGTTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGGCCTACCGGAACAGCCACCAGACTGTTGTAGCCTTTCTCTTTCAGACCa                                                                                                          <  2:333184/251‑1 (MQ=255)                                                                                                                   atatTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGAcaaa                                                                                                         <  2:279401/251‑4 (MQ=255)                                                                                                                                   aaC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGAt                                                                                         >  1:206822/1‑251 (MQ=255)                                                                                                                                             tttATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTc                                                                                >  1:165707/1‑251 (MQ=255)                                                                                                                                                 taCAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGCCTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCaaa                                                                             >  2:169221/1‑250 (MQ=255)                                                                                                                                                       agcagcCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCg                                                                      >  2:85283/1‑251 (MQ=255)                                                                                                                                                               atgatgACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGt                                                                <  1:378749/249‑1 (MQ=255)                                                                                                                                                                               cAGGGTTATACAAATCAGGAGAATCTTAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAATGTAATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCTGAACATCCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCaga                                              <  1:192327/251‑1 (MQ=255)                                                                                                                                                                                     tATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAgccgcc                                        >  2:52916/1‑251 (MQ=255)                                                                                                                                                                                                   gaATCTGAGAGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGa                           >  2:196077/1‑250 (MQ=255)                                                                                                                                                                                                         gagaGGAAATAGCCGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGtt                       <  1:93089/248‑1 (MQ=255)                                                                                                                                                                                                                      cGGGCAGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATa        >  2:14898/1‑250 (MQ=255)                                                                                                                                                                                                                           aGATGCCCAGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTg  >  2:113028/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   CCACCCAATGCAAAATAATATGCACCACGACGGCGGTCAGAAAAATAATACCGAAAATCATCACCATCGAGGTGGTGTGATCGATCTCAATACCCGCTAAATCTTCAACCTGGGATATTAAATCCTGCATAAC‑GTCTCCTTTATACAACAGCAGCCTATGATGACGGCTGAAACAGGGTTATGCAAATCAGGAGAATCTGAGAGGAAATAGCCGGGCAGATGCCCGGCAAGAGAGAATTACACTTCGGTTAAGGTGATGTTTTGCGGCAGACGAGATTTCGGCAGCGTAGCGTTAAAATCTTCAACGCTGTGATGACCTACCGGAACAACCACCAGACTGGTGTAGCCTTTCTCTTTCAGACCAAATTCTGCATCGAGGATGGCGGCGTCAAAACCTTCGATGGGTACCGCGTCCAGACCCAGAGCCGCCACGCCGAGCAGGAAGTTACCGACGTTGAGATAAACCTG  >  NC_000913/604534‑605001

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