Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,591,415 G→T K203N (AAG→AAT yfeS → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,591,4150GT100.0% 341.2 / NA 100K203N (AAG→AATyfeSconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (68/32);  total (68/32)

Only 100 of 122 total aligned reads displayed.

ATGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGAC  >  minE/1591347‑1591475
                                                                    |                                                            
atGATTACCTACGCCACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:907791/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:3774666/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:2752596/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:2749382/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:2725443/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:276735/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:2816468/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:2891705/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:3098983/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:3179085/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:3465400/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:2353428/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:2138549/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:3877909/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:3984424/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:1907784/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:823655/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:4200601/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:742613/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:1537553/1‑71 (MQ=255)
atGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAatat                                                            >  1:1696741/1‑71 (MQ=255)
 tGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATg                                                           >  1:3581896/1‑71 (MQ=255)
 tGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATg                                                           >  1:1482993/1‑71 (MQ=255)
 tGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATg                                                           >  1:2220678/1‑71 (MQ=255)
 tGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATg                                                           >  1:161199/1‑71 (MQ=255)
  gATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATgg                                                          >  1:1252196/1‑71 (MQ=255)
  gATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATgg                                                          >  1:3563883/1‑71 (MQ=255)
  gATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATgg                                                          >  1:578619/1‑71 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:1055945/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:81187/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:327032/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:3392967/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:2488033/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:3489299/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:2287801/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:2145440/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:482176/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:4001233/1‑70 (MQ=255)
       ccTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACt                                                      >  1:409631/1‑70 (MQ=255)
              aaCTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              <  1:2426304/71‑1 (MQ=255)
              aaCTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              <  1:4211122/71‑1 (MQ=255)
              aaCTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              <  1:1859867/71‑1 (MQ=255)
               aCTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              <  1:1060839/70‑1 (MQ=255)
                           cAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATgg                                                          >  1:2909030/1‑46 (MQ=255)
                           cAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCTCTAGCGCGtcat                                   >  1:3813111/1‑69 (MQ=255)
                           cAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAg                                           >  1:20990/1‑61 (MQ=255)
                           cAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtcat                                   >  1:619761/1‑69 (MQ=255)
                           cAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGtcat                                   >  1:1841767/1‑69 (MQ=255)
                           cAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGt                                      >  1:2794010/1‑66 (MQ=255)
                                 aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATgg                                                          >  1:4241544/1‑40 (MQ=255)
                                 aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGctc                                              >  1:1130245/1‑50 (MQ=255)
                                 aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              >  1:528118/1‑52 (MQ=255)
                                 aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              >  1:687075/1‑52 (MQ=255)
                                 aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              >  1:2085155/1‑52 (MQ=255)
                                 aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              >  1:793862/1‑52 (MQ=255)
                                 aaaGTCATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              >  1:1768558/1‑52 (MQ=255)
                                      cATGGGTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCAc                                              <  1:3777286/47‑1 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCa                                               >  1:2492143/1‑43 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa                   >  1:4001624/1‑71 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa                   >  1:1983305/1‑71 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa                   >  1:1364191/1‑71 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa                   >  1:2045663/1‑71 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa                   >  1:3512203/1‑71 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa                   >  1:2135304/1‑71 (MQ=255)
                                         gggTAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCaaaaa                   >  1:3736211/1‑71 (MQ=255)
                                            tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCtt                             >  1:208353/1‑58 (MQ=255)
                                            tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAtt                 >  1:2942540/1‑70 (MQ=255)
                                            tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAtt                 >  1:1726458/1‑70 (MQ=255)
                                            tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAtt                 >  1:3030669/1‑70 (MQ=255)
                                            tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAtt                 >  1:2507704/1‑70 (MQ=255)
                                            tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAtt                 >  1:2511700/1‑70 (MQ=255)
                                            tAAAATCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAt                  >  1:3992945/1‑69 (MQ=255)
                                              aaaTCTGCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttctt              >  1:564307/1‑71 (MQ=255)
                                              aaaTCTCCCATAAACTTAAAAATATGGGACTCAAt                                                  >  1:1859218/1‑35 (MQ=255)
                                              aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttctt              >  1:665094/1‑71 (MQ=255)
                                              aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttctt              >  1:4321608/1‑71 (MQ=255)
                                              aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttctt              >  1:3188483/1‑71 (MQ=255)
                                              aaaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttctt              >  1:3258005/1‑71 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:1420394/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:1531417/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:1355381/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:2659753/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:731049/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:722996/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:2543665/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:2701882/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:4381994/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:2559140/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:2967306/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:3597560/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:354294/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:2333654/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:3440572/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:3276584/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:3207490/71‑1 (MQ=255)
                                               aaTCTCCCATAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAAttcttc             <  1:2527653/71‑1 (MQ=255)
                                                        tAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTg    <  1:3751482/71‑1 (MQ=255)
                                                        tAAACTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTg    <  1:3366374/71‑1 (MQ=255)
                                                          aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGAc  >  1:43258/1‑71 (MQ=255)
                                                          aaCTTAAAAATATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGAc  >  1:4195200/1‑71 (MQ=255)
                                                                    |                                                            
ATGATTACCTACGCAACTGCATTCGGTCAGATTAAAGTCATGGGTAAAATCTCCCATAAACTTAAAAAGATGGGACTCAATGCACTAGCGCGTCATCAGCTTACCGCAAAAATTCTTCAATGGGGTGAC  >  minE/1591347‑1591475

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: