Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,082 T→C 9.0% I23I (ATT→ATC yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0820TC 9.0% 414.5 / 26.2 145I23I (ATT→ATCyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (71/61);  new base C (8/5);  total (79/66)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 157 total aligned reads displayed.

ATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCA  >  minE/218014‑218152
                                                                    |                                                                      
atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                                                      >  1:889177/1‑71 (MQ=255)
atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                                                      >  1:3557409/1‑71 (MQ=255)
  gaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTa                                                                    >  1:1154347/1‑71 (MQ=255)
   aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                                                   >  1:1118129/1‑71 (MQ=255)
   aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                                                   >  1:247012/1‑71 (MQ=255)
   aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                                                   >  1:2222208/1‑71 (MQ=255)
   aatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAt                                                                   >  1:2690485/1‑71 (MQ=255)
    atgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt                                                                   >  1:2960934/1‑70 (MQ=255)
     tgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt                                                                   >  1:884783/1‑69 (MQ=255)
     tgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt                                                                  <  1:3743650/70‑1 (MQ=255)
       aaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                                                                <  1:1125490/70‑1 (MQ=255)
        aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTAAGGGATACGGCACATCGCTAtt                                                                  <  1:971237/67‑1 (MQ=255)
         gacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                                                             <  1:3859721/71‑1 (MQ=255)
               tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                                                                     >  1:3523939/1‑57 (MQ=255)
               tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                                                          <  1:3558693/68‑1 (MQ=255)
               tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                                                       >  1:1412998/1‑71 (MQ=255)
               tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                                                       >  1:1151283/1‑71 (MQ=255)
               tACAAAATAAGAAGAGGAAACCCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                                                       >  1:3452223/1‑71 (MQ=255)
                 cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                                                      <  1:4025700/70‑1 (MQ=255)
                 cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAata                                                              >  1:827738/1‑62 (MQ=255)
                 cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                                                      >  1:3335661/1‑70 (MQ=255)
                 cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                                                      >  1:1544234/1‑70 (MQ=255)
                 cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                                                      >  1:2184564/1‑70 (MQ=255)
                  aaaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGa                                                     <  1:1138848/70‑1 (MQ=255)
                   aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                                                   >  1:3252562/1‑71 (MQ=255)
                      tAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGTTATTaa                                                                <  1:3987488/55‑1 (MQ=255)
                      tAAGAAGAGGAAACGCAGCAGAATAATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                                                   >  1:2261194/1‑68 (MQ=255)
                       aagaagAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt                                                                  <  1:2382889/52‑1 (MQ=255)
                        agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgata                                               <  1:3397028/70‑3 (MQ=255)
                        agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgat                                                >  1:2761210/1‑68 (MQ=255)
                        agaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga                                               <  1:3523578/70‑1 (MQ=255)
                          aagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg                                            <  1:716141/71‑1 (MQ=255)
                           agagGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAgg                                            >  1:1454612/1‑70 (MQ=255)
                           agagGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg                                            <  1:1554768/70‑1 (MQ=255)
                             agGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTCGACGAATGATAAGGTAt                                         >  1:530055/1‑71 (MQ=255)
                             agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt                                         <  1:1766001/71‑1 (MQ=255)
                             agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt                                         <  1:2711932/71‑1 (MQ=255)
                              ggAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg                                            <  1:2449708/67‑1 (MQ=255)
                              ggAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTAt                                         <  1:1519628/70‑1 (MQ=255)
                              ggAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGat                                                    <  1:3280607/59‑2 (MQ=255)
                               gAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc                                       <  1:1069743/71‑1 (MQ=255)
                               gAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc                                       <  1:1858237/71‑1 (MQ=255)
                               gAAACGCAGCAGAATTATTTTCAGGGATACGGAACAtt                                                                        <  1:846615/38‑1 (MQ=255)
                                aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc                                       <  1:2087874/70‑1 (MQ=255)
                                aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc                                       >  1:2414334/1‑70 (MQ=255)
                                 aaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAgg                                            >  1:4120883/1‑64 (MQ=255)
                                  aCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGaa                                              >  1:2303880/1‑61 (MQ=255)
                                    gcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGa                                                     <  1:1337199/52‑1 (MQ=255)
                                      agcagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                                                   >  1:3331619/1‑52 (MQ=255)
                                      agcagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAg                                <  1:3959747/71‑1 (MQ=255)
                                          gAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGt                             >  1:3909818/1‑70 (MQ=255)
                                          gAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGt                             >  1:795900/1‑70 (MQ=255)
                                            attattTTCAGGGATACGGCACATTGCTATTAATAtttt                                                          <  1:2316262/39‑1 (MQ=255)
                                            attattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCGGGGtt                            >  1:3636554/1‑69 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGaa                                                    >  1:561752/1‑44 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGaa                                                    >  1:1262472/1‑44 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGa                                                     >  1:1842020/1‑43 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGtt                            >  1:68016/1‑68 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGt                             >  1:2042152/1‑67 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTAAAGGCAGGGtt                            <  1:2115687/68‑1 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTa                                          >  1:3322569/1‑54 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTa                                          >  1:428243/1‑54 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTc                                       <  1:2369212/57‑1 (MQ=255)
                                             ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAg                         <  1:895905/71‑1 (MQ=255)
                                               attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAg                                             >  1:1490770/1‑49 (MQ=255)
                                               attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGCTAAGa                        <  1:1009832/70‑1 (MQ=255)
                                                ttttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTaa                          <  1:1928564/67‑1 (MQ=255)
                                                tttaCAGGGATAGCGCACATTGCTATTAATATTTTGACGAATGAGaact                                            >  1:3174574/1‑47 (MQ=37)
                                                   tCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAggg                              >  1:2709848/1‑60 (MQ=255)
                                                   tCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTa                    <  1:1200782/70‑1 (MQ=255)
                                                    cAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTaa                   <  1:2037064/70‑1 (MQ=255)
                                                    cAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGtt                            <  1:1082554/61‑1 (MQ=255)
                                                    cAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTaa                   <  1:1493411/70‑1 (MQ=255)
                                                     aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAg                                             <  1:4058393/43‑1 (MQ=255)
                                                     aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAgg                                   <  1:2914731/53‑1 (MQ=255)
                                                       ggATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc                                       <  1:2658287/47‑1 (MQ=255)
                                                         aTACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc                                       >  1:3922228/1‑45 (MQ=255)
                                                         aTACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc                                       >  1:259644/1‑45 (MQ=255)
                                                          tACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTaa                   <  1:1219158/64‑1 (MQ=255)
                                                           aCGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGa            >  1:1944323/1‑70 (MQ=255)
                                                           aCGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGa            >  1:168842/1‑70 (MQ=255)
                                                            cGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTaa                          >  1:894209/1‑55 (MQ=255)
                                                            cGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaa           >  1:3809139/1‑70 (MQ=255)
                                                            cGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTaa                          >  1:924637/1‑55 (MQ=255)
                                                            cGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTaa                          >  1:2223755/1‑55 (MQ=255)
                                                            cGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTaa                   >  1:2026056/1‑62 (MQ=255)
                                                            cGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaa           >  1:3708275/1‑70 (MQ=255)
                                                            cGGCACATTGCTATTAATATTTCGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaa           >  1:2273019/1‑70 (MQ=255)
                                                                acaTTGCTATTAATATTTTGACGAATGATAAGGTATTCa                                      <  1:1805061/39‑1 (MQ=255)
                                                                acaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagca      <  1:238345/71‑1 (MQ=255)
                                                                acaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag        >  1:471657/1‑69 (MQ=255)
                                                                acaTCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagca      <  1:3211342/71‑1 (MQ=255)
                                                                 caTCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagcag     <  1:1212147/71‑1 (MQ=255)
                                                                  aTTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagcagc    >  1:1240525/1‑71 (MQ=255)
                                                                  aTTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagcagc    >  1:974694/1‑71 (MQ=255)
                                                                  aTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGtt                            <  1:393297/47‑1 (MQ=255)
                                                                    tGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCa  >  1:1233439/1‑71 (MQ=255)
                                                                    tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGa                        >  1:1432645/1‑49 (MQ=255)
                                                                    tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGa                        >  1:1659491/1‑49 (MQ=255)
                                                                    tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag        >  1:348068/1‑65 (MQ=255)
                                                                    |                                                                      
ATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCA  >  minE/218014‑218152

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: