Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,904,593 C→T 27.7% K443K (AAG→AAA)  ptr ← protease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,904,5930CT27.7% 115.1 / 42.1 104K443K (AAG→AAA) ptrprotease III
Reads supporting (aligned to +/- strand):  ref base C (17/58);  new base T (14/15);  total (31/73)
Fisher's exact test for biased strand distribution p-value = 1.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 109 total aligned reads displayed.

CTTTCGGGCTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGA  >  minE/1904535‑1904656
                                                          |                                                               
ctTTCGGGCTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCttt                                                               >  1:310758/1‑61 (MQ=255)
ctTTCGGGCTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCttt                                                               >  1:1420755/1‑61 (MQ=255)
ctTTCGGGCTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCttt                                                               >  1:1364804/1‑61 (MQ=255)
     gggCTGATATACCAGATACGCGCATTCTGCGTCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAg                                                    <  1:1359335/67‑1 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:520513/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:2999678/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:663523/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:275668/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:750429/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:2150232/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:900444/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:210073/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:2011210/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:1963509/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:1142884/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGCATCg                                               >  1:1142883/1‑69 (MQ=255)
        cTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCTTTTACTGCTTTAGAATCg                                               >  1:2795146/1‑69 (MQ=255)
         tGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGt                                              >  1:3229343/1‑69 (MQ=255)
              taCCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGa                                         >  1:1196815/1‑69 (MQ=255)
              taCCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGa                                         >  1:1432070/1‑69 (MQ=255)
              taCCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGa                                         >  1:2242171/1‑69 (MQ=255)
              taCCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGa                                         >  1:145597/1‑69 (MQ=255)
              aaCCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGa                                         >  1:3071423/2‑69 (MQ=255)
                ccAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATc                                       >  1:1060979/1‑69 (MQ=255)
                ccAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATc                                       >  1:1690022/1‑69 (MQ=255)
                ccAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATc                                       >  1:618923/1‑69 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCg                                          >  1:1923359/1‑61 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:3187475/69‑1 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:1731996/69‑1 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:1433540/69‑1 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:1907299/69‑1 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:2189628/69‑1 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:2165234/69‑1 (MQ=255)
                     tACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:139651/69‑1 (MQ=255)
                     aaCGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:1890682/68‑1 (MQ=255)
                      aCGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCaa                                  <  1:677441/68‑1 (MQ=255)
                          gcATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTg                             <  1:3076798/69‑1 (MQ=255)
                          gcATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTg                             <  1:441744/69‑1 (MQ=255)
                          gcATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTg                             <  1:2659820/69‑1 (MQ=255)
                          gcATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTg                             <  1:3299639/69‑1 (MQ=255)
                          gcATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTg                             <  1:1411299/69‑1 (MQ=255)
                          gcATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTg                             <  1:1665735/69‑1 (MQ=255)
                               cTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGc                        <  1:3198278/69‑1 (MQ=255)
                                tGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGc                        <  1:3086216/68‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGGTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:195200/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:1173370/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:616156/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:3362960/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:1444342/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:2207518/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:1001537/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:2701996/69‑1 (MQ=255)
                                       tcatcaTCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCg                <  1:2677696/69‑1 (MQ=255)
                                         atcatcGCCAGACGTTCCTTTACTGCTTTAGCATCg                                               >  1:1348568/1‑36 (MQ=255)
                                         atcatcGCCAGACGTTCCTTTACTGCTTTAGCATCg                                               >  1:1813369/1‑36 (MQ=255)
                                         atcatcGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGt               >  1:1059595/1‑68 (MQ=255)
                                         atcatcGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTa              <  1:2063489/69‑1 (MQ=255)
                                           catcGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGAGCGGCAATATTGACTGCATCCAGCGTa              <  1:2927912/67‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:2564530/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:2111430/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:2056110/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:541446/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:1213303/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:3024567/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:1044853/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:2656669/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:1805471/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:2560398/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:2558007/69‑1 (MQ=255)
                                              cGCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:1675183/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:1850288/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:2573738/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:385488/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:1984706/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:2204279/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:1481960/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:1540019/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:1089928/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:1698009/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:2770252/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:2938944/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:2794330/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:2853195/69‑1 (MQ=255)
                                              cGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTc         <  1:2840186/69‑1 (MQ=255)
                                               gCCAGACGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:659055/68‑1 (MQ=255)
                                                    aCGTTCTTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGTGTATGCTc         <  1:521076/63‑1 (MQ=255)
                                                    aCGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAgg   <  1:905692/69‑1 (MQ=255)
                                                    aCGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAgg   <  1:912055/69‑1 (MQ=255)
                                                    aCGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAgg   >  1:2317849/1‑69 (MQ=255)
                                                    aCGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAgg   >  1:3406132/1‑69 (MQ=255)
                                                    aCGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAgg   >  1:3312512/1‑69 (MQ=255)
                                                    aCGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAgg   <  1:2440891/69‑1 (MQ=255)
                                                    aCGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAgg   <  1:2135563/69‑1 (MQ=255)
                                                     cGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGa  <  1:454852/69‑1 (MQ=255)
                                                     cGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGa  <  1:2885854/69‑1 (MQ=255)
                                                     cGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGa  <  1:2656396/69‑1 (MQ=255)
                                                     cGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGa  <  1:1189296/69‑1 (MQ=255)
                                                     cGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGa  <  1:1035454/69‑1 (MQ=255)
                                                         ccTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAAc      >  1:2810378/1‑61 (MQ=255)
                                                         ccTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAAc      >  1:666548/1‑61 (MQ=255)
                                                          |                                                               
CTTTCGGGCTGATATACCAGATACGCGCATTCTGCGGCGTCATCATCGCCAGACGTTCCTTTACTGCTTTAGCATCGTACCGATCGGCAATATTGACTGCATCCAGCGTATGCTCAACAGGA  >  minE/1904535‑1904656

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