Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,968,005	(A)7→8	coding (90/1104 nt)	rfe →	UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,998	1	.	A	100.0%	89.6 / NA	25	A28E (GCA→GAA)	rfe	UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (13/12);  total (13/12)

AATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGC‑AAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGG  >  NC_000913/3967919‑3968070                                                                                 |                                                                         aaTTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGc                                                                           >  2:725610/1‑80 (MQ=255) aaTTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGc                                                                           <  2:725562/80‑1 (MQ=255)     tACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGC‑aaaa                                                                      <  1:725564/80‑1 (MQ=255)     tACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGC‑aaaa                                                                      <  2:725564/80‑1 (MQ=255)      aCTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGC‑aaaaa                                                                     >  1:725620/1‑80 (MQ=255)               gAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGtt                                                            >  1:725616/1‑81 (MQ=255)               gAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGtt                                                            >  1:725604/1‑81 (MQ=255)                 gTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGtt                                                            >  2:725640/1‑79 (MQ=255)                 gTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTa                                                          <  1:725566/81‑1 (MQ=255)                    cTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGt                                                        <  1:725567/80‑1 (MQ=255)                           catcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAAc                                                >  1:725618/1‑81 (MQ=255)                            atcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGAt                                                    <  1:725568/76‑1 (MQ=255)                            atcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGAt                                                    <  2:725569/76‑1 (MQ=255)                            atcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATa                                                   <  2:725568/77‑1 (MQ=255)                                   tttttttATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACtt                                         <  1:725571/80‑1 (MQ=255)                                      ttttATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCg                                      >  2:725612/1‑80 (MQ=255)                                              aCGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTc                              <  2:725571/80‑1 (MQ=255)                                               cGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCAc                            >  1:725626/1‑81 (MQ=255)                                                gACACTGTTTCTGTCTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCAc                            >  2:725618/1‑80 (MQ=255)                                                 acacTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCAc                            >  1:725629/1‑79 (MQ=255)                                                   acTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAgg                        >  1:725630/1‑81 (MQ=255)                                                   acTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAgg                        >  1:725627/1‑81 (MQ=255)                                                            gtttTTTGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATAc               >  1:725623/1‑81 (MQ=255)                                                             ttttttGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATAcc              <  1:725573/81‑1 (MQ=255)                                                               ttttGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACct             <  1:725574/80‑1 (MQ=255)                                                                  tGCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCtctc          <  1:725575/80‑1 (MQ=255)                                                                   gCCCGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGt        >  1:725621/1‑81 (MQ=255)                                                                    cccGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGtt       <  1:725576/81‑1 (MQ=255)                                                                     ccGTAAGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTg      <  1:725577/81‑1 (MQ=255)                                                                          aGGTGGCAAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTggggg  >  2:725638/1‑80 (MQ=255)                                                                                 |                                                                         AATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGC‑AAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGG  >  NC_000913/3967919‑3968070

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