Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,585,072 C→T 35.1% C184C (TGC→TGT cysZ → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,585,0720CT35.1% 94.3 / 71.7 125C184C (TGC→TGTcysZpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (66/15);  new base T (15/29);  total (81/44)
Fisher's exact test for biased strand distribution p-value = 1.88e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.30e-01

Only 100 of 125 total aligned reads displayed.

TTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAA  >  minE/1585005‑1585138
                                                                   |                                                                  
ttGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCg                                                                   >  1:1131003/1‑69 (MQ=255)
ttGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCg                                                                   >  1:1641926/1‑69 (MQ=255)
  gggCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGAt                                                                 >  1:1588741/1‑69 (MQ=255)
  gggCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGAt                                                                 >  1:940042/1‑69 (MQ=255)
         ccGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTAcccc                                                           >  1:3308598/1‑68 (MQ=255)
         ccGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCt                                                          >  1:2176217/1‑69 (MQ=255)
         ccGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCt                                                          >  1:3322329/1‑69 (MQ=255)
         ccGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCt                                                          >  1:318452/1‑69 (MQ=255)
         ccGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCt                                                          >  1:1453341/1‑69 (MQ=255)
         ccGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCt                                                          >  1:3036939/1‑69 (MQ=255)
         ccGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCt                                                          >  1:2017845/1‑69 (MQ=255)
              gcgcCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                                                     >  1:2364835/1‑69 (MQ=255)
              gcgcCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                                                     >  1:2336302/1‑69 (MQ=255)
              gcgcCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                                                     >  1:1350672/1‑69 (MQ=255)
              gcgcCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                                                     >  1:1641430/1‑69 (MQ=255)
              gcgcCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                                                     >  1:3208341/1‑69 (MQ=255)
              gcgcCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                                                     >  1:1478629/1‑69 (MQ=255)
              gcgcCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                                                     >  1:1047002/1‑69 (MQ=255)
                      aCTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACaa                                             >  1:2742857/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAAcc                                                 >  1:1692855/1‑62 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:1843183/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:3299681/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:1784720/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:178856/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:399581/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:425672/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:1157959/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:2875304/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:1137968/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:2405497/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:2461547/1‑69 (MQ=255)
                         gtgGTTCCTGTTTAGCGCCTGGATGTTAGCAATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGt                                          >  1:2161029/1‑69 (MQ=255)
                          ttgTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:1518863/3‑69 (MQ=255)
                          tgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:3098977/1‑69 (MQ=255)
                          tgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:3026401/1‑69 (MQ=255)
                          tgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:2658212/1‑69 (MQ=255)
                          tgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:1019707/1‑69 (MQ=255)
                          tgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:537915/1‑69 (MQ=255)
                          tgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:118854/1‑69 (MQ=255)
                          tgGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTg                                         >  1:1132144/1‑69 (MQ=255)
                               ccTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTa                                                               >  1:199648/1‑42 (MQ=255)
                               ccTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGt                                     >  1:1143951/1‑68 (MQ=255)
                               ccTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGt                                     >  1:1965072/1‑68 (MQ=255)
                               ccTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGt                                     >  1:2895140/1‑68 (MQ=255)
                               ccTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCAATTACCCCTTCGATAACCACAAAGTGCCGt                                     >  1:44354/1‑68 (MQ=255)
                                  gTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGtt                                    >  1:2075557/1‑66 (MQ=255)
                                  gTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGtt                                    >  1:2661034/1‑66 (MQ=255)
                                  atgtAGCGCCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:1065851/66‑1 (MQ=255)
                                     tAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:2765050/1‑65 (MQ=255)
                                     tAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:2316069/1‑65 (MQ=255)
                                     tAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:783434/1‑65 (MQ=255)
                                     tAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:1987707/1‑65 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:967360/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:282980/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:862977/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:3003015/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:655109/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:556689/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:520498/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:330787/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:3246443/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:1771730/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:1030783/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:116935/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:1267226/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:1373899/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:163215/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:2089010/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:2130346/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:2174583/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:2542051/62‑1 (MQ=255)
                                         gcCTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:2535825/62‑1 (MQ=255)
                                          ccTGGATGTTAGCCATCCAGTATTGTGATTACCCCTTCGATAACCACAAAGTGCCGTTTaa                                 <  1:2622875/61‑1 (MQ=255)
                                            tGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATgcgc                       <  1:1821668/69‑1 (MQ=255)
                                                  ttAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGcc                 >  1:3243980/1‑69 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:1618082/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:1114475/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:90966/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:2692772/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:817099/1‑50 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:2541770/1‑50 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:1364463/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:1556956/1‑50 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:2795921/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:2288402/1‑50 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:2371761/1‑50 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:2381673/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:3054943/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:2099698/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:2096533/50‑1 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  >  1:2810363/1‑50 (MQ=255)
                                                    aGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTa                                  <  1:1618726/50‑1 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCa                      >  1:1453427/1‑48 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTg              >  1:1438735/1‑56 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCaaa  >  1:1618305/1‑68 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCaaa  >  1:706479/1‑68 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCaaa  >  1:1566396/1‑68 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCaaa  >  1:408740/1‑68 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCaaa  >  1:2347853/1‑68 (MQ=255)
                                                                  gCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCaaa  >  1:388658/1‑68 (MQ=255)
                                                                   |                                                                  
TTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGTTTAAAGAGATGCGCACCGCCCTGCGCACACGCAAA  >  minE/1585005‑1585138

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: