Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 871,512 T→C 10.4% I504T (ATT→ACT)  topA → DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE871,5120TC10.4% 175.8 / 10.6 67I504T (ATT→ACT) topADNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base T (20/40);  new base C (2/5);  total (22/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

GCTGGTTAAAGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTC  >  minE/871444‑871579
                                                                    |                                                                   
gCTGGTTAAAGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCAtt                                                                    >  1:705937/1‑70 (MQ=255)
      tAAAGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATc                                                             >  1:292143/1‑71 (MQ=255)
       aaaGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCACTCAGGATc                                                             <  1:240537/70‑1 (MQ=255)
       aaaGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCAACCTATGCGTCGATCATTTCGACCATTCAGGATc                                                             >  1:2239606/1‑70 (MQ=255)
        aaGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGt                                                           <  1:3539840/71‑1 (MQ=255)
        aaGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGt                                                           <  1:3165828/71‑1 (MQ=255)
            gCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGc                                                        >  1:577191/1‑70 (MQ=255)
            gCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCt                                                       <  1:1287335/71‑1 (MQ=255)
                 aaaaaCGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTg                                                  >  1:778427/1‑71 (MQ=255)
                  aaaaCGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGc                                                 <  1:3086369/71‑1 (MQ=255)
                  aaaaCGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCACTCAGGATCGTGGCTACGTGc                                                 <  1:196884/71‑1 (MQ=255)
                   aaaCGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCg                                                <  1:3395628/71‑1 (MQ=255)
                    aaCGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGa                                               <  1:2458617/71‑1 (MQ=255)
                     tcgcgGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCg                                                <  1:2512756/68‑1 (MQ=255)
                     aCGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAg                                              >  1:3395682/1‑71 (MQ=255)
                     aCGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAg                                              >  1:2937217/1‑71 (MQ=255)
                      cgcgGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAg                                              >  1:3659015/1‑70 (MQ=255)
                          gTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCACTCAGGATCg                                                            <  1:3749629/52‑1 (MQ=255)
                            aTCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCACTCAGGATCGTg                                                          >  1:1325819/1‑52 (MQ=255)
                            aTCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCACTCAGGATCGTg                                                          >  1:1321122/1‑52 (MQ=255)
                             tCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGaaaa                                       <  1:1881674/70‑1 (MQ=255)
                             tCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAAt                                      >  1:2555946/1‑71 (MQ=255)
                               ggTCGTCCGTTTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAAtcg                                    >  1:1526481/1‑71 (MQ=255)
                               ggTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAAtcg                                    >  1:347542/1‑71 (MQ=255)
                               ggTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAAtcg                                    <  1:3635669/71‑1 (MQ=255)
                                 tcgtcCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAAtcgtc                                  <  1:1918246/71‑1 (MQ=255)
                                 tcgtcCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAAtcgtc                                  <  1:1396576/71‑1 (MQ=255)
                                   gtccgtcTACCTATGCGTCGATCATTTCGACCACTCAGGATCGTGGCTACGTGCGAGTAGAAAATcgtcgt                                <  1:1622344/71‑1 (MQ=255)
                                    tccgtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGt                                                   <  1:1359483/51‑1 (MQ=255)
                                    tccgtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGa                                          <  1:3654075/60‑1 (MQ=255)
                                    tccgtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATcgtcgt                                <  1:2103126/70‑1 (MQ=255)
                                      cgtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGTCTACGTGCGAGTAGAAAATCGTCGttt                              <  1:1409827/70‑1 (MQ=255)
                                       gtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGc                                                 <  1:2835657/50‑1 (MQ=255)
                                       gtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATcgtcgt                                >  1:1950373/1‑67 (MQ=255)
                                       gtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCt                            <  1:658056/71‑1 (MQ=255)
                                       gtcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCt                            <  1:3626620/71‑1 (MQ=255)
                                        tcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGc                                                 <  1:791843/49‑1 (MQ=255)
                                        tcTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTa                           >  1:662853/1‑71 (MQ=255)
                                        tcTACCTATGCGTCGATCATTTCGACCACTCAGGATCg                                                            <  1:1269626/38‑1 (MQ=255)
                                          tACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATg                         >  1:1376223/1‑71 (MQ=255)
                                           aCCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGc                        >  1:1848318/1‑71 (MQ=255)
                                           aCCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGc                        >  1:1850424/1‑71 (MQ=255)
                                            ccTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATcgtcgt                                <  1:2683224/62‑1 (MQ=255)
                                            ccTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTAt                          <  1:1526517/68‑1 (MQ=255)
                                            ccTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATg                         <  1:646660/69‑1 (MQ=255)
                                            ccTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATg                         <  1:2637198/69‑1 (MQ=255)
                                                 gCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAAtcgtc                                  >  1:409145/1‑55 (MQ=255)
                                                  cGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGaaaa                  <  1:907079/70‑1 (MQ=255)
                                                   gTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGa                                          <  1:1923900/45‑1 (MQ=255)
                                                   gTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGaaaaaa                <  1:405341/71‑1 (MQ=255)
                                                    tCGATCATTTCGACCATTCAGGATCGTGGCTACGt                                                   <  1:487850/35‑1 (MQ=255)
                                                    tCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAAt               <  1:2601900/71‑1 (MQ=255)
                                                    tCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAAt               <  1:3795223/71‑1 (MQ=255)
                                                    tCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAAt               <  1:106694/71‑1 (MQ=255)
                                                       aTCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGaaaaa                 <  1:1123033/66‑1 (MQ=255)
                                                        tCATTTCGACCATTCAGGATCGTGGCTACGTGCGAg                                              <  1:3939040/36‑1 (MQ=255)
                                                         cATTTCGACCATTCAGGATCGTGGCTACGTGCGAg                                              >  1:2851118/1‑35 (MQ=255)
                                                          aTTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATcgtcgt                                <  1:2708568/48‑1 (MQ=255)
                                                          aTTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATcgtcgt                                <  1:1280791/48‑1 (MQ=255)
                                                          aTTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCg          <  1:3709316/70‑1 (MQ=255)
                                                           tttCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGa         >  1:2594939/1‑70 (MQ=255)
                                                           tttCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGa         >  1:2453245/1‑70 (MQ=255)
                                                           tttCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGa         >  1:92984/1‑70 (MQ=255)
                                                               gACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCg    <  1:3659459/71‑1 (MQ=255)
                                                               gACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCCGAAAAAATGGGCGAAATCg    <  1:1199626/71‑1 (MQ=255)
                                                               gACCATTCAGGATCGTGGCTACGTACGTGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCg    <  1:2498344/71‑1 (MQ=255)
                                                                 ccATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTc  <  1:3159026/71‑1 (MQ=255)
                                                                 ccATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTc  <  1:786180/71‑1 (MQ=255)
                                                                    ttCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGaaaa                  <  1:1548237/52‑1 (MQ=255)
                                                                    ttCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGaaaa                  <  1:1348925/52‑1 (MQ=255)
                                                                    |                                                                   
GCTGGTTAAAGAGCTGGAAAAACGCGGTATCGGTCGTCCGTCTACCTATGCGTCGATCATTTCGACCATTCAGGATCGTGGCTACGTGCGAGTAGAAAATCGTCGTTTCTATGCGGAAAAAATGGGCGAAATCGTC  >  minE/871444‑871579

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: