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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0820TC10.8% 300.7 / 22.4 111I23I (ATT→ATCyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (52/47);  new base C (5/7);  total (57/54)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 120 total aligned reads displayed.

TAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCA  >  minE/218012‑218152
                                                                      |                                                                      
tAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAtt                                                                        >  1:2290569/1‑71 (MQ=255)
 aatgaatgAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                                                       >  1:3231337/1‑71 (MQ=255)
  atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                                                      <  1:1563068/71‑1 (MQ=255)
      atgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGAACATTGCTAt                                                                   >  1:602004/1‑70 (MQ=255)
           gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                                                             >  1:3281020/1‑71 (MQ=255)
           gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                                                             >  1:121127/1‑71 (MQ=255)
           gacgacTACAAAATAAGAAGAGGGAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                                                             >  1:492523/1‑71 (MQ=255)
            acgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                                                             >  1:1630863/1‑70 (MQ=255)
             cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAttt                                                           >  1:125265/1‑71 (MQ=255)
              gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                                                          <  1:3392271/71‑1 (MQ=255)
              gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                                                          >  1:1274865/1‑71 (MQ=255)
                   cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGGACATTGCTATTAATATTTTGACg                                                      <  1:1249133/70‑1 (MQ=255)
                   cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                                                      >  1:2679737/1‑70 (MQ=255)
                   cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                                                      <  1:26962/70‑1 (MQ=255)
                   cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                                                      >  1:1025623/1‑70 (MQ=255)
                   cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                                                      >  1:1962411/1‑70 (MQ=255)
                   cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                                                      >  1:1587024/1‑70 (MQ=255)
                    aaaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGaa                                                    >  1:2943575/1‑71 (MQ=255)
                     aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                                                   >  1:1282403/1‑71 (MQ=255)
                     aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                                                   >  1:1162456/1‑71 (MQ=255)
                     aaaTAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt                                                   >  1:2080023/1‑71 (MQ=255)
                     aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                                                                     >  1:3248726/1‑53 (MQ=255)
                     aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                                                                     >  1:3244053/1‑53 (MQ=255)
                      aaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                                                   >  1:714974/1‑70 (MQ=255)
                       aTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATga                                                 >  1:746633/1‑71 (MQ=255)
                         aagaagAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgag                                                >  1:1398727/1‑70 (MQ=255)
                          agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgata                                               <  1:3797051/70‑3 (MQ=255)
                          agaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga                                               <  1:35170/70‑1 (MQ=255)
                           gaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAat                                                               <  1:3589853/53‑1 (MQ=255)
                           gaagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAg                                             <  1:3688200/71‑1 (MQ=255)
                            aagaGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACAtt                                                                        <  1:2948928/43‑1 (MQ=255)
                             agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                                                       <  1:3535227/43‑1 (MQ=255)
                             agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                                                       <  1:83172/43‑1 (MQ=255)
                             agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg                                            >  1:1311194/1‑70 (MQ=255)
                             agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg                                            >  1:3544892/1‑70 (MQ=255)
                             agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGt                                           <  1:2237575/71‑1 (MQ=255)
                              gagGAAATGCAGAAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTa                                          >  1:468475/1‑71 (MQ=255)
                               agGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                                                       <  1:2126842/41‑1 (MQ=255)
                               agGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAt                                         <  1:148946/71‑1 (MQ=255)
                               agGAAATGCAGCAGAATTATTATCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTAt                                         >  1:2042261/1‑71 (MQ=255)
                               agGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAt                                         >  1:2218775/1‑71 (MQ=255)
                                ggAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTAtt                                        >  1:590098/1‑71 (MQ=255)
                                 gAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTa                                          >  1:2488924/1‑68 (MQ=255)
                                  aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt                                                                  <  1:1234284/43‑1 (MQ=255)
                                  aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATaa                                              >  1:3004111/1‑63 (MQ=255)
                                  aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt                                                   <  1:3268993/58‑1 (MQ=255)
                                  aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGt                                           <  1:2092416/66‑1 (MQ=255)
                                   aatgcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTc                                       >  1:1713405/4‑69 (MQ=255)
                                     tgcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAgg                                   <  1:397642/70‑1 (MQ=255)
                                      gcagcaGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGc                                  <  1:357132/71‑1 (MQ=255)
                                        agcagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAg                                >  1:3221078/1‑71 (MQ=255)
                                          cagAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                                                          <  1:2420735/43‑1 (MQ=255)
                                          cagAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAggg                              >  1:3282911/1‑71 (MQ=255)
                                           agAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGt                             >  1:3746727/1‑71 (MQ=255)
                                           agAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTAtt                                        <  1:2851185/60‑1 (MQ=255)
                                            gAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTg                                                         <  1:3005782/42‑1 (MQ=255)
                                            gAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                                                       >  1:514160/1‑44 (MQ=255)
                                            gAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                                                      >  1:3291412/1‑45 (MQ=255)
                                            gAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGt                             >  1:3750623/1‑70 (MQ=255)
                                             aaTTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGtt                            >  1:106149/1‑70 (MQ=255)
                                               ttattTTCAGGGATACGGCACATTGCTATTAatat                                                             <  1:1912676/35‑1 (MQ=255)
                                               ttattTTCAGGGATACGGCACATTGCTATTAATAtttt                                                          >  1:360415/1‑38 (MQ=255)
                                               ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                                                       >  1:910615/1‑41 (MQ=255)
                                               ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGAc                                                       >  1:3902723/1‑41 (MQ=255)
                                               ttattTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgaga                                               <  1:2267441/49‑1 (MQ=255)
                                               ttatcTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGtt                            >  1:1096022/1‑68 (MQ=255)
                                                 attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGAc                       >  1:3921119/1‑71 (MQ=255)
                                                 attTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGAc                       >  1:2504320/1‑71 (MQ=255)
                                                  ttttCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGAc                       <  1:1024937/70‑1 (MQ=255)
                                                    ttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                     <  1:1830179/70‑1 (MQ=255)
                                                    ttCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGt                     <  1:3534068/70‑1 (MQ=255)
                                                    ttCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                     <  1:618225/70‑1 (MQ=255)
                                                     tCAGGGATACGGCACATTGCTATTAATATTTTGACg                                                      >  1:3424510/1‑36 (MQ=255)
                                                     tCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGt                     >  1:814337/1‑69 (MQ=255)
                                                     tCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTa                    <  1:709753/70‑1 (MQ=255)
                                                       aGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCa                                      <  1:346821/50‑1 (MQ=255)
                                                        gggATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTaa                          >  1:1177408/1‑61 (MQ=255)
                                                        gggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCa                                 <  1:1375320/54‑1 (MQ=255)
                                                         ggATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAg                         >  1:1663374/1‑61 (MQ=255)
                                                         ggATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGc                                  <  1:1538305/52‑1 (MQ=255)
                                                          gATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGc              <  1:3494260/71‑1 (MQ=255)
                                                           aTACGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCg             <  1:2167799/71‑1 (MQ=255)
                                                             aCGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCATGGTTAAGACGTAAGATGCGa            <  1:905533/70‑1 (MQ=255)
                                                             aCGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGa            >  1:62942/1‑70 (MQ=255)
                                                              cGGCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaa           <  1:2105306/70‑1 (MQ=255)
                                                              cGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaaa          <  1:1962751/71‑1 (MQ=255)
                                                               ggCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaa           <  1:715268/69‑1 (MQ=255)
                                                                gCACATTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag        <  1:3701228/71‑1 (MQ=255)
                                                                gCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag        >  1:2652520/1‑71 (MQ=255)
                                                                 cacaTTGCTATTATTATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAag        >  1:91430/1‑70 (MQ=255)
                                                                 cacaTTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagc       <  1:2067929/71‑1 (MQ=255)
                                                                   caTCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTaa                   <  1:1048425/57‑1 (MQ=255)
                                                                   caTCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTaa                   <  1:1898900/57‑1 (MQ=255)
                                                                    aTTGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACg                      <  1:3178609/53‑1 (MQ=255)
                                                                     ttGCTATTAATATTTTGACGAATGATAAGGTATTCAAgg                                   <  1:1036586/39‑1 (MQ=255)
                                                                     ttGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGcc   <  1:1390508/71‑1 (MQ=255)
                                                                     ttGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACg                      >  1:3125694/1‑52 (MQ=255)
                                                                      tGCTATTAATATTTTGACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCa  >  1:3720397/1‑71 (MQ=255)
                                                                      tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGAt                >  1:845740/1‑57 (MQ=255)
                                                                      tGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGAt                >  1:2014333/1‑57 (MQ=255)
                                                                      |                                                                      
TAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCA  >  minE/218012‑218152

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: