Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 868,858 T→G 6.5% V191V (GTT→GTG sohB → predicted inner membrane peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported868,8580TG 6.5% 160.9 / 6.0 62V191V (GTT→GTGsohBpredicted inner membrane peptidase
Reads supporting (aligned to +/- strand):  ref base T (25/33);  new base G (2/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCT  >  Exported/868789‑868926
                                                                     |                                                                    
aaCTGTTACGGTAGCCAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGttt                                                                     >  1:1549437/1‑71 (MQ=255)
  cTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTcc                                                                   <  1:1261991/71‑1 (MQ=255)
     ttACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCg                                                                  <  1:1173116/69‑1 (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:2281860/1‑71 (MQ=255)
           tAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGttt                                                          <  1:2924684/71‑1 (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:3508469/48‑1 (MQ=25)
                          ccAGCGGCGGTTACATGATGGCCTGTGTGGCGGAAACAATTTTGTCCg                                                                  <  1:1728659/48‑1 (MQ=25)
                            aGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCAt                                         <  1:562443/71‑1 (MQ=255)
                            aGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCAt                                         <  1:2645047/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:2927328/71‑1 (MQ=255)
                                            tgGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCaa                         <  1:1073524/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:3550952/1‑70 (MQ=255)
                                                 tgtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGc                     >  1:328421/1‑70 (MQ=255)
                                                  gtgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTggg                                               >  1:1267936/1‑43 (MQ=255)
                                                   tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGccc                   <  1:753912/70‑1 (MQ=255)
                                                   tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGccc                   >  1:488293/1‑70 (MQ=255)
                                                   tgtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCa                  <  1:255261/71‑1 (MQ=255)
                                                    gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                                       <  1:2409160/49‑1 (MQ=255)
                                                    gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                                       <  1:3481199/49‑1 (MQ=255)
                                                    gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                                       <  1:3130057/49‑1 (MQ=255)
                                                    gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                                       <  1:3063692/49‑1 (MQ=255)
                                                    gtgGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTg                                       <  1:3000332/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)
                                                                aTTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCtt     <  1:921644/71‑1 (MQ=255)
                                                                aTTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCtt     <  1:1442797/71‑1 (MQ=255)
                                                                 ttGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAAt                       <  1:2015193/52‑1 (MQ=255)
                                                                 ttGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAAt                       <  1:2113583/52‑1 (MQ=255)
                                                                 ttGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAAt                       <  1:1707392/52‑1 (MQ=255)
                                                                 ttGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAAt                       <  1:785823/52‑1 (MQ=255)
                                                                   gTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAAc          >  1:2109392/1‑63 (MQ=255)
                                                                    tttCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGc       >  1:3200046/1‑65 (MQ=255)
                                                                    tttCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCt  <  1:3378767/70‑1 (MQ=255)
                                                                    tttCCGAACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCt      >  1:2246367/1‑66 (MQ=255)
                                                                     ttCCGCACCTTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACttt             <  1:1106419/58‑1 (MQ=255)
                                                                     ttCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACttt             <  1:26809/58‑1 (MQ=255)
                                                                     ttCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACttt             <  1:499925/58‑1 (MQ=255)
                                                                     |                                                                    
AACTGTTACGGTAGACAAAGTCGCTGCCAGCGGCGGTTACATGATGGCCTGTGTGGCGGACAAAATTGTTTCCGCACCGTTTGCTATTGTGGGTTCCATTGGCGTGGTGGCGCAAATGCCCAACTTTAACCGCTTCCT  >  Exported/868789‑868926

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: