Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,049 C→T 48.1% N12N (AAC→AAT yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,0490CT48.1% 11.8 / 145.1 106N12N (AAC→AATyafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base T (30/21);  total (57/49)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 106 total aligned reads displayed.

TAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATT  >  minE/217986‑218114
                                                               |                                                                 
   gCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                                         <  1:2540746/71‑1 (MQ=255)
    cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                                         >  1:1681841/1‑70 (MQ=255)
     tGCACTGGCGTCTGGACGTGGTAATGAATGAAGCCGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatta                                                       >  1:2413801/1‑71 (MQ=255)
     tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAatta                                                       >  1:3952469/1‑71 (MQ=255)
     tGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatt                                                        >  1:1169153/1‑70 (MQ=255)
       cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattat                                                      <  1:616730/70‑1 (MQ=255)
         ctggCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTAtttt                                                   >  1:299572/3‑71 (MQ=255)
          aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGcagcag                                                            <  1:3086400/61‑1 (MQ=255)
          aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                   >  1:2538903/1‑70 (MQ=255)
            gCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                   >  1:309878/1‑68 (MQ=255)
            gCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt                                                   >  1:3619811/1‑68 (MQ=255)
              gTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa                                                          >  1:3437745/1‑59 (MQ=255)
                cTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGa                                             >  1:334325/1‑70 (MQ=255)
                 tGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                                  >  1:2199733/1‑64 (MQ=255)
                 tGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                           >  1:3495409/1‑71 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATgcagc                                                              >  1:3611270/1‑50 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACg                                         >  1:535588/1‑71 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAat                                                         <  1:2514967/55‑1 (MQ=255)
                   gACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATa                                           <  1:1624961/69‑1 (MQ=255)
                    aCGTGTTAATTAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACgg                                        <  1:750625/71‑1 (MQ=255)
                     cGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                                                 >  1:1290340/1‑61 (MQ=255)
                       tGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGca                                      >  1:3189581/1‑70 (MQ=255)
                        ggTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGcaca                                    >  1:3813831/1‑71 (MQ=255)
                        ggTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGcaca                                    >  1:907077/1‑71 (MQ=255)
                         gTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAt                                   <  1:41141/71‑1 (MQ=255)
                          tAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACAtt                                  >  1:2290569/1‑71 (MQ=255)
                          aaatgaatgAAGATGCCAGCCGGATAAGAAGAGGAAACgca                                                                <  1:1021342/40‑1 (MQ=255)
                           aatgaatgAAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTg                                 >  1:3231337/1‑71 (MQ=255)
                            atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCa                                                 <  1:757557/54‑1 (MQ=255)
                            atgaatgaAGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                >  1:3676027/1‑55 (MQ=255)
                            atgaatgaAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGc                                <  1:1563068/71‑1 (MQ=255)
                             tgaatgaaGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                               <  1:1235152/71‑1 (MQ=255)
                               aatgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                >  1:1358576/1‑52 (MQ=255)
                                atgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGAACATTGCTAt                             >  1:602004/1‑70 (MQ=255)
                                atgaaGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAatta                                                       <  1:2802961/44‑1 (MQ=38)
                                 tgaaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAg                                                <  1:24334/50‑1 (MQ=39)
                                  gaaGACGACTACAAAATAAGAAGAGGAAACGcagcag                                                            <  1:470995/37‑1 (MQ=255)
                                   aaGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTaa                          >  1:133170/1‑70 (MQ=255)
                                    cgacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGAt                                            <  1:567610/50‑1 (MQ=38)
                                    cgacgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTc                                                  <  1:2962830/44‑1 (MQ=38)
                                    cgacgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATAc                                          <  1:1144971/52‑1 (MQ=255)
                                    aGACGACTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAgg                                               <  1:814196/48‑1 (MQ=39)
                                    aGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACg                                         <  1:3830853/54‑1 (MQ=255)
                                    aGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAattatt                                                     <  1:1573422/42‑1 (MQ=255)
                                     gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                       >  1:3281020/1‑71 (MQ=255)
                                     gacgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                       >  1:121127/1‑71 (MQ=255)
                                     gacgacTACAAAATAAGAAGAGGGAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                       >  1:492523/1‑71 (MQ=255)
                                      acgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTc                                                  <  1:2079541/43‑1 (MQ=38)
                                      acgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGa                                             >  1:204121/1‑48 (MQ=255)
                                      acgacTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAatat                       >  1:1630863/1‑70 (MQ=255)
                                       cgacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAttt                     >  1:125265/1‑71 (MQ=255)
                                       cgacTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCa                                                 <  1:1548989/43‑1 (MQ=38)
                                        gacTGCAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                    >  1:1274865/1‑71 (MQ=255)
                                         acTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATAtttt                    <  1:2097155/70‑1 (MQ=255)
                                           tgcAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGca                                      <  1:1160048/48‑1 (MQ=39)
                                           tACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCa                                                 >  1:863718/1‑39 (MQ=255)
                                             cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGGACATTGCTATTAATATTTTGACg                <  1:1249133/70‑1 (MQ=255)
                                             cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                >  1:3446679/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                >  1:2679737/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGAt                                            <  1:1649140/42‑1 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                >  1:1025623/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                <  1:26962/70‑1 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                >  1:1587024/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACg                >  1:1962411/1‑70 (MQ=255)
                                             cAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGa               <  1:3178782/71‑1 (MQ=255)
                                              aaaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGaa              >  1:2943575/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCa                                                 <  1:2566579/35‑1 (MQ=255)
                                               aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt             >  1:1162456/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt             >  1:1282403/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt             >  1:2080023/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt             >  1:2635144/1‑71 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                               >  1:3244053/1‑53 (MQ=255)
                                               aaaTAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCt                               >  1:3248726/1‑53 (MQ=255)
                                                aaTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt             >  1:714974/1‑70 (MQ=255)
                                                 aTAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATga           >  1:746633/1‑71 (MQ=255)
                                                  tAAGAAGAGGAAATGCAGCAGAATTATTTTCAGGGATa                                           <  1:2759227/38‑1 (MQ=255)
                                                  tAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgag          <  1:1561891/71‑1 (MQ=255)
                                                   aagaagAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgag          >  1:1398727/1‑70 (MQ=255)
                                                    agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgata         <  1:3797051/70‑3 (MQ=255)
                                                    agaagaGGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATgata         <  1:629809/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:3544892/1‑70 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg      >  1:1311194/1‑70 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAgg      >  1:3544019/1‑70 (MQ=255)
                                                       agagGAAATGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATAAGGt     <  1:2237575/71‑1 (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:1810703/43‑1 (MQ=255)
                                                            aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTAtt                            <  1:1234284/43‑1 (MQ=255)
                                                            aaaTGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGATaa        >  1:3004111/1‑63 (MQ=255)
                                                            aaaCGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAAt             <  1:3268993/58‑1 (MQ=255)
                                                               |                                                                 
TAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTTCAGGGATACGGCACATTGCTATTAATATTTTGACGAATGAGAAGGTATT  >  minE/217986‑218114

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: