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.

GCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGA  >  minE/1270323‑1270426
                                                               |                                          
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCt                                             >  1:115736/1‑62 (MQ=255)
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCt                                             >  1:121417/1‑62 (MQ=255)
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCt                                             >  1:221121/1‑62 (MQ=255)
gCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCt                                             >  1:586972/1‑62 (MQ=255)
                    ctttCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:168713/59‑1 (MQ=255)
                    cTATCAACAATTGTCGGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1123339/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATTGGCCAAACGAt                         <  1:1000987/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:213373/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:221809/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:544628/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:263134/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:275655/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:276004/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:281159/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:708829/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:315786/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:194640/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:326240/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:36094/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:366731/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:386029/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:664954/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:427744/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:437950/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:446596/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:647481/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:607142/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:566712/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:49500/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:501681/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1021059/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:879721/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1087068/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:109084/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:854898/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:109628/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:810419/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1053403/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1112745/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:112073/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1131718/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1135178/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:794403/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:787859/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1192860/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1020673/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:769925/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:905840/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:138572/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:170213/62‑1 (MQ=255)
                    cTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:764659/62‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:983482/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:547108/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:879791/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:684798/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:80246/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:697792/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:712559/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:715955/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:744355/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:766002/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1008893/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:490829/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1190574/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1135616/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:179532/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1104415/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1056971/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:309582/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1065463/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:1095788/61‑1 (MQ=255)
                     tATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:40386/61‑1 (MQ=255)
                     tATCAACAATTGTCAGACAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:837092/61‑1 (MQ=255)
                     tATCAACAATTGTCACGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:121188/61‑1 (MQ=255)
                      aTCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:273976/60‑1 (MQ=255)
                       tCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:608257/59‑1 (MQ=255)
                                gTCAGGCAAACGCCCAACCACACTTGGCCTGGGGCGATGGGCCAAACGAt                         <  1:951989/50‑1 (MQ=255)
                                     ccAAACGCCCAACCACACTTGGCCGGGGGCGATGGGCCAAACGAt                         <  1:298294/44‑1 (MQ=255)
                                        aaCGCCCAACCACACTTGGCCGGGGGCGATGGGCCAAACGAt                         <  1:496333/42‑1 (MQ=255)
                                          cGCCCAACCATACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:636999/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:241620/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:819919/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:1107999/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:1012164/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:906473/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:170475/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:765493/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:268539/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:70652/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:401922/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:681190/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:661590/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:466292/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:631065/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:477627/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:580042/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:556569/62‑1 (MQ=255)
                                          cGCCCAACCACACTTGGCCGGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:808970/62‑1 (MQ=255)
                                           gCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGa  <  1:489765/61‑1 (MQ=255)
                                                       ttGGCCTGGGCGATGGGCCAAACGGTGCGCCCTTACTGGAGGTAATGGa  <  1:944690/49‑1 (MQ=255)
                                                               |                                          
GCCAACTGGATTATCGCGACCTATCAACAATTGTCAGGCAAACGCCCAACCACACTTGGCCTGGGCGATGGGCCAAACGATGCGCCCTTACTGGAGGTAATGGA  >  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: