Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,590,902 A→G 45.9% V32V (GTA→GTG yfeS → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,590,9020AG45.9% 9.9 / 95.5 73V32V (GTA→GTGyfeSconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (13/26);  new base G (12/21);  total (25/48)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAGTAAA  >  minE/1590834‑1590971
                                                                    |                                                                     
gAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtggg                                                                     >  1:4209920/1‑71 (MQ=255)
gAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtagg                                                                     >  1:4467478/1‑71 (MQ=255)
 aaaaaTCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtagg                                                                     <  1:3579489/70‑1 (MQ=255)
  aaaaTCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtggg                                                                     >  1:3271782/1‑69 (MQ=255)
  aaaaTCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtggg                                                                     >  1:4502171/1‑69 (MQ=255)
  aaaaTCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtaggt                                                                    <  1:3449448/70‑1 (MQ=255)
   aaaTCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGgtgggtag                                                                  <  1:4227057/71‑1 (MQ=255)
      tCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTAt                                                               <  1:3950855/71‑1 (MQ=255)
       cGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTAtt                                                              >  1:4176437/1‑71 (MQ=255)
          aTAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTgg                                                            >  1:4084111/1‑70 (MQ=255)
             aaTTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTaaa                                                        >  1:1646546/1‑71 (MQ=255)
              aTTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAAt                                                       >  1:2321281/1‑71 (MQ=255)
                 ttGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTc                                                     >  1:4015058/1‑70 (MQ=255)
                  tggtggATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTc                                                     <  1:3634587/69‑1 (MQ=255)
                   ggtggATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCa                                                   <  1:924084/70‑1 (MQ=255)
                   ggtggATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCa                                                   <  1:2622406/70‑1 (MQ=255)
                    gtggATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAg                                                  <  1:1643976/70‑1 (MQ=255)
                     tggATCGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGgtagt                                                                 <  1:3577136/54‑1 (MQ=255)
                     tggATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCa                                                   <  1:590634/68‑1 (MQ=255)
                     tggATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGgtagt                                                                 <  1:501285/54‑1 (MQ=255)
                       gATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaa                                              <  1:2621989/71‑1 (MQ=255)
                       gATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaa                                              <  1:2526369/71‑1 (MQ=255)
                       gATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaa                                              <  1:2352675/71‑1 (MQ=255)
                       gATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaa                                              <  1:1624782/71‑1 (MQ=255)
                        aTAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaaa                                             <  1:1521097/71‑1 (MQ=255)
                        aTAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaa                                              >  1:3326123/1‑70 (MQ=255)
                         tAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaaaa                                            <  1:69325/71‑1 (MQ=255)
                         tAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaaaa                                            <  1:3702172/71‑1 (MQ=255)
                         tAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaaaa                                            <  1:1376950/71‑1 (MQ=255)
                         tAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaaaa                                            <  1:4457593/71‑1 (MQ=255)
                          aGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaa                                              >  1:2116697/1‑68 (MQ=255)
                           gATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAag                                           <  1:3403504/70‑1 (MQ=255)
                            aTTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaaaa                                            <  1:1154887/68‑1 (MQ=255)
                            aTTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAagag                                         <  1:3758853/71‑1 (MQ=255)
                            aTTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAagag                                         <  1:3275689/71‑1 (MQ=255)
                            aTTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAagag                                         >  1:1427849/1‑71 (MQ=255)
                             ttACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtggg                                                                     <  1:974542/42‑1 (MQ=255)
                             ttACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaaaa                                            <  1:741056/67‑1 (MQ=255)
                              tACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtagg                                                                     >  1:1282880/1‑41 (MQ=255)
                               aCGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCa                                                   >  1:2474176/1‑58 (MQ=255)
                               aCGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaa                                              >  1:371800/1‑63 (MQ=255)
                               aCGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACa                                               >  1:91967/1‑62 (MQ=255)
                                cGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTa                                                          <  1:1556951/50‑1 (MQ=255)
                                cGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTc                                      <  1:4646176/70‑1 (MQ=255)
                                 gAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTa                                                          <  1:3581872/49‑1 (MQ=255)
                                 gAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGa                                    >  1:3649960/1‑71 (MQ=255)
                                 gAAGGGGATAGTTTAGCTCTCAACTATGGCAAGGTGGGTAGTATTGGTa                                                          <  1:2058844/49‑1 (MQ=255)
                                  aaGGGGATAGTTTAGCTGTCAAATATGGCAAGGTTTGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGAt                                   <  1:2612532/71‑1 (MQ=255)
                                    ggggATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaaa                                             >  1:3784249/1‑59 (MQ=255)
                                       gATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGt                                                           <  1:126258/42‑1 (MQ=255)
                                       gATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATg                               >  1:2763074/1‑70 (MQ=255)
                                        aTAGTTTAGCTGTCAATTATGGCTTGGTAGGTAGTATTTGTa                                                          <  1:3364096/42‑1 (MQ=255)
                                         tAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaat                            <  1:454959/71‑2 (MQ=255)
                                         tAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaag                            <  1:1494811/71‑1 (MQ=255)
                                         tAGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaag                            <  1:4477492/71‑1 (MQ=255)
                                         tAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCa                                                   <  1:4674221/48‑1 (MQ=255)
                                         tAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCa                                                   <  1:3920091/48‑1 (MQ=255)
                                         tAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaag                            <  1:115241/71‑1 (MQ=255)
                                          aGTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaag                            >  1:3252030/1‑70 (MQ=255)
                                           gTTTAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGAt                                   <  1:361394/62‑1 (MQ=255)
                                             ttAGCTGTCAACTATGGCAAGGTGGGTAGTATTgg                                                            >  1:4545354/1‑35 (MQ=255)
                                             ttAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaagaa                          >  1:3460046/1‑69 (MQ=255)
                                              tAGCTGTCAACTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATaa                                 >  1:240854/1‑61 (MQ=255)
                                              tAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATa                                  >  1:1900422/1‑60 (MQ=255)
                                                gCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAgtg                     <  1:1709845/71‑1 (MQ=255)
                                                    tCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAAc                         <  1:4672833/63‑1 (MQ=255)
                                                      aaCTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAAc                         <  1:3245803/61‑1 (MQ=255)
                                                      aaCTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGa                <  1:459473/70‑1 (MQ=255)
                                                      aaCTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGa                <  1:3233615/70‑1 (MQ=255)
                                                        cTATGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaaaa                                            <  1:2049410/40‑1 (MQ=255)
                                                        cTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaa                             >  1:1261745/1‑55 (MQ=255)
                                                        cTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAaag             <  1:2686660/71‑1 (MQ=255)
                                                         tATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAaga                                          >  1:1978274/1‑41 (MQ=255)
                                                         tATGGCAAGGTAGGTAGTATAGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTCTCTGAaag             >  1:3562681/1‑70 (MQ=255)
                                                          aTGGCAAGGTGGGTAGTATTGGTAAATTCCAGACaaaa                                            >  1:1308384/1‑38 (MQ=255)
                                                          aTGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAAc                         >  1:1419550/1‑57 (MQ=255)
                                                           tGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaa                             <  1:236505/52‑1 (MQ=255)
                                                             gCAAGGTGGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGc         <  1:3905329/70‑1 (MQ=255)
                                                             gCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGcc        <  1:2478986/71‑1 (MQ=255)
                                                             gCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGc         <  1:1621966/70‑1 (MQ=255)
                                                              cAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAaagaag          >  1:4136706/1‑68 (MQ=255)
                                                              cAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGcc        <  1:2712810/70‑1 (MQ=255)
                                                                aggtgggtAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAg      >  1:2162952/1‑70 (MQ=255)
                                                                aggtgggtAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAg      >  1:3219947/1‑70 (MQ=255)
                                                                 ggtgggtaGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAGt     >  1:1589186/1‑70 (MQ=255)
                                                                 ggtgggtaGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAGTa    <  1:2748029/71‑1 (MQ=255)
                                                                  gtaggtagTATTGGTAAATTCCAGACAAAAGAGTTCGATAATga                              <  1:4257264/44‑1 (MQ=255)
                                                                  gtaggtagTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAaag             <  1:3113942/61‑1 (MQ=255)
                                                                  gtaggtagTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAaag             <  1:1012731/61‑1 (MQ=255)
                                                                  gtaggtagTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAGTa    >  1:3647574/1‑70 (MQ=255)
                                                                   taggtagTATTGGTAAATTCCAGACAAAAGAGTTCGATAATgaagaa                          >  1:680207/1‑47 (MQ=255)
                                                                   taggtagTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAaag             <  1:4036019/60‑1 (MQ=255)
                                                                   taggtagTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAGTaaa  >  1:2171084/1‑71 (MQ=255)
                                                                    |                                                                     
GAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACAGTGTCTGAAAGAAGCCAGTAAA  >  minE/1590834‑1590971

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: