Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	442,961	G→A	E471E (GAG→GAA)	thiI →	tRNA uridine 4‑sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	442,961	0	G	A	100.0%	74.8 / NA	23	E471E (GAG→GAA)	thiI	tRNA uridine 4‑sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/13);  total (10/13)

ATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGA  >  NC_000913/442833‑443095                                                                                                                                 |                                                                                                                                       aTTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTaa                                                                                                                             <  2:1501167/140‑1 (MQ=255)   tGATGTGGTTTCTCTGCCGTTCTATAAATTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAAct                                                                                                                           >  1:1275396/1‑139 (MQ=255)          gTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGaa                                                                                                                    >  2:1599532/1‑140 (MQ=255)                      ttCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGccc                                                                                                        <  2:1504187/140‑1 (MQ=255)                         tataAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTa                                                                                                     >  2:1597030/1‑140 (MQ=255)                              aCTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTtg                                                                                                <  2:1429558/140‑1 (MQ=255)                                  aGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGt                                                                                            >  1:1496521/1‑140 (MQ=255)                                  aGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGt                                                                                            >  2:1541647/1‑140 (MQ=255)                                  aGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGt                                                                                            >  1:1245221/1‑140 (MQ=255)                                    cACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAAC‑‑TGTGAAGGTATATCGCCCGTAATTTGTGGttttt                                                                                        <  2:1275396/140‑1 (MQ=255)                                          aTTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGt                                                                                    <  2:1525137/140‑1 (MQ=255)                                                                 aaaCCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGtt                                                             >  1:1659137/1‑140 (MQ=255)                                                                   aCCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTg                                                           <  1:1597030/140‑1 (MQ=255)                                                                     cTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGcc                                                         <  2:1245221/140‑1 (MQ=255)                                                                       ggCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGc                                                       >  2:1372911/1‑140 (MQ=255)                                                                             tgtgGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGa                                                 <  1:1599532/140‑1 (MQ=255)                                                                               tggtgTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGaca                                               <  2:1591371/140‑1 (MQ=255)                                                                                                  tgaGCCGTCTGCAGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATaa                            >  1:1407701/1‑140 (MQ=255)                                                                                                              aGGCGCTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTgccgcc                <  1:1400798/140‑1 (MQ=255)                                                                                                                gcgcTCTATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCAt              >  2:1745065/1‑140 (MQ=255)                                                                                                                     ctATCTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCa         <  2:1659137/140‑1 (MQ=255)                                                                                                                         cTGCGCGAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTtg     <  1:1541647/140‑1 (MQ=255)                                                                                                                            cgcgAACAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGa  <  2:1496521/140‑1 (MQ=255)                                                                                                                                 |                                                                                                                                       ATTGATGTGGTTTCTCTGCCGTTCTATAAACTGAGCACCAAATTTGGCGATCTCGACCAGAACAAAACCTGGCTGCTGTGGTGTGAGCGCGGGGTGATGAGCCGTCTGCAGGCGCTCTATCTGCGCGAGCAGGGCTTTAACAATGTGAAGGTATATCGCCCGTAATTTGTGGTTTTTACGTCGCATCTGGTCAGATGCGACGTTTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGA  >  NC_000913/442833‑443095

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