Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,102,851 G→A 8.6% W143* (TGG→TAG)  yhaV → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,102,8510GA 8.6% 178.6 / 6.7 58W143* (TGG→TAG) yhaVconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (3/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

ACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATT  >  minE/2102782‑2102917
                                                                     |                                                                  
aCGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTgg                                                                   >  1:2423373/1‑71 (MQ=255)
    tAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaa                                                                <  1:601298/70‑1 (MQ=255)
    tAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaa                                                                <  1:1553830/70‑1 (MQ=255)
    tAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaa                                                                <  1:2758825/70‑1 (MQ=255)
      aaaaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAc                                                              >  1:204795/1‑70 (MQ=255)
        aaaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCt                                                           >  1:252492/1‑71 (MQ=255)
         aaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaa                                                                <  1:542023/65‑1 (MQ=255)
         aaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaa                                                                <  1:1883186/65‑1 (MQ=255)
         aaaCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaa                                                                <  1:2995421/65‑1 (MQ=255)
              gATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCg                                                     <  1:1965871/71‑1 (MQ=255)
                      taCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGaaa                                                >  1:1557748/1‑68 (MQ=255)
                       aCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaaa                                                               >  1:1155165/1‑52 (MQ=255)
                       aCCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaaa                                                               >  1:2448683/1‑52 (MQ=255)
                           tATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCt                                                           <  1:566815/52‑1 (MQ=255)
                           tATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCt                                                           <  1:215711/52‑1 (MQ=255)
                             ttCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaaa                                                               <  1:1331427/46‑1 (MQ=255)
                             ttCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaaa                                                               <  1:880503/46‑1 (MQ=255)
                             ttCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAc                                      >  1:1451344/1‑71 (MQ=255)
                             ttCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAc                                      >  1:2378928/1‑71 (MQ=255)
                             ttCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAc                                      >  1:1932497/1‑71 (MQ=255)
                                aGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAccc                                    <  1:1069744/70‑1 (MQ=255)
                                aGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAccc                                    <  1:678934/70‑1 (MQ=255)
                                 gCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTAGGAAACCCTCACCCGAGAAACag                                             <  1:772152/60‑1 (MQ=255)
                                  cAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCAtt                                 <  1:528949/71‑1 (MQ=255)
                                   aaaaTGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCAtt                                 <  1:2212097/70‑1 (MQ=255)
                                   aaaaTGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCAt                                  >  1:1489631/1‑69 (MQ=255)
                                   aaaaTGTTAAAAAGAGGACATCCTCCTGCCGACTAGGAAACCCTCACCCGAGAAACagaa                                           >  1:481074/1‑60 (MQ=255)
                                   aaaaTGTTAAAAAGAGGACATCCTCCTGCCGACTAGGAAACCCTCACCCGAGAAACagaa                                           >  1:826017/1‑60 (MQ=255)
                                   aaaaTGTTAAAAAGAGGACATCCTCCTGCCGACTAGGAAACCCTCACCCGAGAAACagaa                                           >  1:2743748/1‑60 (MQ=255)
                                    aaaTGTTAAAAAGAGGACATCCTCCTGCCGACTgg                                                                   >  1:73491/1‑35 (MQ=255)
                                     aaTGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGa                               <  1:2726665/70‑1 (MQ=255)
                                      aTGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATg                             >  1:1496259/1‑71 (MQ=255)
                                       tGTTAAAAAGAGGACATCCTCCTGCCGACTAGGAAACCCTCACCCGAGAAACag                                             <  1:2891295/54‑1 (MQ=255)
                                         ttAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGgtg                          >  1:2699562/1‑71 (MQ=255)
                                         ttAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGgtg                          <  1:798054/71‑1 (MQ=255)
                                         ttAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGgtg                          >  1:2234209/1‑71 (MQ=255)
                                         ttAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATCgtg                          >  1:2703592/1‑71 (MQ=255)
                                          tAAAAAGAGGACATCCTCCTGCCGACTGGGAACCCTCACCCGAGAAACaga                                            >  1:799881/1‑51 (MQ=255)
                                           aaaaaGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGctgt                         <  1:1065803/70‑1 (MQ=255)
                                            aaaaGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCAt                                  <  1:2331928/60‑1 (MQ=255)
                                            aaaaGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGtt                        <  1:917087/70‑1 (MQ=255)
                                              aaGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTcac                     <  1:3408026/71‑1 (MQ=255)
                                              aaGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTcac                     <  1:1917132/71‑1 (MQ=255)
                                                   gACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGc                <  1:3032344/71‑1 (MQ=255)
                                                   gACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGc                <  1:2531292/71‑1 (MQ=255)
                                                     cATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGccc              >  1:2255891/1‑71 (MQ=255)
                                                     cATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGccc              >  1:1796030/1‑71 (MQ=255)
                                                      aTCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACag                                             <  1:473187/39‑1 (MQ=255)
                                                      aTCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccg             >  1:2957786/1‑71 (MQ=255)
                                                       tcctccTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccgc            >  1:3024110/1‑71 (MQ=255)
                                                       tcctccTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccgc            >  1:539575/1‑71 (MQ=255)
                                                       tcctccTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccgc            >  1:1137785/1‑71 (MQ=255)
                                                         ctcctGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccgccg          <  1:1156734/71‑1 (MQ=255)
                                                          tcctGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCccgccg          >  1:611167/1‑70 (MQ=255)
                                                             tGCCGACTGGGAAACCCTCACCCGAGAAACagaag                                          >  1:1908798/1‑35 (MQ=255)
                                                             tGCCGACTGGGAAACCCTCACCCGAGAAACagaag                                          >  1:18188/1‑35 (MQ=255)
                                                               ccGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACg      <  1:654688/69‑1 (MQ=255)
                                                               ccGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACg      <  1:3114121/69‑1 (MQ=255)
                                                                cGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCAt   <  1:1090713/71‑1 (MQ=255)
                                                                 gACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCAtt  >  1:196534/1‑71 (MQ=255)
                                                                     |                                                                  
ACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATTGATGGTGTTCACAATGCCCGCCGAACGCATT  >  minE/2102782‑2102917

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: