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	100.0%	120.2 / NA	36	intergenic (‑437/‑238)	yhcC/gltB	radical SAM family oxidoreductase YhcC/glutamate synthase subunit GltB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/20);  total (16/20)

ACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGG  >  NC_000913/3354362‑3354613                                                                                                                              |                                                                                                                               aCTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTc                                                                                                                         >  2:88010/1‑133 (MQ=255)     tttCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGa                                                                                                                     >  2:407877/1‑133 (MQ=255)     tttCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGa                                                                                                                     <  2:274234/133‑1 (MQ=255)       tctcTAACAACTTCCGTTGTAGGCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATa                                                                                                                   <  1:253363/133‑1 (MQ=255)               gACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTCACCTGTAGTAGCTTCCCCGATAAGTTGGaa                                                                                                           >  1:174311/1‑133 (MQ=255)                aCTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGaaa                                                                                                          >  2:348863/1‑133 (MQ=255)                         ttAGCCTTAAGGATAAAATCCATTTTAAATTCCGTCATTTAATAAAGAATTTTGTGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGa                                                                                                 <  2:213456/133‑1 (MQ=255)                            gCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCCTTCCCGAGAAGGTGGAAATCCGCTGGAAGc                                                                                              >  1:160932/1‑133 (MQ=255)                              cTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCtt                                                                                            <  1:88010/133‑1 (MQ=255)                                       aaaaTCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGa                                                                                   >  1:417959/1‑133 (MQ=255)                                         aaTCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGagc                                                                                 >  2:498721/1‑133 (MQ=255)                                              aTTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGa                                                                                   <  1:274902/126‑1 (MQ=255)                                              aTTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGa                                                                                   >  2:274902/1‑126 (MQ=255)                                               ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                                           <  1:405962/133‑1 (MQ=255)                                               ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                                           <  1:348863/133‑1 (MQ=255)                                               ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                                           >  1:58854/1‑133 (MQ=255)                                               ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                                           <  1:498293/133‑1 (MQ=255)                                               ttttAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                                           <  1:498721/133‑1 (MQ=255)                                               ttttAATTACAGTCATTTCATATAGAATTTCGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTg                                                                           <  1:145228/133‑1 (MQ=255)                                                            cATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTa                                                              <  2:58854/133‑1 (MQ=255)                                                            cATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTa                                                              <  1:462875/133‑1 (MQ=255)                                                              tttAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATg                                                            <  1:114926/133‑1 (MQ=255)                                                              tttAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATg                                                            <  1:423657/133‑1 (MQ=255)                                                                 aataaAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAg                                                         <  2:499724/133‑1 (MQ=255)                                                                         aTTTTGCGCTTAAGCTCATTTCTGAACCCATATGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTgag                                                 <  2:160932/133‑1 (MQ=255)                                                                         aTTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTgag                                                 >  1:441379/1‑133 (MQ=255)                                                                                gcTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCct                                          <  2:402039/133‑1 (MQ=255)                                                                                     aGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                                     <  1:407877/133‑1 (MQ=255)                                                                                     aGCACATTTCTGTACAAATAAGCTTGCCCTTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCa                                     <  2:174311/133‑1 (MQ=255)                                                                                                 tACCAATAAGCTTGCCATTTGACCTGGATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCg                         >  2:473375/1‑133 (MQ=255)                                                                                                    cAATAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATg                      >  1:404976/1‑133 (MQ=255)                                                                                                      aTAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCg                    >  2:394955/1‑133 (MQ=255)                                                                                                      aTAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCg                    >  2:464827/1‑133 (MQ=255)                                                                                                      aTAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCg                    >  2:258037/1‑133 (MQ=255)                                                                                                      aTAAGCTTGCCATTTGACCTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCg                    >  2:12880/1‑133 (MQ=255)                                                                                                                        cTGTATTAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAggggg  <  1:473375/133‑1 (MQ=255)                                                                                                                              |                                                                                                                               ACTGTTTCTCTAACGACTTCCCTTTTAGCCTTAAAGATAAAATCCATTTTAATTTCAGTCATTTAATAAAGAATTTTGCGCTAAAGCACATTTCTGTACCAATAAGCTTGCCATTTGACCTGTATCAGCTTTCCCGATAAGTTGGAAATCCGCTGGAAGCTTTCTGGATGAGCAGCCTGCTCATCATATTTATGCAGTAATTGAGATCCCCTCTTCACCGTATTAACCGATGCGAAAAGGACAACAAGGGGG  >  NC_000913/3354362‑3354613

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: ‑