Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,299,535 T→C G40G (GGA→GGG)  ccmC ← heme exporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,299,5350TC98.7% 578.6 / ‑2.5 152G40G (GGA→GGG) ccmCheme exporter subunit
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (150/0);  total (152/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 152 total aligned reads displayed.

TGCAGGTAGATAATGCGGTAGCTATTTCCCTGCTGATAATCAGCCGGAGCAAAGCCGAATCCCCAGATCCAGCCGACGGTAAGCACGACCACACT  >  W3110S.gb/2299476‑2299570
                                                           |                                   
tGCAGGTAGATAATGCGGTAGCTATTTCCCTGCTGATAATCAGCCGGAGCAAAGCCGAATc                                    >  1:864334/1‑61 (MQ=255)
tGCAGGTAGATAATGCGGTAGCTATTTCCCTGCTGATAATCAGCCGGAGCAAAGCCGAATc                                    >  1:83964/1‑61 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGCACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:108550/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGACCCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:227192/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATcc                           >  1:179836/1‑40 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCa                          >  1:1678224/1‑41 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCa                          >  1:123424/1‑41 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACg                   >  1:400787/1‑48 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTa                >  1:1245251/1‑51 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTACGCACGACCACACt  >  1:596389/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACg          >  1:1206780/1‑57 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGCCCACACt  >  1:749101/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACcaca    >  1:6972/1‑63 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCacac   >  1:526554/1‑64 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCacac   >  1:1885752/1‑64 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:244236/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:630038/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:506000/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:504232/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:226074/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:477534/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:473128/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:42173/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:357427/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:271720/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:383713/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:381011/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:277379/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:371119/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:346062/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:352633/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:766295/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:957313/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:945935/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:941648/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:895726/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:844559/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:826318/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:79659/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:787735/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:781984/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:775989/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:534265/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:745822/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:73602/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:716383/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:680949/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:651646/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:640761/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1870977/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:6161/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:590971/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:559399/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1233039/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1465565/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1371496/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1360860/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1353272/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1338085/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1336822/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1325837/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1319267/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1275863/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:127000/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1256828/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1469783/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1213509/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1187576/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1173738/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1115222/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1108488/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1078325/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1063388/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:106012/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1026833/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1025557/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1007781/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1715151/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1877254/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:100100/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1819956/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1818519/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1804680/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1765303/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1751363/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1746558/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1731940/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1728581/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:204764/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1704171/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1693676/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1692309/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1679807/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1625892/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1610942/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1563412/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1516737/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1502041/1‑65 (MQ=255)
                              tgctgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGACCACACt  >  1:1480178/1‑65 (MQ=255)
                              tgatgATAATCAGCCGGAGCAAAGCCGAACCCCCAGATCCAGCCGACGGTAAGCACGCCCACACt  >  1:243120/4‑65 (MQ=255)
                                                           |                                   
TGCAGGTAGATAATGCGGTAGCTATTTCCCTGCTGATAATCAGCCGGAGCAAAGCCGAATCCCCAGATCCAGCCGACGGTAAGCACGACCACACT  >  W3110S.gb/2299476‑2299570

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