Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 913,779 A→G L153P (CTG→CCG)  hcp ← hybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb913,7790AG100.0% 302.8 / NA 81L153P (CTG→CCG) hcphybrid‑cluster [4Fe‑2S‑2O] protein in anaerobic terminal reductases
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (0/81);  total (0/81)

ATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGA  >  W3110S.gb/913714‑913781
                                                                 |  
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCGCAGCGGa  <  1:1325414/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:409296/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:531463/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:495748/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:481295/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:47952/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:459414/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:434853/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:433015/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:422638/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:531680/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:388189/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:367405/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:356089/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:312233/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:301224/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:294078/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:248535/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:236564/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1000485/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:794650/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:975141/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:970933/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:928078/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:905255/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:883889/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:87523/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:856144/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:852890/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:843506/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1596718/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:764963/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:757848/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:679556/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:656057/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:604690/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:562091/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:560293/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:552790/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:552317/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1113472/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1595110/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1269684/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1262260/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1257001/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1254075/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1211262/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1205278/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1188862/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1149196/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1270303/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:111329/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1082684/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:106673/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:106313/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1052206/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1039703/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1029645/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:102952/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1408206/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1501945/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1500167/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1540268/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1497983/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1472432/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1541878/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1559126/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1430387/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1393266/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1393017/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1320792/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1307437/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1306077/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1276287/68‑1 (MQ=255)
aTTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1300811/68‑1 (MQ=255)
 ttGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:22776/67‑1 (MQ=255)
 ttGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:108825/67‑1 (MQ=255)
       gagaACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1167773/61‑1 (MQ=255)
                        ccATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:541445/44‑1 (MQ=255)
                          atatAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:1462611/42‑1 (MQ=255)
                              aGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCGGa  <  1:243312/38‑1 (MQ=255)
                                                                 |  
ATTGACCGAGAACATGCGCGTGTTCCATATAGGCCGCCGCACCTTTCAGGCCATACAGGCACAGCAGA  >  W3110S.gb/913714‑913781

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