Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 868,851 A→C 7.1% K189T (AAA→ACA)  sohB → predicted inner membrane peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported868,8510AC 7.1% 145.1 / 6.2 56K189T (AAA→ACA) sohBpredicted inner membrane peptidase
Reads supporting (aligned to +/- strand):  ref base A (28/24);  new base C (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

CCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTT  >  Exported/868784‑868915
                                                                   |                                                                
ccTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACaaaa                                                                <  1:773491/70‑1 (MQ=255)
 cTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAAtt                                                              >  1:944664/1‑71 (MQ=255)
     aaCTGTTACGGTAGCCAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGttt                                                          >  1:1549437/1‑71 (MQ=255)
       cTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTcc                                                        <  1:1261991/71‑1 (MQ=255)
          ttACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCg                                                       <  1:1173116/69‑1 (MQ=255)
            aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACaa                                                                  >  1:3440219/1‑56 (MQ=255)
            aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAAt                                                               >  1:644012/1‑59 (MQ=255)
            aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAAt                                                               >  1:2420967/1‑59 (MQ=255)
            aCGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTATGTGGCGGACaa                                                                  >  1:1075562/1‑56 (MQ=255)
             cGGTAGACAAAGTCGCTGCCAGCGGCGGTTATATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCg                                                  >  1:1018428/1‑71 (MQ=255)
             cGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCg                                                  >  1:149972/1‑71 (MQ=255)
               gTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGtt                                                <  1:3323325/71‑1 (MQ=255)
                tAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGttt                                               <  1:2924684/71‑1 (MQ=255)
                tAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGttt                                               >  1:2281860/1‑71 (MQ=255)
                tAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGttt                                               >  1:2203802/1‑71 (MQ=255)
                       aGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtg                                        >  1:3525435/1‑71 (MQ=255)
                               ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTc                                 >  1:3316596/1‑70 (MQ=255)
                               ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTc                                 >  1:492027/1‑70 (MQ=255)
                               ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTc                                 >  1:928853/1‑70 (MQ=255)
                               ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGAAACAATTTTGTCCg                                                       <  1:1728659/48‑1 (MQ=25)
                               ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGAAACAATTTTGTCCg                                                       <  1:3508469/48‑1 (MQ=25)
                                 aGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCAt                              <  1:2645047/71‑1 (MQ=255)
                                 aGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCAt                              <  1:562443/71‑1 (MQ=255)
                                  gcggcgGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtg                                        <  1:2329611/60‑1 (MQ=255)
                                     gcggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGtt                                  <  1:2093305/63‑1 (MQ=255)
                                      cggTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCg                         <  1:2392541/71‑1 (MQ=255)
                                         ttACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtgg                      >  1:2977573/1‑71 (MQ=255)
                                          tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCCCCGTTTGCTATtgtg                                      >  1:26578/1‑54 (MQ=255)
                                          tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg                                      >  1:2386917/1‑54 (MQ=255)
                                          tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg                                      >  1:2110890/1‑54 (MQ=255)
                                          tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg                                      >  1:1827127/1‑54 (MQ=255)
                                          tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATtgtg                                      >  1:1294628/1‑54 (MQ=255)
                                          tACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTgg                           >  1:3551744/1‑65 (MQ=255)
                                           aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtggtg                    >  1:2495669/1‑71 (MQ=255)
                                           aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtggtg                    >  1:1736901/1‑71 (MQ=255)
                                           aCATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCgtggt                     <  1:2234813/70‑1 (MQ=255)
                                                 tgGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa              <  1:1073524/71‑1 (MQ=255)
                                                 tgGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa              <  1:2927328/71‑1 (MQ=255)
                                                 tgGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCa               >  1:1263836/1‑70 (MQ=255)
                                                  ggCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaaa             >  1:1160323/1‑71 (MQ=255)
                                                  ggCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa              <  1:1078019/70‑1 (MQ=255)
                                                      tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGc          >  1:328421/1‑70 (MQ=255)
                                                      tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGc          >  1:3550952/1‑70 (MQ=255)
                                                       gtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTggg                                    >  1:1267936/1‑43 (MQ=255)
                                                        tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGccc        >  1:488293/1‑70 (MQ=255)
                                                        tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGccc        <  1:753912/70‑1 (MQ=255)
                                                        tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCa       <  1:255261/71‑1 (MQ=255)
                                                         gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                            <  1:3481199/49‑1 (MQ=255)
                                                         gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                            <  1:3000332/49‑1 (MQ=255)
                                                         gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                            <  1:2409160/49‑1 (MQ=255)
                                                         gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                            <  1:3130057/49‑1 (MQ=255)
                                                         gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                            <  1:3063692/49‑1 (MQ=255)
                                                         gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                            <  1:3059359/49‑1 (MQ=255)
                                                           tgCGGAAACAATTTTGTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACt    >  1:133263/2‑71 (MQ=38)
                                                           tgCGGAAACAATTTTGTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACt    >  1:2686665/2‑71 (MQ=38)
                                                            gCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACtt   <  1:971128/71‑1 (MQ=255)
                                                             cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGgcgc                <  1:955983/57‑1 (MQ=255)
                                                             cGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACttt  <  1:2163775/71‑1 (MQ=255)
                                                               gACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa              >  1:1108831/1‑57 (MQ=255)
                                                                  aaaaTTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa              <  1:1191850/54‑1 (MQ=255)
                                                                  aaaaTTGTTTCCGCACCGTTTGCTATTGTGGGGTCCATTGGCGTGGTGGCGCaa              <  1:397670/54‑1 (MQ=255)
                                                                   aaaTTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa              <  1:992382/53‑1 (MQ=255)
                                                                   |                                                                
CCTTTAACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTT  >  Exported/868784‑868915

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: