Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,154,666 A→G 83.3% I23V (ATT→GTT)  yhfA → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,154,6660AG83.3% 102.6 / 12.9 42I23V (ATT→GTT) yhfAconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (19/16);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GCAAGCGCGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATGGACGGCAACTCAGGCGATAAAGCACCAAGTCCGATGGAAATGGTGTT  >  W3110S.gb/4154602‑4154721
                                                                |                                                       
gCAAGCGCGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGtttt                                                     <  1:2138947/69‑1 (MQ=255)
 cAAGCGCGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGATTTTa                                                    >  1:2715604/1‑69 (MQ=255)
   aGCGCGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTaa                                                   >  1:127804/1‑68 (MQ=255)
   aGCGCGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTaa                                                   >  1:3348964/1‑68 (MQ=255)
   aGCGCGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTaa                                                   <  1:2781535/68‑1 (MQ=255)
    gcgcGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATg                                                 <  1:2636865/69‑1 (MQ=255)
     cgcgAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTaa                                                   >  1:1310238/1‑66 (MQ=255)
     cgcgAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTaa                                                   >  1:182721/1‑66 (MQ=255)
     cgcgAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTaa                                                   >  1:1837427/1‑66 (MQ=255)
     cgcgAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTaa                                                   >  1:1405288/1‑66 (MQ=255)
      gcgAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATg                                                 <  1:1329758/67‑1 (MQ=255)
      gcgAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATg                                                 <  1:2557150/67‑1 (MQ=255)
      gcgAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATg                                                 >  1:3659977/1‑67 (MQ=255)
          gTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGtttt                                                     >  1:742389/1‑59 (MQ=255)
          gTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGtttt                                                     <  1:209926/59‑1 (MQ=255)
          gTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGtttt                                                     <  1:373988/59‑1 (MQ=255)
               gTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTc                                      <  1:2854387/69‑1 (MQ=255)
                 gggTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATGGACGGCAACTCa                                     >  1:434787/1‑68 (MQ=255)
                        aGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATg                                                 <  1:1763118/49‑1 (MQ=255)
                            ttAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGc                          >  1:90579/1‑68 (MQ=255)
                            ttAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGc                          >  1:2725654/1‑68 (MQ=255)
                            ttAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGc                          >  1:3016160/1‑68 (MQ=255)
                            ttAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGCa                         <  1:202953/69‑1 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATaaa                            >  1:2020778/1‑64 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATaa                             >  1:3207737/1‑63 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATaa                             >  1:3453963/1‑63 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAg                           >  1:3297524/1‑65 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAg                           >  1:912721/1‑65 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAg                           >  1:1748509/1‑65 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAg                           <  1:3296541/65‑1 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAg                           >  1:2025038/1‑65 (MQ=255)
                              aaCTTTTCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGCa                         <  1:1015893/67‑1 (MQ=255)
                                  tttCTGGGCGAATCCGCCTCTGGTCATCAGGtttt                                                     <  1:875998/35‑1 (MQ=255)
                                  tttCTGGGCGAATCCGCCTCTGGTCATCAGGtttt                                                     <  1:1851456/35‑1 (MQ=255)
                                  tttCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGAt                               >  1:451554/1‑57 (MQ=255)
                                  tttCTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATg                                                 <  1:2248643/39‑1 (MQ=255)
                                   ttCTGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGCACCAAGt                   >  1:2495110/1‑68 (MQ=255)
                                     cTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATTGACGGCAACTCAGGCGATAAAGCACCAAGTcc                 <  1:2842155/68‑1 (MQ=255)
                                      tGGGCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGCACCAAGTCCg                <  1:2919243/68‑1 (MQ=255)
                                       gggCGAATCCGCCTCTGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGCACCAAGTCCg                <  1:956458/67‑1 (MQ=255)
                                                   ctctGGTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGCACCAAGTCCGATGGAAATGGTGtt  <  1:2116210/69‑1 (MQ=255)
                                                        gTCATCAGGTTTTAATGGACGGCAACTCAGGCGATAAAGCACCAAGTCCGATGGAAATGgtg    <  1:1459670/62‑1 (MQ=255)
                                                                |                                                       
GCAAGCGCGAGTGAAGTGGGTCGAAGGGTTAACTTTTCTGGGCGAATCCGCCTCTGGTCATCAGATTTTAATGGACGGCAACTCAGGCGATAAAGCACCAAGTCCGATGGAAATGGTGTT  >  W3110S.gb/4154602‑4154721

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