Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 821,445 T→C 38.5% S7G (AGC→GGC)  ldrA ← toxic polypeptide, small

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE821,4450TC38.5% 30.3 / 38.7 39S7G (AGC→GGC) ldrAtoxic polypeptide, small
Reads supporting (aligned to +/- strand):  ref base T (13/11);  new base C (7/8);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGAC  >  minE/821395‑821511
                                                  |                                                                  
cGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGt                                                <  1:2078650/71‑1 (MQ=255)
 gCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTc                                               <  1:673433/71‑1 (MQ=34)
     ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCa                                              >  1:2619757/1‑68 (MQ=35)
     ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCa                                              >  1:829238/1‑68 (MQ=35)
     ccaGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATg                                            >  1:1722756/1‑70 (MQ=35)
       aGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGc                                          >  1:2178173/1‑70 (MQ=35)
        gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCa                                         >  1:3702546/1‑70 (MQ=35)
        gCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCa                                         >  1:973467/1‑70 (MQ=35)
          cgACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAAc                                      <  1:2772487/70‑1 (MQ=19)
           gACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACt                                     <  1:3640022/71‑1 (MQ=21)
             cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgc                                   >  1:807236/1‑71 (MQ=34)
             cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTg                                    >  1:2803204/1‑70 (MQ=35)
             cAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgc                                   >  1:591217/1‑71 (MQ=21)
               aTCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGt                                                <  1:933498/56‑1 (MQ=33)
               aTCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgcgc                                 <  1:1447074/71‑1 (MQ=21)
                  gctgctGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGa                               <  1:3317541/70‑1 (MQ=25)
                   cTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGc                             >  1:467586/1‑71 (MQ=34)
                     gCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGt                           >  1:358348/1‑71 (MQ=255)
                      cGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGt                                                >  1:3782639/1‑49 (MQ=255)
                        gTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCat                        <  1:2510512/71‑1 (MQ=34)
                         taataaTTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCata                       <  1:3015779/71‑1 (MQ=34)
                         taataaTTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCata                       <  1:1158666/71‑1 (MQ=34)
                             aatTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATg                                            >  1:3704152/1‑46 (MQ=38)
                             aatTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATg                                            >  1:1257121/1‑46 (MQ=38)
                              atTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGcgcg                                <  1:233042/57‑1 (MQ=35)
                               ttCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgc                                   <  1:1125714/53‑1 (MQ=35)
                               ttCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAgcc                  >  1:1607271/1‑70 (MQ=35)
                                   cGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGTCAAACTGCGCGAGCGTCATATAGCCGCCtt             <  1:3484011/71‑1 (MQ=21)
                                   cGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCtt             <  1:1269293/71‑1 (MQ=34)
                                        ggATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCTCGAGCGTCATATAGCCGCCTtgttgt        <  1:1617364/71‑1 (MQ=255)
                                        ggATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCttgttg         >  1:1801539/1‑70 (MQ=35)
                                         gATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTtgttgt        >  1:2965086/1‑70 (MQ=25)
                                         gATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTtgttgt        >  1:2130256/1‑70 (MQ=25)
                                         gATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTa       <  1:1069279/71‑1 (MQ=21)
                                          aTCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTaa      <  1:383146/71‑1 (MQ=21)
                                            cGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTgcgc                                 >  1:3932904/1‑42 (MQ=255)
                                            cGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCtt             >  1:542941/1‑62 (MQ=25)
                                             ggTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATg    <  1:135627/70‑1 (MQ=25)
                                             ggTGCCGCCAGGTCGTGCCAGAAAGACATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTCAtt    <  1:3483857/70‑2 (MQ=11)
                                               tgctgcCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGac  >  1:1841920/1‑70 (MQ=35)
                                                  |                                                                  
CGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGAC  >  minE/821395‑821511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: