Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 458,686 G→A R192H (CGT→CAT)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb458,6860GA86.3% 126.3 / 12.4 51R192H (CGT→CAT) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (0/7);  new base A (0/44);  total (0/51)
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.75e-01

ACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCG  >  W3110S.gb/458639‑458717
                                               |                               
agtGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2914875/60‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:1835050/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2326629/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2154114/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2567645/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2596487/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:1969280/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:1960071/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2625435/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2921932/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:786960/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:802381/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:861539/62‑1 (MQ=255)
aCTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2858038/62‑1 (MQ=255)
 cTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:2120723/61‑1 (MQ=255)
 cTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:325467/61‑1 (MQ=255)
 cTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:866888/61‑1 (MQ=255)
           aaCAGTCTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGa                   <  1:3005020/51‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTTGAATATCTGATGGCAATGATGGAATCg  <  1:745212/62‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATGGAATCg  <  1:2148901/62‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATGGAATCg  <  1:2241991/62‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATGGAATCg  <  1:219858/62‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATGGAATCg  <  1:448623/62‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATGGAATCg  <  1:70592/62‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATGGAATCg  <  1:1332114/62‑1 (MQ=255)
                 cTGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATGGAATCg  <  1:1207719/62‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:359839/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2575144/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2995067/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:575279/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:630173/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2830620/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2139113/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:1146632/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:145653/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:1628052/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:1691767/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:1727766/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2815795/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2392450/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2452309/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:1022113/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:2694898/56‑1 (MQ=255)
                  tGTTCTGGAGATGTCCGACGTTAACGAACATCTGGAATATCTGATGGCAATGATgg       <  1:273316/56‑1 (MQ=255)
                         gagaTGTCCGACGTTAACGAACGTCTGGAATTTCt                     <  1:1620166/35‑1 (MQ=255)
                         gagaTGTCCGACGTTAACGAACGTCTGGAATATCt                     <  1:247444/35‑1 (MQ=255)
                         gagaTGTCCGACGTTAACGAACGTCTGGAATATCt                     <  1:1985293/35‑1 (MQ=255)
                         gagaTGTCCGACGTTAACGAACGTCTGGAATATCt                     <  1:1769512/35‑1 (MQ=255)
                         gagaTGTCCGACGTTAACGAACGTCTGGAATATCt                     <  1:695516/35‑1 (MQ=255)
                         gagaTGTCCGACGTTAACGAACGTCTGGAATATCt                     <  1:2929589/35‑1 (MQ=255)
                         gagaTGTCCGACGTTAACGAACGTCTGGAATATCt                     <  1:2664692/35‑1 (MQ=255)
                                               |                               
ACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCG  >  W3110S.gb/458639‑458717

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