Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,405,008 C→T 22.2% E406E (GAG→GAA nuoF ← NADH:ubiquinone oxidoreductase, chain F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,405,0080CT22.2% 70.4 / 15.5 36E406E (GAG→GAAnuoFNADH:ubiquinone oxidoreductase, chain F
Reads supporting (aligned to +/- strand):  ref base C (10/18);  new base T (5/3);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAGAATCGACACAGTTGCT  >  W3110S.gb/2404940‑2405072
                                                                    |                                                                
ttGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGTTc                                                                <  1:4003792/71‑1 (MQ=255)
  gCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCa                                                              >  1:853934/1‑71 (MQ=255)
  gCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCa                                                              >  1:2482847/1‑71 (MQ=255)
   cTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCa                                                              <  1:4500757/70‑1 (MQ=255)
   cTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCa                                                              <  1:228507/70‑1 (MQ=255)
   cTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCa                                                              <  1:164809/70‑1 (MQ=255)
   cTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCAc                                                             >  1:1898450/1‑71 (MQ=255)
     gAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTc                                                                <  1:1321708/66‑1 (MQ=255)
        cGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCAc                                                        >  1:1302374/1‑71 (MQ=255)
              tttGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAg                                                     <  1:1910070/68‑1 (MQ=255)
              tttGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAg                                                     <  1:1904684/68‑1 (MQ=255)
              tttGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTg                                                  >  1:1201278/1‑71 (MQ=255)
               ttGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCAc                                                             >  1:3295688/1‑59 (MQ=255)
                   ttCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGTTCCACTGCAc                                                        <  1:2246123/60‑1 (MQ=255)
                     cccGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGTTCCACTGCACCa                                                      >  1:3601893/1‑60 (MQ=255)
                     cccGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGAGTggg                                            <  1:296495/70‑1 (MQ=255)
                           tCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCa                                                              >  1:2644365/1‑46 (MQ=255)
                                ttcttcGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGTGCACAGAAAGttt                                <  1:1008204/71‑1 (MQ=255)
                                ttcttcGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGttt                                <  1:394533/71‑1 (MQ=255)
                                  cttcGCGGAAATATTTGATGGCGCTCTGTAACGGTTCCACt                                                            >  1:1287880/1‑41 (MQ=255)
                                  cttcGCGGAAATATTTGATGGCGCTCTGTAACGGTTCCACt                                                            >  1:3781570/1‑41 (MQ=255)
                                  cttcGCGGAAATATTTGATGGCGCTCTGTAACGGTTCCACt                                                            >  1:2990731/1‑41 (MQ=255)
                                     cgcgGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTa                              <  1:3560205/68‑1 (MQ=255)
                                         gAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCAcc                                                       <  1:2341212/39‑1 (MQ=255)
                                           aaTATTTGATGGCGCTCTGTAACGGTTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCt                     <  1:2161054/71‑1 (MQ=255)
                                            atatTTGATGGCGCTCTGTAACGGTTCCACTGCACCAGGt                                                   >  1:4099175/1‑40 (MQ=255)
                                            atatTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTAccc                           >  1:2385450/1‑64 (MQ=255)
                                               tttGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTAc                             >  1:832158/1‑59 (MQ=255)
                                               tttGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAg                  <  1:2119242/70‑1 (MQ=255)
                                                 tGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGc                                                 <  1:1310716/37‑1 (MQ=255)
                                                   aTGGCTCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAGAATc              <  1:1461139/70‑1 (MQ=255)
                                                   aTGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAGAATc              <  1:655954/70‑1 (MQ=255)
                                                       cgcTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAGAATCGacac         <  1:504671/71‑1 (MQ=255)
                                                            tGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGt                                  <  1:1364078/41‑1 (MQ=255)
                                                             gTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAGa                 <  1:1622177/57‑1 (MQ=255)
                                                               aaCGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAGAATCGACACAGTTGCt  >  1:3720850/1‑70 (MQ=255)
                                                                    |                                                                
TTGCTGAACGGCTGTTTGATTCCCGCCTCAAATTCTTCGCGGAAATATTTGATGGCGCTCTGTAACGGCTCCACTGCACCAGGTGCGTGGGCACAGAAAGTTTTACCCGGGCCTAAGAATCGACACAGTTGCT  >  W3110S.gb/2404940‑2405072

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: