Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,270,387 +G coding (638/816 nt) yedP → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,270,3841.G76.3% 323.6 / 93.2 114coding (635/816 nt)yedPconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/27);  new base G (0/87);  total (0/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

GCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGA  >  minE/1270323‑1270426
                                                               |                                          
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑t                                             >  1:115736/1‑62 (MQ=255)
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑t                                             >  1:121417/1‑62 (MQ=255)
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑t                                             >  1:221121/1‑62 (MQ=255)
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑t                                             >  1:586972/1‑62 (MQ=255)
                    ctttCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:168713/59‑1 (MQ=255)
                    cTATCAACAATTGTCGGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1123339/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATTGGCCAAACGAt                         <  1:1000987/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:213373/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:221809/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:544628/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:263134/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:275655/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:276004/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:281159/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:708829/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:315786/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:194640/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:326240/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:36094/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:366731/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:386029/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:664954/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:427744/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:437950/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:446596/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:647481/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:607142/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:566712/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:49500/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:501681/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1021059/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:879721/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1087068/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:109084/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:854898/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:109628/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:810419/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1053403/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1112745/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:112073/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1131718/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1135178/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:794403/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:787859/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1192860/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1020673/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:769925/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:905840/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:138572/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:170213/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:764659/62‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:983482/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:547108/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:879791/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:684798/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:80246/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:697792/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:712559/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:715955/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:744355/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:766002/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1008893/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:490829/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1190574/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1135616/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:179532/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1104415/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1056971/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:309582/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1065463/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:1095788/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:40386/61‑1 (MQ=255)
                     tATCAACAATTGTCAGACAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:837092/61‑1 (MQ=255)
                     tATCAACAATTGTCACGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:121188/61‑1 (MQ=255)
                      aTCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:273976/60‑1 (MQ=255)
                       tCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:608257/59‑1 (MQ=255)
                                gTCAGGCAAACGCCCAACCACACTTGGCC‑TGGGGCGATGGGCCAAACGAt                         <  1:951989/50‑1 (MQ=255)
                                     ccAAACGCCCAACCACACTTGGCCGG‑GGGCGATGGGCCAAACGAt                         <  1:298294/44‑1 (MQ=255)
                                        aaCGCCCAACCACACTTGGCCGG‑GGGCGATGGGCCAAACGAt                         <  1:496333/42‑1 (MQ=255)
                                          cGCCCAACCATACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:636999/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:241620/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:819919/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:1107999/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:1012164/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:906473/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:170475/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:765493/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:268539/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:70652/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:401922/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:681190/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:661590/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:466292/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:631065/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:477627/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:580042/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:556569/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCC‑G‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:808970/62‑1 (MQ=255)
                                           gCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:489765/61‑1 (MQ=255)
                                                       ttGGCC‑T‑GGGCGATGGGCCAAACGGTGCGCCCTTACTGGAGGTAATGGa  <  1:944690/49‑1 (MQ=255)
                                                               |                                          
GCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCC‑T‑GGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGA  >  minE/1270323‑1270426

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: ‑