Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,467,238 G→T 16.6% L234L (CTG→CTT atpA → F1 sector of membrane‑bound ATP synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,467,2380GT16.6% 98.3 / 5.3 36L234L (CTG→CTTatpAF1 sector of membrane‑bound ATP synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (23/7);  new base T (3/3);  total (26/10)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GGAAGAGCACGGCGCACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATGCGCTGA  >  minE/2467172‑2467308
                                                                  |                                                                      
ggAAGAGCACGGCGCACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCaa                                                                     >  1:1270926/1‑70 (MQ=255)
 gAAGAGCACGGCGCACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTTCACTTCAAt                                                                    <  1:1761504/70‑1 (MQ=255)
 gAAGAGCACGGCGCACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTTCAATa                                                                   <  1:518883/71‑1 (MQ=255)
           gcgcACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTTCAATACCTGGCACCg                                                         >  1:1111233/1‑71 (MQ=255)
           gcgcACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCg                                                         >  1:1982997/1‑71 (MQ=255)
           gcgcACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCg                                                         <  1:37255/71‑1 (MQ=255)
                        cacCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTgcg                                            >  1:169852/1‑71 (MQ=255)
                        cacCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTgcg                                            >  1:1167503/1‑71 (MQ=255)
                        cacCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTgcg                                            >  1:1758674/1‑71 (MQ=255)
                         acCATCGTTGTGGTAGCAACCGCGTCTTAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTgcgc                                           <  1:2142054/71‑1 (MQ=255)
                          ccATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTTCAATACCTGGCACCGTATGCCGGTTGCGCa                                          <  1:1820151/71‑1 (MQ=255)
                          ccATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTgcgc                                           <  1:1329941/70‑1 (MQ=255)
                            aTCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTTCAATACCTGGCACCGTATGCCGGTTGCGCAAt                                        >  1:193760/1‑71 (MQ=255)
                            aTCGTTGTGGTAGCAACCGCGTCAGAATCCGCTGCACTGCAATAc                                                                  >  1:1113052/1‑45 (MQ=255)
                             tCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTg                                              >  1:1242036/1‑64 (MQ=255)
                             tCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCCTATGCCGGTTGCGCAAtt                                       >  1:1790039/1‑70 (MQ=255)
                              cGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATgg                                      >  1:1010536/1‑71 (MQ=255)
                               gTTGTGGTAGCAACCGCGTCTGAATCCTCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATggg                                     >  1:992936/1‑71 (MQ=255)
                                ttGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTa                                                       >  1:390201/1‑52 (MQ=255)
                                  gtgGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCg                                   >  1:1702667/1‑70 (MQ=255)
                                   tgGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGaa                                 <  1:13102/71‑1 (MQ=255)
                                     gTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTg                                                               >  1:1079173/1‑39 (MQ=255)
                                          aaCCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGtt                                               >  1:297523/1‑50 (MQ=255)
                                          aaCCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTcc                          >  1:1350908/1‑71 (MQ=255)
                                          aaCCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTcc                          <  1:1875194/71‑1 (MQ=255)
                                            ccGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCTGTTGCGCAATGGGCGAATACTTCCGt                        <  1:696996/71‑1 (MQ=255)
                                            ccGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGt                        <  1:1969641/71‑1 (MQ=255)
                                              gcgTCTGAATCCGCTGCACTTCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGa                      >  1:810450/1‑71 (MQ=255)
                                              gcgTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGa                      >  1:64501/1‑71 (MQ=255)
                                               cgTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGAc                     >  1:1699021/1‑71 (MQ=255)
                                                 tCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACcg                   >  1:699777/1‑71 (MQ=255)
                                                   tGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACcgcg                 >  1:1652749/1‑71 (MQ=255)
                                                   tGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACcgcg                 >  1:1707727/1‑71 (MQ=255)
                                                    gAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACcgcg                 >  1:2201594/1‑70 (MQ=255)
                                                     aaTCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGt               >  1:1727806/1‑71 (MQ=255)
                                                          gCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGt                        >  1:14291/1‑57 (MQ=255)
                                                           ctgcactgcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGAt         <  1:1654620/71‑1 (MQ=255)
                                                            tgcactgcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGaa            <  1:170733/67‑1 (MQ=255)
                                                            tgcactgcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGaa            <  1:1243379/67‑1 (MQ=255)
                                                               actgcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATgcgc     >  1:159225/1‑71 (MQ=255)
                                                                ctgcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATGCGCt    >  1:1777966/1‑71 (MQ=255)
                                                                 tgcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGa          <  1:164052/64‑1 (MQ=255)
                                                                 tgcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGa          <  1:121828/64‑1 (MQ=255)
                                                                  gcaATACCTGGCACCGTATGCCGTTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATGCGCTGa  >  1:2151513/1‑71 (MQ=255)
                                                                  gcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATGCGCTGa  >  1:150039/1‑71 (MQ=255)
                                                                  gcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATGCGCTGa  >  1:1400615/1‑71 (MQ=255)
                                                                  gcaATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATGCGCTGa  >  1:1044182/1‑71 (MQ=255)
                                                                  |                                                                      
GGAAGAGCACGGCGCACTGGCTAACACCATCGTTGTGGTAGCAACCGCGTCTGAATCCGCTGCACTGCAATACCTGGCACCGTATGCCGGTTGCGCAATGGGCGAATACTTCCGTGACCGCGGTGAAGATGCGCTGA  >  minE/2467172‑2467308

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: