Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 457,813 C→T 33.3% Q41* (CAG→TAG)  ybgL → predicted lactam utilization protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE457,8130CT33.3% 68.4 / 49.8 57Q41* (CAG→TAG) ybgLpredicted lactam utilization protein
Reads supporting (aligned to +/- strand):  ref base C (22/16);  new base T (9/10);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AACGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGAC  >  minE/457752‑457878
                                                             |                                                                 
aaCGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGc                                                                   >  1:3835427/1‑62 (MQ=255)
aaCGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTtgc                                                           >  1:3272348/1‑70 (MQ=255)
       gTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGcgtgcgtg                                                    <  1:2213861/70‑1 (MQ=255)
          tCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAg                                                 >  1:3983823/1‑70 (MQ=255)
             tctGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTtg                                                            <  1:1598323/56‑1 (MQ=255)
               tGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGGGTGAAGCAATa                                            >  1:4112101/1‑70 (MQ=255)
               tGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATa                                            >  1:3934280/1‑70 (MQ=255)
                 ccAATATTGCCTGTGGATTTCATGCAGGCTATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATa                                            <  1:2061929/68‑1 (MQ=255)
                  cAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAg                                                                 <  1:1119477/46‑1 (MQ=255)
                   aaTATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAg                                                                 >  1:1786738/1‑45 (MQ=255)
                    atatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAAt                                      <  1:3014271/71‑1 (MQ=255)
                     tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCtt                                                             <  1:2508663/47‑1 (MQ=255)
                     tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa                                       <  1:1820242/69‑1 (MQ=255)
                     tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaa                                          <  1:3835156/66‑1 (MQ=255)
                     tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaa                                          <  1:3839800/66‑1 (MQ=255)
                     tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaa                                          <  1:3540747/66‑1 (MQ=255)
                     tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATg                                     >  1:3016769/1‑71 (MQ=255)
                     tatTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATg                                     >  1:214500/1‑71 (MQ=255)
                      atTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTtgcg                                                          >  1:3723583/1‑49 (MQ=255)
                      atTGACTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGc                                                >  1:2681523/1‑59 (MQ=255)
                       ttGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGgt                                   >  1:2005888/1‑71 (MQ=255)
                       ttGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCtt                                                             >  1:3367990/1‑45 (MQ=255)
                           cTGTGGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATaaaa                                         <  1:2064656/61‑1 (MQ=255)
                           cTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGAGTGAAGCAATAAAAAATGGTGTcg                               >  1:3955938/1‑71 (MQ=255)
                            tgtgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaa                                         <  1:4056036/60‑1 (MQ=255)
                              tgGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgcg                             <  1:1724838/70‑1 (MQ=255)
                              tgGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgcg                             <  1:2541379/70‑1 (MQ=255)
                              tgGATTTCATGCAGGCGATGCGCAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgc                              <  1:1854743/69‑1 (MQ=255)
                              tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGcgtgcgtg                                                    <  1:3044354/47‑1 (MQ=255)
                              tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTtgc                              <  1:1780446/69‑1 (MQ=255)
                              tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgc                              <  1:3212319/69‑1 (MQ=255)
                              tgGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGa                            >  1:1371059/1‑71 (MQ=255)
                                gATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGAtt                          >  1:2041976/1‑71 (MQ=255)
                                gATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGAtt                          >  1:1809768/1‑71 (MQ=255)
                                 aTTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaa                                           <  1:44830/53‑1 (MQ=255)
                                   ttCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa                                       >  1:2675740/1‑55 (MQ=255)
                                   ttCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa                                       >  1:2497920/1‑55 (MQ=255)
                                         cAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCAc                 >  1:987112/1‑71 (MQ=255)
                                            gCGATGCGCAAATCATGTAGGCTTGCATGCGTGAAGCAATAAAAAATg                                     <  1:1945943/48‑1 (MQ=39)
                                             cGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCAccc               >  1:333278/1‑69 (MQ=255)
                                              gATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa                                       >  1:2525955/1‑44 (MQ=255)
                                                 gcgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCTCTCACCCGAGtt          >  1:1979731/1‑70 (MQ=255)
                                                 gcgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGttt         >  1:555552/1‑71 (MQ=255)
                                                 gcgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGtt          >  1:617474/1‑70 (MQ=255)
                                                  cgcAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCCCCCGAGtttt        >  1:427799/1‑71 (MQ=255)
                                                  cgcAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGttt         <  1:2792444/70‑1 (MQ=255)
                                                    cAAATCATGTAGGCTTGCGTGCGTGAAGCAATaaaaa                                        >  1:3270021/1‑37 (MQ=255)
                                                    cAAATCATGTAGGCTTGCGTGCGTGAAGCAATaaaaa                                        >  1:935573/1‑37 (MQ=255)
                                                    cAAATCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTcgc                              <  1:3222428/47‑1 (MQ=255)
                                                    cAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa                                       <  1:3273294/38‑1 (MQ=255)
                                                    cAAATCATGCAGGCTTGCGTGCGTGAAGCAATaaaaaa                                       <  1:806176/38‑1 (MQ=255)
                                                    cAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTc       <  1:104190/70‑1 (MQ=255)
                                                     aaaTCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCa                  >  1:1088487/1‑58 (MQ=255)
                                                     aaaTCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCa                  >  1:2678156/1‑58 (MQ=255)
                                                        tCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTccc     >  1:1578746/1‑68 (MQ=255)
                                                        tCATGTAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGAc  <  1:1764423/71‑1 (MQ=255)
                                                        tCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGa                            >  1:1249778/1‑45 (MQ=255)
                                                        tCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGAt                           >  1:1597783/1‑46 (MQ=255)
                                                             |                                                                 
AACGCTGGTTTCCTCTGCCAATATTGCCTGTGGATTTCATGCAGGCGATGCGCAAATCATGCAGGCTTGCGTGCGTGAAGCAATAAAAAATGGTGTCGCGATTGGCGCTCACCCGAGTTTTCCCGAC  >  minE/457752‑457878

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: