Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,113,386	G→A	A32T (GCT→ACT)	uspD →	universal stress protein D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,386	0	G	A	100.0%	113.1 / NA	33	A32T (GCT→ACT)	uspD	universal stress protein D
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/22);  total (11/22)

TTGGCGTGGCAATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAA  >  NC_000913/4113309‑4113460                                                                              |                                                                           ttGGCGTGGCAATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTc                                                                         <  1:814145/81‑1 (MQ=255)  tGGCGTGGCAATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCa                                                                        >  2:814190/1‑81 (MQ=255)   ggCGTGGCAATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCAc                                                                       >  2:814188/1‑81 (MQ=255)       tGGCAATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGa                                                                   <  2:814146/81‑1 (MQ=255)          cAATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGt                                                                <  2:814147/81‑1 (MQ=255)           aaTTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGAcct                                                                >  2:814183/1‑78 (MQ=255)              ttCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAAt                                                            >  2:814185/1‑81 (MQ=255)               tCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAAtt                                                           >  1:814186/1‑81 (MQ=255)               tCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAAtt                                                           >  1:814196/1‑81 (MQ=255)               tCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAAt                                                            <  2:814148/80‑1 (MQ=255)                 cGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTc                                                          <  1:814149/80‑1 (MQ=255)                 cGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTc                                                          <  2:814150/80‑1 (MQ=255)                  gggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCat                                                        <  1:814151/81‑1 (MQ=255)                  gggAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCa                                                         >  2:814198/1‑80 (MQ=255)                     aaTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATAtt                                                     <  2:814152/81‑1 (MQ=255)                              gaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGc                                            >  2:814202/1‑81 (MQ=255)                              gaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGc                                            <  2:814153/81‑1 (MQ=255)                                tGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACTTGACGTTAATTCATATTGATGATGGCtt                                          <  1:814155/81‑1 (MQ=255)                                tGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCtt                                          <  1:814154/81‑1 (MQ=255)                                tGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCtt                                          <  2:814156/81‑1 (MQ=255)                                  ccTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTaa                                        <  1:814157/81‑1 (MQ=255)                                       cTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGa                                    <  1:814159/80‑1 (MQ=255)                                       cTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAg                                   <  2:814158/81‑1 (MQ=255)                                          gTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGt                                  <  2:814160/79‑1 (MQ=255)                                                 aaGCCCTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCg                          <  1:814161/80‑1 (MQ=255)                                                     ccTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTAt                     >  2:814193/1‑81 (MQ=255)                                                      cTGGAGCTCGCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATc                    <  2:814162/81‑1 (MQ=255)                                                               gCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTcc             <  2:814164/79‑1 (MQ=255)                                                               gCCAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCt           <  2:814163/81‑1 (MQ=255)                                                                ccAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTg          <  2:814165/81‑1 (MQ=255)                                                                 cAGACATAATGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTg          >  2:814192/1‑80 (MQ=255)                                                                        aaTGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACag    >  2:814199/1‑79 (MQ=255)                                                                         aTGACACTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACagaa  <  2:814166/80‑1 (MQ=255)                                                                              |                                                                           TTGGCGTGGCAATTTCCGGGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAA  >  NC_000913/4113309‑4113460

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