Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,200,454	(A)8→7	coding (176/3633 nt)	yehI →	DUF4132 domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,200,447	0	A	.	100.0%	119.3 / NA	29	coding (169/3633 nt)	yehI	DUF4132 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (13/16);  total (13/16)

CCCAACTGGAACTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCT  >  NC_000913/2200313‑2200573                                                                                                                                       |                                                                                                                               cccAACTGGAACTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaa                                                                                                                            >  1:811857/1‑139 (MQ=255)     aCTGGAACTGAAATATAAAAAAGGCAAAACGCAGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaaca                                                                                                                        >  1:927891/1‑137 (MQ=255)        ggAACTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGT‑CGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGt                                                                                                                   <  1:298162/139‑1 (MQ=255)               aaaTATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAATATTGTCTGGCg                                                                                                             >  2:866619/1‑139 (MQ=255)                    taAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATg                                                                                                        >  2:511654/1‑139 (MQ=255)                         aaGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAAt                                                                                                   >  1:928110/1‑139 (MQ=255)                            gCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAAtttt                                                                                                >  2:222577/1‑139 (MQ=255)                              aaaaCGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTc                                                                                              <  2:69342/139‑1 (MQ=255)                                aaCGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCAc                                                                                            <  1:222577/139‑1 (MQ=255)                                aaCGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCAc                                                                                            >  1:584729/1‑139 (MQ=255)                                                      tGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaaca                                                                                                                        >  2:748525/1‑88 (MQ=255)                                                      tGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGaaaaaaaca                                                                                                                        <  1:748525/90‑3 (MQ=255)                                                      tGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCAc                                                                      >  2:195883/1‑139 (MQ=255)                                                       ggCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACg                                                                     <  1:789851/139‑1 (MQ=255)                                                         ccGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                                   <  2:927891/139‑1 (MQ=255)                                                            ggCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGc                                                                >  2:856889/1‑139 (MQ=255)                                                                      tgtCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                                   >  1:473725/1‑126 (MQ=255)                                                                      tgtCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCa                                                                   <  2:473725/126‑1 (MQ=255)                                                                      tgtCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTg                                                      >  2:598662/1‑139 (MQ=255)                                                                           gTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATTAACTTGATCgg                                                 <  2:584729/139‑1 (MQ=255)                                                                            ttATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAACCATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGc                                                <  1:951632/139‑1 (MQ=255)                                                                             tATCACTGGAAGTCCCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAgg                                                                 <  2:678115/121‑1 (MQ=255)                                                                             tATCACTGGAAGTCCCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAgg                                                                 >  1:678115/1‑121 (MQ=255)                                                                                 aCTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAaca                                           >  2:383530/1‑139 (MQ=255)                                                                                   tGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAAcaca                                         <  1:856889/139‑1 (MQ=255)                                                                                     gAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACaaa                                       >  2:938538/1‑139 (MQ=255)                                                                                               ccAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATcc                             <  1:598662/139‑1 (MQ=255)                                                                                                                   tGCTACAAAAAGCTGAGAAG‑AAAAAAATATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg         <  1:866619/139‑1 (MQ=255)                                                                                                                   tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg         <  1:195883/139‑1 (MQ=255)                                                                                                                   tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg         <  1:643686/139‑1 (MQ=255)                                                                                                                   tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg         <  1:511654/139‑1 (MQ=255)                                                                                                                   tGCTACAAAAAGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGgcggcg         <  2:928110/139‑1 (MQ=255)                                                                                                                          aaaaGCTGAGAAG‑AAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCt  >  2:38524/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                               CCCAACTGGAACTGAAATATAAAAAAGGCAAAACGCCGCTCAGTCATCGTCGCTGGCCGGGCGAACCAGTGTCCGTTATCACTGGAAGTCTCATCCAGACATTGGGTGATGAATTGCTACAAAAAGCTGAGAAGAAAAAAAACATTGTCTGGCGTTATGAGAATTTTTCACTGGAGTGGCAGTCCGCCATCACGCAGGCCATCAACTTGATCGGCGAACACAAACCCTCAATCCCGGCCCGGACAATGGCGGCGCTAGCCT  >  NC_000913/2200313‑2200573

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