Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,354,487	C→T	intergenic (‑437/‑238)	yhcC ← / → gltB	radical SAM family oxidoreductase YhcC/glutamate synthase subunit GltB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,354,487	0	C	T	97.3%	116.8 / NA	37	intergenic (‑437/‑238)	yhcC/gltB	radical SAM family oxidoreductase YhcC/glutamate synthase subunit GltB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (16/20);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAG  >  NC_000913/3354366‑3354600                                                                                                                          |                                                                                                                  tttCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGa                                                                                                        >  1:11567/1‑133 (MQ=255) tttCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGa                                                                                                        >  2:426816/1‑133 (MQ=255)  ttCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGAt                                                                                                       >  2:158765/1‑133 (MQ=255)         aCGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATCAGCTTGCCATTTGACCCGTATTAGCGTTCCCGATAAGTTgg                                                                                                >  2:277992/1‑133 (MQ=255)                      tAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGaa                                                                                   >  1:614382/1‑133 (MQ=255)                      tAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTAACCTGTATTAGCTTTCACGATAAGTTGGAATTCCGCTGGaa                                                                                   >  2:212112/1‑133 (MQ=255)                            tAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCGCGATAAGTTGGGAAACCGCTGGAAGCTTTc                                                                             >  1:186892/1‑133 (MQ=255)                            tAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTc                                                                             >  2:480066/1‑133 (MQ=255)                                           ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                              >  2:92748/1‑133 (MQ=255)                                           ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                              <  2:614382/133‑1 (MQ=255)                                           ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                              <  2:94970/133‑1 (MQ=255)                                           ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                              <  1:526096/133‑1 (MQ=255)                                                aTTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCtcat                                                         <  1:28440/133‑1 (MQ=255)                                                aTTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCtcat                                                         <  1:480066/133‑1 (MQ=255)                                                 tttCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCtcatc                                                        >  1:56914/1‑133 (MQ=255)                                                     aGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCatat                                                    >  1:63046/1‑133 (MQ=255)                                                      gTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                              <  1:97728/122‑1 (MQ=255)                                                      gTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                              >  2:97728/1‑122 (MQ=255)                                                        cATTTAATGAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTa                                                 <  1:277992/133‑1 (MQ=255)                                                        cATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTa                                                 <  2:376568/133‑1 (MQ=255)                                                          tttAATAATGAATTTTGCTCTAGAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATg                                               <  1:212112/133‑1 (MQ=255)                                                          tttAATAAAGAATTTGGCGCTAAATCACCTTTCTGTACCAAAAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATg                                               <  2:186892/133‑1 (MQ=255)                                                           ttAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGc                                              >  1:8419/1‑133 (MQ=255)                                                                 aaGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAAt                                        >  2:594845/1‑133 (MQ=255)                                                                     aTTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTgag                                    <  2:56914/133‑1 (MQ=255)                                                                                 aGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                        <  1:426816/133‑1 (MQ=255)                                                                                 aGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                        <  1:417050/133‑1 (MQ=255)                                                                                 aGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                        <  2:8419/133‑1 (MQ=255)                                                                                 aGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                        <  1:321461/133‑1 (MQ=255)                                                                                 aGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                        <  2:63046/133‑1 (MQ=255)                                                                                 aGCACATTTCCGTACCAATAAGCTTGCCATTTTACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                        <  2:177148/133‑1 (MQ=255)                                                                                      aTTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTa                   <  2:452756/133‑1 (MQ=255)                                                                                             tACCAATAAGCTTGCGATTTGACCTGTATTTGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGATTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAAcag            <  1:311865/133‑3 (MQ=255)                                                                                               ccAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGAt          <  1:455235/133‑1 (MQ=255)                                                                                                  aTAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCg       <  1:92748/133‑1 (MQ=255)                                                                                                  aTAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCg       >  1:152123/1‑133 (MQ=255)                                                                                                       cTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGCAATTGAGACCCCCTCTTCACCGTATTAACCGGTGCGaaagg  >  2:296977/1‑131 (MQ=255)                                                                                                                          |                                                                                                                  TTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAG  >  NC_000913/3354366‑3354600

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