Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,632,142 C→T 22.7% M195I (ATG→ATA guaB ← IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,632,1420CT22.7% 41.0 / 8.6 22M195I (ATG→ATAguaBIMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (9/8);  new base T (3/2);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGC  >  W3110S.gb/2632073‑2632210
                                                                     |                                                                    
cAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCa                                                                     <  1:1802164/71‑1 (MQ=255)
      gCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGa                                                               <  1:2506810/71‑1 (MQ=255)
          tGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGcc                                                                 >  1:853310/1‑65 (MQ=255)
          tGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGcc                                                                 >  1:1387930/1‑65 (MQ=255)
           gCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATTATGCCGATCAgg                                                          >  1:511074/1‑71 (MQ=255)
                      cAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGt                                                <  1:1090804/70‑1 (MQ=255)
                      cAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGt                                                <  1:1089657/70‑1 (MQ=255)
                            ttCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAac                                         <  1:93993/71‑1 (MQ=255)
                              cGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATTATGCCGATCAGGTGGAATTCGTCATCAacca                                       >  1:1680575/1‑71 (MQ=255)
                              cGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAacc                                        >  1:365783/1‑70 (MQ=255)
                               ggTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATTATGCCGAt                                                              <  1:2543565/47‑1 (MQ=255)
                                  ttACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAg                                   >  1:1972240/1‑71 (MQ=255)
                                     cGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAgcgc                                <  1:2223863/71‑1 (MQ=255)
                                     cGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAgcgc                                >  1:1732669/1‑71 (MQ=255)
                                      gTTCCGCTTTCTGGAAGTCTTTCACGGTGATTATGCCGATCAGGTGGAATTCGTCATCAac                                         <  1:1173488/61‑1 (MQ=255)
                                            cTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAac                                         <  1:2588453/55‑1 (MQ=255)
                                                          ttCACGGTGATTATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTg           >  1:1245141/1‑71 (MQ=255)
                                                           tCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCaa                                          >  1:409363/1‑39 (MQ=255)
                                                           tCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCaa                                          >  1:999514/1‑39 (MQ=255)
                                                                   aTCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCtt                              >  1:1695077/1‑43 (MQ=255)
                                                                   aTCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGttttt               <  1:1705017/58‑1 (MQ=255)
                                                                   aTCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGc  >  1:2801415/1‑71 (MQ=255)
                                                                     |                                                                    
CAGACGGCCTTGCTCGTCTTTACAGGCGTTCGGTTTACGTTCCGCTTTCTGGAAGTCTTTCACGGTGATCATGCCGATCAGGTGGAATTCGTCATCAACCACCAGCGCTTTTTCAACGCGTTTTTCGTGCATTTTTGC  >  W3110S.gb/2632073‑2632210

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: