Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,158,125 A→T F111Y (TTT→TAT)  astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,158,1250AT77.8% 191.8 / 50.6 90F111Y (TTT→TAT) astAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (20/0);  new base T (0/70);  total (20/70)
Fisher's exact test for biased strand distribution p-value = 1.96e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 107 total aligned reads displayed.

TTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTG  >  minE/1158058‑1158171
                                                                   |                                              
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:477108/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:205758/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:645089/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:1229449/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:1266347/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:23283/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:1310078/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:511748/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:787061/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:261125/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:1415008/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:1494063/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:1086365/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:446206/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:444375/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGaaa                                               >  1:1541537/1‑69 (MQ=255)
tttGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCG‑‑GTGATCGTTACTGAGAAACa                                             >  1:414401/1‑69 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGc                                     >  1:630047/1‑48 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCTATGCGTTATAGACATTCAGc                >  1:570073/1‑69 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCg                               >  1:448696/1‑54 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTc                   >  1:1262807/1‑66 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAg                 >  1:1041872/1‑68 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:160182/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:1043286/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:895820/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:820316/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:459003/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:1423080/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:1130019/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:1311879/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:1176939/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:553957/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGctc              >  1:580915/1‑71 (MQ=255)
                               cacaGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGct               >  1:785753/1‑70 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGTTACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:412245/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGTGGCGTg  <  1:200941/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:329659/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:298278/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:324731/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:352659/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:93026/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:365947/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:427905/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:442299/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:482807/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:492280/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:507927/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:514181/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:619492/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:649351/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:705138/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:762280/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:764568/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:771289/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:794724/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:897352/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:116232/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1327838/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1300144/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:283333/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1266512/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1245192/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1244881/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1222137/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1220557/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1181933/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:127622/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1157445/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1146412/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1142893/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1133251/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:11220/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1096913/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1085008/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1056371/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1479807/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1535125/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1514063/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1513249/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:15458/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:193124/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:210372/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:230888/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:235518/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1460598/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1375569/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:249567/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1363637/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1348842/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1344436/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:276521/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGGCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:744343/71‑1 (MQ=255)
                                           ctgcCGGTGTGATCGTTACTGAGATACAGCGCCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1500930/71‑1 (MQ=255)
                                            tgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGGCATTCAGCTCTTTTGAGGCGTg  <  1:1330932/70‑1 (MQ=255)
                                            tgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:659233/70‑1 (MQ=255)
                                            tgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:653954/70‑1 (MQ=255)
                                            tgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:1247110/70‑1 (MQ=255)
                                            tgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:914515/70‑1 (MQ=255)
                                            tgcCGGTGTGATCGTTACTGAGATACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:52358/70‑1 (MQ=255)
                                                            aCTGAGATACAGCGTCGGCATTGCGTTATAGACATTCAGCTCTTTTGAGGCGTg  <  1:755330/54‑1 (MQ=255)
                                                                   |                                              
TTTGCGCCAGTCCGGGTCGAGAAACAGCGTGCACAGCTCGCTGCTGCCGGTGTGATCGTTACTGAGAAACAGCGTCGGCAATGCGTTATAGACATTCAGCTCTTTTGAGGCGTG  >  minE/1158058‑1158171

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: