Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,289,208 C→T 17.8% G279G (GGC→GGT amn → AMP nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,289,2080CT17.8% 133.2 / 20.7 56G279G (GGC→GGTamnAMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (19/27);  new base T (6/4);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

AAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCGGATGTCTGGTTGATGA  >  minE/1289138‑1289278
                                                                      |                                                                      
aaaaaCATCAGATGCCAGCATGGTATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGc                                                                        >  1:182990/1‑71 (MQ=255)
aaaaaCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGc                                                                        >  1:447110/1‑71 (MQ=255)
   aaCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTg                                                                     >  1:97760/1‑71 (MQ=255)
   aaCATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTg                                                                     >  1:2424822/1‑71 (MQ=255)
         aGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCg                                                               <  1:3236950/71‑1 (MQ=255)
            tGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGTGTGGGACCGTCa                                                            <  1:1350749/71‑1 (MQ=255)
            tGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGTGTGGGACCGTCa                                                            <  1:186888/71‑1 (MQ=255)
                     ggCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCaaa                                                          >  1:2863304/1‑64 (MQ=255)
                     ggCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTaa                                                    <  1:3464336/70‑1 (MQ=255)
                        aTTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCaaa                                                          >  1:1204020/1‑61 (MQ=255)
                        aTTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAAcc                                                <  1:2545938/71‑1 (MQ=255)
                        aTTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAAcc                                                <  1:1610222/71‑1 (MQ=255)
                          ttAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGtgtg                                                                     <  1:3470915/48‑1 (MQ=255)
                           tAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTg                                                                     >  1:850124/1‑47 (MQ=255)
                           tAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCAt                                              >  1:2241740/1‑70 (MQ=255)
                           tAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATc                                             >  1:1006269/1‑71 (MQ=255)
                             aTTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTg                                           <  1:867634/71‑1 (MQ=255)
                             aTTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGAc                                                                 <  1:1556312/49‑1 (MQ=255)
                             aTTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGAc                                                                 <  1:3945695/49‑1 (MQ=255)
                             aTTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGAc                                                                 <  1:598075/49‑1 (MQ=255)
                                aCCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTGCGa                                        >  1:1186662/1‑71 (MQ=255)
                                aCCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTc                                                             >  1:1406270/1‑50 (MQ=255)
                                aCCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCg                                         <  1:3411921/70‑1 (MQ=255)
                                  cgccgATGGTCAGGGTATTACTCTGGTGAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTGCGatc                                      <  1:3961532/71‑1 (MQ=255)
                                  cgccgATGGTCAGGGTATTACTCTGGTGAATATTGGCGTggg                                                                   <  1:125812/42‑1 (MQ=255)
                                    ccgATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGatcat                                    <  1:471331/71‑1 (MQ=255)
                                      gATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCt                                  <  1:2924282/71‑1 (MQ=255)
                                       aTGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCt                                  <  1:332225/70‑1 (MQ=255)
                                           tCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAg                             <  1:3793893/71‑1 (MQ=255)
                                            cAGGGTATTACTCTGGTGAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGt                            >  1:2498164/1‑71 (MQ=255)
                                             aGGGTATTACTCTGGTGAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTg                           >  1:686818/1‑71 (MQ=255)
                                              gggTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCt                                  <  1:2972279/63‑1 (MQ=255)
                                              gggTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGc                          >  1:2798200/1‑71 (MQ=255)
                                               ggTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGc                          >  1:2303443/1‑70 (MQ=255)
                                                   ttACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGc                                          <  1:749184/50‑1 (MQ=255)
                                                   ttACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTgg                                <  1:1322142/60‑1 (MQ=255)
                                                   ttACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGc                          <  1:1116043/66‑1 (MQ=255)
                                                     aCTCTGGTGAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTGCGatcat                                    >  1:905057/1‑54 (MQ=255)
                                                     aCTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCt                                  >  1:2018615/1‑56 (MQ=255)
                                                     aCTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGccc                   >  1:2397635/1‑71 (MQ=255)
                                                     aCTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGccc                   >  1:2058775/1‑71 (MQ=255)
                                                     aCTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCAACTGCGATCATCTGGCAGTGCTACGccc                   >  1:2095537/1‑71 (MQ=255)
                                                      ctctGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCg                  >  1:1020816/1‑71 (MQ=255)
                                                       tctGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCg                                         <  1:1028300/47‑1 (MQ=255)
                                                       tctGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGatc                                      <  1:2918721/50‑1 (MQ=255)
                                                       tctGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGatc                                      <  1:191806/50‑1 (MQ=255)
                                                       tctGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGatc                                      <  1:344665/50‑1 (MQ=255)
                                                        ctGGTGAATATTGGCGTGGGACCGTCAAATGCTaaaa                                                  >  1:2774925/1‑37 (MQ=255)
                                                        ctGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCg                                         <  1:2444818/46‑1 (MQ=255)
                                                         tGGTGAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTGCGat                                       >  1:1961115/1‑47 (MQ=255)
                                                         tGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCg                  <  1:4105022/68‑1 (MQ=255)
                                                          ggTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGc                          <  1:3867007/59‑1 (MQ=255)
                                                           gTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAg                             <  1:3482283/55‑1 (MQ=255)
                                                            tGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTa                        >  1:255736/1‑59 (MQ=255)
                                                             gAATATTGGTGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCGGATGTCt           >  1:420944/1‑71 (MQ=255)
                                                             gAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCGGATGTCt           >  1:3354228/1‑71 (MQ=255)
                                                              aaTATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCg                  <  1:2839997/63‑1 (MQ=255)
                                                                tatTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGt                            >  1:3912141/1‑51 (MQ=255)
                                                                  ttGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCGGATGTCTGGTtg      <  1:1338004/71‑1 (MQ=255)
                                                                      cGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCGGATGTCTGGTtgatga  <  1:3894311/71‑1 (MQ=255)
                                                                      |                                                                      
AAAAACATCAGATGCCAGCATGGCATTTAATTACCGCCGATGGTCAGGGTATTACTCTGGTGAATATTGGCGTGGGACCGTCAAATGCTAAAACCATCTGCGATCATCTGGCAGTGCTACGCCCGGATGTCTGGTTGATGA  >  minE/1289138‑1289278

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: