Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 217,968 C→T 24.4% intergenic (‑60/‑46) yafU ← / → yafF predicted inner membrane protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported217,9680CT24.4% 78.9 / 24.0 45intergenic (‑60/‑46)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/12);  new base T (7/4);  total (29/16)
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.94e-01

TGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTAC  >  Exported/217900‑218031
                                                                    |                                                               
tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                               >  1:3764827/1‑71 (MQ=255)
tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCa                                                               >  1:3517786/1‑71 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                              <  1:4325023/70‑1 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                              <  1:2679150/70‑1 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                              <  1:3644044/70‑1 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                              <  1:3821054/70‑1 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                             >  1:259609/1‑71 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                             >  1:4511118/1‑71 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                             >  1:1708473/1‑71 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                             >  1:1639437/1‑71 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                             >  1:837663/1‑71 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt                                                             >  1:1398361/1‑71 (MQ=255)
  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTAGCCACAGCAATCCGAAACCACt                                                             >  1:2279135/1‑71 (MQ=255)
   aCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                            >  1:1790189/1‑71 (MQ=255)
   aCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                            >  1:2108639/1‑71 (MQ=255)
   aCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg                                                            >  1:919237/1‑71 (MQ=255)
        taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc                                                              <  1:3882919/64‑1 (MQ=255)
        taataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                       <  1:1283294/71‑1 (MQ=255)
        aaataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACg                                                       <  1:127667/70‑1 (MQ=255)
         aataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                                 <  1:2468057/60‑1 (MQ=255)
         aataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                                 <  1:4770672/60‑1 (MQ=255)
          ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc                                                                 >  1:2325212/1‑59 (MQ=255)
          ataaATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTg                                                     >  1:4016527/1‑71 (MQ=255)
                    tGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATa                                              <  1:3891386/68‑1 (MQ=255)
                     ggTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCt                                          >  1:197515/1‑71 (MQ=255)
                            ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAgg                                   >  1:2937070/1‑71 (MQ=255)
                             tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                  >  1:2751679/1‑71 (MQ=255)
                             tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                  >  1:3919022/1‑71 (MQ=255)
                               aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGt                                <  1:3954901/71‑1 (MQ=255)
                                      cGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                         <  1:2151764/71‑1 (MQ=255)
                                                ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa               >  1:4029888/1‑71 (MQ=255)
                                                ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa               >  1:975920/1‑71 (MQ=255)
                                                ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa               >  1:907351/1‑71 (MQ=255)
                                                ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa               >  1:3499276/1‑71 (MQ=255)
                                                ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa               >  1:3326733/1‑71 (MQ=255)
                                                ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGGaatgaa               >  1:3503472/1‑71 (MQ=255)
                                                  cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaat              >  1:4132721/1‑70 (MQ=255)
                                                    ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatga                <  1:3809680/66‑1 (MQ=255)
                                                       caGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                                  >  1:1450720/1‑45 (MQ=255)
                                                        aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                             <  1:3250859/49‑1 (MQ=255)
                                                        aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                             <  1:43725/49‑1 (MQ=255)
                                                        aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                             <  1:2109642/49‑1 (MQ=255)
                                                        aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                             <  1:2076933/49‑1 (MQ=255)
                                                             tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACCGGCGTCTGGACGTGGTAATGAATGAAGACGACTAc  <  1:1041421/71‑1 (MQ=255)
                                                              ccGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTAc  >  1:1504268/1‑70 (MQ=255)
                                                               cGAAACCACTGGCACGTGGAGAATAAGCTTCACAGGCGTCTGGACGTgg                      >  1:883595/1‑49 (MQ=255)
                                                               cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg                      >  1:73805/1‑49 (MQ=255)
                                                               cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg                      >  1:1376534/1‑49 (MQ=255)
                                                                    |                                                               
TGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTAC  >  Exported/217900‑218031

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: