Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,887,887 T→A 8.2% intergenic (+52/+27) ytfK → / ← ytfL conserved hypothetical protein/predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,887,8870TA 8.2% 132.6 / 4.3 49intergenic (+52/+27)ytfK/ytfLconserved hypothetical protein/predicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base A (2/2);  total (22/27)
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.78e-01

AAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTT  >  Exported/2887819‑2887955
                                                                    |                                                                    
aaaCAGAGATGGCTTACCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAtt                                                                     >  1:3039140/1‑70 (MQ=255)
 aaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                    <  1:2808253/70‑1 (MQ=255)
  aCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                    >  1:3436325/1‑69 (MQ=255)
  aCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTc                                                                    >  1:3209703/1‑69 (MQ=255)
     gagaTGGCTAAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAg                                                               <  1:671163/71‑1 (MQ=255)
         tGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCa                                                                   <  1:3492950/63‑1 (MQ=255)
         tGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCa                                                                   <  1:2237517/63‑1 (MQ=255)
          ggCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCa                                                                   <  1:3055370/62‑1 (MQ=255)
           gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTg                                                         <  1:3052610/71‑1 (MQ=255)
            cTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCa                                                                   <  1:2899090/60‑1 (MQ=255)
            cTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGa                                                        <  1:2526104/71‑1 (MQ=255)
                  cAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTc                                                  <  1:2480938/71‑1 (MQ=255)
                    aagtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTg                                                <  1:2999396/71‑1 (MQ=255)
                     agtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgg                                               >  1:2096857/1‑71 (MQ=255)
                        gctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATg                                                      >  1:2334233/1‑61 (MQ=255)
                         ctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTg                                                <  1:1839646/66‑1 (MQ=255)
                          tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTgg                                               <  1:3409448/66‑1 (MQ=255)
                            tgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTa                                         <  1:1739226/70‑1 (MQ=255)
                             gcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcg                                       >  1:3695203/1‑71 (MQ=255)
                               agtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                      >  1:1246275/1‑70 (MQ=255)
                                  aaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACg                                  >  1:1135449/1‑71 (MQ=255)
                                  aaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACg                                  >  1:677175/1‑71 (MQ=255)
                                   aaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACg                                  <  1:2621520/70‑1 (MQ=255)
                                   aaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGc                                 >  1:2718489/1‑71 (MQ=255)
                                     aagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCtt                               >  1:1403656/1‑71 (MQ=255)
                                      agtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCt                                <  1:630470/69‑1 (MQ=255)
                                             tgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttt                       <  1:867860/71‑1 (MQ=255)
                                              gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGAc                                   >  1:1981551/1‑58 (MQ=255)
                                              gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                      >  1:1842707/1‑71 (MQ=255)
                                               cagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttta                      >  1:2830728/1‑70 (MQ=255)
                                                agtaaTAAGACGGCGCCTGATTCAGGAGCCGTTGATGTTTCTggg                                              <  1:1941298/45‑1 (MQ=255)
                                                 gtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCttt                              <  1:1909355/60‑1 (MQ=255)
                                                 gtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatct                   >  1:3169918/1‑71 (MQ=255)
                                                 gtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatct                   >  1:623067/1‑71 (MQ=255)
                                                      aaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCt                                                 >  1:629546/1‑36 (MQ=255)
                                                      aaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCttt                              <  1:1913534/55‑1 (MQ=255)
                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCt                                                 <  1:3543375/35‑1 (MQ=255)
                                                       aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCt                                                 <  1:1161571/35‑1 (MQ=255)
                                                          cGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGt            <  1:3293707/69‑1 (MQ=255)
                                                          cGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGt            >  1:1978962/1‑69 (MQ=255)
                                                          cGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGt            >  1:1201682/1‑69 (MQ=255)
                                                          cGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTc           <  1:1255849/70‑1 (MQ=255)
                                                           ggCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg             <  1:960345/67‑1 (MQ=255)
                                                           ggCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttta                >  1:1537436/1‑64 (MQ=255)
                                                           ggCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatcttta                >  1:1336525/1‑64 (MQ=255)
                                                           ggCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAg               <  1:2442913/65‑1 (MQ=255)
                                                             cTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGc       <  1:430154/71‑1 (MQ=255)
                                                              tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGc       <  1:3051359/70‑1 (MQ=255)
                                                               ccTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGt            >  1:1856684/1‑64 (MQ=255)
                                                                  gATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGttt  >  1:976555/1‑71 (MQ=255)
                                                                    |                                                                    
AAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTT  >  Exported/2887819‑2887955

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: