Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,226,993 G→A 5.8% L22L (CTG→CTA yhdY → predicted amino‑acid transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,226,9930GA 5.8% 213.5 / 2.4 69L22L (CTG→CTAyhdYpredicted amino‑acid transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AGGACATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGC  >  minE/2226924‑2227060
                                                                     |                                                                   
aGGACATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCctgc                                                                    >  1:2140974/1‑71 (MQ=255)
aGGACATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCctg                                                                     >  1:550492/1‑70 (MQ=255)
       aaCTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTAt                                                             >  1:3403381/1‑71 (MQ=255)
        aCTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTAtt                                                            >  1:1133014/1‑71 (MQ=255)
        aCTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATATGTTCTCCAGCTGGAGCAATAGCCTGCTGACTAtt                                                            <  1:252504/71‑1 (MQ=255)
        aCTTAAGACGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTAtt                                                            >  1:84233/1‑71 (MQ=255)
            aaGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCt                                                        >  1:4093639/1‑71 (MQ=255)
            aaGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCt                                                        >  1:1834933/1‑71 (MQ=255)
             aGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCt                                                        <  1:1022956/70‑1 (MQ=255)
               ccGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCa                                                     <  1:4558359/71‑1 (MQ=255)
                cgcgcgATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGc                                                      >  1:1713963/1‑69 (MQ=255)
                cgcgcgATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCat                                                    <  1:1367537/71‑1 (MQ=255)
                   gcgATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACt                                                               <  1:2634916/57‑1 (MQ=255)
                   gcgATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATg                                                 >  1:737229/1‑71 (MQ=255)
                    cgATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATgg                                                >  1:3308894/1‑71 (MQ=255)
                       tGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCctgctg                                                                  >  1:4222179/1‑50 (MQ=255)
                       tGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGtt                                              >  1:87167/1‑70 (MQ=255)
                        ggTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGTTGACTa                                                              <  1:602557/53‑1 (MQ=255)
                        ggTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCa                                                     >  1:287715/1‑62 (MQ=255)
                        gctgtgGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGa                                            >  1:567663/3‑71 (MQ=255)
                            tgGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGa                                            <  1:3371731/67‑1 (MQ=255)
                             gggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTAtt                                                            >  1:4186508/1‑50 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:2409310/71‑1 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:2993524/71‑1 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:289379/71‑1 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:4505528/71‑1 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:1218425/71‑1 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTgg                                       >  1:2329580/1‑70 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTgg                                       >  1:3450285/1‑70 (MQ=255)
                              ggTGCGAAAAAATCTGTTCTCCAGCCGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:3160029/71‑1 (MQ=255)
                                tGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGAtgt                                         <  1:401344/66‑1 (MQ=255)
                                 gCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGGTGCata                                                   >  1:3038875/1‑55 (MQ=255)
                                 gCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAAc                                   >  1:3860027/1‑71 (MQ=255)
                                  cGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACt                                  <  1:4674784/71‑1 (MQ=255)
                                   gAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTg                                 <  1:3889991/71‑1 (MQ=255)
                                   gAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTg                                 <  1:961718/71‑1 (MQ=255)
                                    aaaaaaTCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTg                                 <  1:2587297/70‑1 (MQ=255)
                                    aaaaaaTCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGa                                >  1:1731171/1‑71 (MQ=255)
                                    aaaaaaTCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTATGTggg                                      <  1:401257/64‑1 (MQ=255)
                                     aaaaaTCTGTTCTCCAGCTGGAGCAATAGCctact                                                                   >  1:1104691/1‑35 (MQ=255)
                                     aaaaaTCTGTTCTCCAGCTGGAGCAATAGCCTGCTGAc                                                                >  1:1875146/1‑38 (MQ=255)
                                      aaaaTCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGc                                                      <  1:1724359/47‑1 (MQ=255)
                                      aaaaTCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATc                              >  1:654310/1‑71 (MQ=255)
                                           cTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGCTCCCAc                          >  1:922906/1‑70 (MQ=255)
                                            tGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGa                                            <  1:3021737/51‑1 (MQ=255)
                                            tGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCg                        >  1:177850/1‑71 (MQ=255)
                                              ttCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCg                        <  1:1967809/69‑1 (MQ=255)
                                                 tcCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCt                   >  1:3503139/1‑71 (MQ=255)
                                                   cAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGa                                >  1:2054257/1‑56 (MQ=255)
                                                   cAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGGTg                  <  1:1082057/70‑1 (MQ=255)
                                                       tGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:32138/46‑1 (MQ=255)
                                                        ggAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTgg            >  1:4496224/1‑71 (MQ=255)
                                                         gAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTgg            >  1:2464308/1‑70 (MQ=255)
                                                         gAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTg             <  1:3758459/69‑1 (MQ=255)
                                                         gAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTg             <  1:3131491/69‑1 (MQ=255)
                                                          aGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGAtgt                                         >  1:1894021/1‑40 (MQ=255)
                                                           gCAATAGCCTGCTGACTATTGGCTGCATATGGTTGAt                                           >  1:4050423/1‑37 (MQ=255)
                                                           gCAATAGCCTACTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCg                        >  1:3591917/1‑56 (MQ=255)
                                                           gCAATAGCCTACTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCa         >  1:3388920/1‑71 (MQ=255)
                                                             aaTAGCCTACTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTg                     <  1:1806962/57‑1 (MQ=255)
                                                             aaTAGCCTACTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGa                 <  1:71722/61‑1 (MQ=255)
                                                              aTAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGTAACTGATCCCACCGTTGCTGAAc               >  1:4474547/1‑62 (MQ=255)
                                                              aTAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGt                       >  1:1149118/1‑54 (MQ=255)
                                                              aTAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGt                       <  1:743581/54‑1 (MQ=255)
                                                               tAGCCTGCTGACTATTGGCTGCATATGGTTGATGTggg                                      <  1:2847924/38‑1 (MQ=255)
                                                               tAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCAtttt     >  1:203853/1‑71 (MQ=255)
                                                                aGCCTGCTGACTATTGGCTGCATATGGTTGAtgtg                                        >  1:1434282/1‑35 (MQ=255)
                                                                 gCCTGCTGACTATTGGCTGCATATGGTTGATGTgg                                       >  1:1809178/1‑35 (MQ=255)
                                                                 gCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGaa                                    >  1:189925/1‑38 (MQ=255)
                                                                 gCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTg                     <  1:2235831/53‑1 (MQ=255)
                                                                 gCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCAttttt    <  1:944179/70‑1 (MQ=255)
                                                                  ccTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGc  <  1:527763/71‑1 (MQ=255)
                                                                  ccTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGc  <  1:4321730/71‑1 (MQ=255)
                                                                   ctgctgACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGc  <  1:2239869/70‑1 (MQ=255)
                                                                    tgctgACTATTGGCTGCATATGGTTGATGTGGGAACTGAt                               <  1:3387361/40‑1 (MQ=255)
                                                                     gctgACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTg   <  1:4583243/67‑1 (MQ=255)
                                                                     |                                                                   
AGGACATAACTTAAGCCGCGCGATGGTGTGGGTGCGAAAAAATCTGTTCTCCAGCTGGAGCAATAGCCTGCTGACTATTGGCTGCATATGGTTGATGTGGGAACTGATCCCACCGTTGCTGAACTGGGCATTTTTGC  >  minE/2226924‑2227060

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: