Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,625,149 G→A 14.3% R16H (CGT→CAT)  ftsX → predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,625,1490GA14.3% 149.7 / 15.4 56R16H (CGT→CAT) ftsXpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base A (2/6);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CACTTGCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGTAACGCACCAAAACG  >  Exported/2625081‑2625212
                                                                    |                                                               
cACTTGCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGt                                                                >  1:978401/1‑70 (MQ=255)
 aCTTGCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCATCt                                                              >  1:4343981/1‑71 (MQ=255)
     gCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTg                                                            <  1:3183533/69‑1 (MQ=255)
     gCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTg                                                            <  1:4594076/69‑1 (MQ=255)
     gCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTg                                                            <  1:4488780/69‑1 (MQ=255)
     gCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTg                                                            <  1:4339193/69‑1 (MQ=255)
     gCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTg                                                            <  1:2537370/69‑1 (MQ=255)
     gCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGAt                                                          >  1:1282369/1‑71 (MQ=255)
     gCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGAt                                                          >  1:1332022/1‑71 (MQ=255)
       aTGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATc                                                         >  1:4346075/1‑70 (MQ=255)
          gAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCt                                                              >  1:650619/1‑62 (MQ=255)
            ggCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCATCTTGATCGCTTcc                                                   >  1:1659923/1‑71 (MQ=255)
             gCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTcc                                                   <  1:50498/70‑1 (MQ=255)
               gTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGt                                                 >  1:3642010/1‑70 (MQ=255)
                 gggCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTaaa                                              >  1:1080442/1‑71 (MQ=255)
                 gggCCATGAATAAGCGCGATGCAATAAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCg                                                  >  1:4426267/1‑67 (MQ=255)
                       tGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTc                                                               <  1:4063123/48‑1 (MQ=255)
                       tGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAAtcggtc                                        >  1:4225322/1‑71 (MQ=255)
                       tGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAAtcggt                                         >  1:4544801/1‑70 (MQ=255)
                         aaTAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCtt                                                             >  1:695010/1‑48 (MQ=255)
                         aaTAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGAt                                                          >  1:2232247/1‑51 (MQ=255)
                         aaTAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGAt                                                          >  1:1893321/1‑51 (MQ=255)
                         aaTAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGAt                                                          >  1:831051/1‑51 (MQ=255)
                          aTAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCATCTTGATCGCTTCCGTAAATCGGTcggc                                     <  1:405395/71‑1 (MQ=255)
                           tAAGCGCGATGCAATCAATCATATTCGGCAGTTTggcgggcg                                                                 <  1:2629866/42‑1 (MQ=255)
                             aGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTc                                                               >  1:2054804/1‑42 (MQ=255)
                             aGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTc                                                               >  1:1344516/1‑42 (MQ=255)
                             aGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCggcggc                                  >  1:4150845/1‑71 (MQ=255)
                             aGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCggcggc                                  >  1:1189390/1‑71 (MQ=255)
                              gcgcGATGCAATCAATCATATTCAGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCt                                 <  1:1956514/71‑1 (MQ=255)
                                gcgATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCa                               <  1:2968556/71‑1 (MQ=255)
                                      caatcaatCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAc                         <  1:4405058/71‑1 (MQ=255)
                                      caatcaatCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAc                         <  1:3235697/71‑1 (MQ=255)
                                      caatcaatCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAc                         <  1:1707400/71‑1 (MQ=255)
                                      caatcaatCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAc                         <  1:3861646/71‑1 (MQ=255)
                                      caatcaatCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAc                         <  1:1417346/71‑1 (MQ=255)
                                       aatcaatcATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAcg                        <  1:141899/71‑1 (MQ=255)
                                        atcaatcaTATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAAt                                             <  1:1510537/49‑1 (MQ=255)
                                        atcaatcaTATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAcgg                       <  1:1300679/71‑1 (MQ=255)
                                        atcaatcaTATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGAcgg                       >  1:3482416/1‑71 (MQ=255)
                                               atatTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGt                <  1:4046702/71‑1 (MQ=255)
                                                tatTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAAt                                             <  1:1361080/41‑1 (MQ=255)
                                                tatTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGt                >  1:2110119/1‑70 (MQ=255)
                                                tatTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGt                >  1:2734930/1‑70 (MQ=255)
                                                tatTCGGCAGTTTGGCGGGCATCTTGATCGCTTCCGTaa                                               <  1:708199/39‑1 (MQ=255)
                                                tatTCGGCAGTTTGGCGGGCATCTTGATCGCTTCCGTaa                                               <  1:22238/39‑1 (MQ=255)
                                                tatTCGGCAGTTTGGCGGGCATCTTGATCGCTTCCGTaa                                               <  1:2645296/39‑1 (MQ=255)
                                                  ttCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAAtcg                                           <  1:2512835/41‑1 (MQ=255)
                                                  ttCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGTCGGCGGTCGTaa              >  1:3626656/1‑70 (MQ=255)
                                                  ttCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGTaa              >  1:4499019/1‑70 (MQ=255)
                                                   tCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAAtcg                                           <  1:1309142/40‑1 (MQ=255)
                                                    cGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGTAACg            <  1:2528308/70‑1 (MQ=255)
                                                          tttGGCGGGCATCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCg                           <  1:2856683/49‑1 (MQ=255)
                                                          tttGGCGGGCATCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCg                           <  1:1536876/49‑1 (MQ=255)
                                                           ttGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGTAACGCACCaaaa    >  1:4526598/1‑71 (MQ=255)
                                                           ttGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGTAACGCACCaaaa    >  1:4361211/1‑71 (MQ=255)
                                                            tGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGTAACGCACCaaaa    >  1:1338447/1‑70 (MQ=255)
                                                              gcgggcgTCTTGATCGCTTCCGTAAATCGGTCGACGGCTCAGGCGACGGCGGTCGTAACGCACCAAAAcg  <  1:3435355/70‑1 (MQ=255)
                                                                    |                                                               
CACTTGCATGGAGGCGTGGGCCATGAATAAGCGCGATGCAATCAATCATATTCGGCAGTTTGGCGGGCGTCTTGATCGCTTCCGTAAATCGGTCGGCGGCTCAGGCGACGGCGGTCGTAACGCACCAAAACG  >  Exported/2625081‑2625212

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: