Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 670,276 A→C 11.5% L18L (CTT→CTG cobC ← predicted alpha‑ribazole‑5'‑P phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb670,2760AC11.5% 106.6 / 11.9 52L18L (CTT→CTGcobCpredicted alpha‑ribazole‑5'‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base C (4/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGA  >  W3110S.gb/670214‑670346
                                                              |                                                                      
tGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATc                                                                >  1:3335330/1‑71 (MQ=255)
tGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATc                                                                <  1:1450066/71‑1 (MQ=255)
tGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATc                                                                >  1:3654181/1‑71 (MQ=255)
  cAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGa                                                              >  1:289362/1‑71 (MQ=255)
    gATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCa                                                                  <  1:2104901/65‑1 (MQ=255)
       tttGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTc                                                         >  1:3706649/1‑71 (MQ=255)
        ttGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCTCATGACCGCTGTAAAGACCATCGATATTc                                                         >  1:941193/1‑70 (MQ=255)
        ttGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTc                                                         >  1:2605643/1‑70 (MQ=255)
          gcgcTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCt                                                      >  1:1487757/1‑71 (MQ=255)
            gcTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAATACCATCGATATTCGCt                                                      <  1:3113738/69‑1 (MQ=255)
              ttGCTCAATACCGCGCGCGGTCAGGGGGGTGGGGGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCg                                                  >  1:101836/1‑71 (MQ=255)
                gCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGata                                                            >  1:1439691/1‑59 (MQ=255)
                    aaTACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTc                                              <  1:2706783/69‑1 (MQ=255)
                    aaTACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCAc                                            <  1:301933/71‑1 (MQ=255)
                      tACCGCGCGCGGTCAGGGTGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTc                                              >  1:686517/1‑67 (MQ=255)
                      tACCGCGCGCGGTCAGGGGTGTGGGCGCATGACCGCTGTACAGACCATCGACATTCGCTTGCGtt                                                >  1:2341359/1‑65 (MQ=255)
                      tACCGCGCGCGGTCAGGGGTGTGGGCGCATGACCGCTGTACAGACCATCGACATTCGCTTGCGTTTCAcaa                                          >  1:2539449/1‑69 (MQ=255)
                      tACCGCGCGCGGTCAGGGGTGTGGGCGCATGACCGCTGTACAGACCATCGACATTCGCTTGCGTTTCACCa                                          >  1:14460/1‑71 (MQ=255)
                      tACCGCGCGCGGTCAGGGGTGTGGGCGCATGACCGCTGTACAGACCATCGACATTCGCTTGCGTTTCACCa                                          >  1:2376128/1‑71 (MQ=255)
                        ccGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCAt                                         >  1:2638779/1‑70 (MQ=255)
                        ccGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCAt                                         >  1:1217987/1‑70 (MQ=255)
                        ccGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCCCTTGCGTTTCACCAt                                         >  1:3422796/1‑70 (MQ=255)
                        ccGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATAGATATTCGCTTGCGTTTCACCAt                                         >  1:3258190/1‑70 (MQ=255)
                           cgcgcgGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                     <  1:1774159/71‑1 (MQ=255)
                               cgGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAAtt                                 <  1:507740/71‑1 (MQ=255)
                               cgGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAAtt                                 <  1:3313006/71‑1 (MQ=255)
                               cgGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAAt                                  <  1:3693912/70‑1 (MQ=255)
                                  tCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGa                                                              <  1:2367696/39‑1 (MQ=255)
                                   cAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAAc                              >  1:1783570/1‑70 (MQ=255)
                                    aGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACca                            >  1:1117042/1‑71 (MQ=255)
                                      gggtgTGGGCGCATGACCGCTGTACAGACCATCGACATTCGCTTGCGTTTCACCATGACGAATTAACcaca                          <  1:1983816/71‑1 (MQ=255)
                                      gggtgTGGGCGCATGACCGCTGTACAGACCATCGACATTCGCTTGCGTTTCACCATGACGAATTAACcaca                          <  1:3453196/71‑1 (MQ=255)
                                       ggggTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAg                         <  1:3800061/71‑1 (MQ=255)
                                          gTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCg                      <  1:1565210/71‑1 (MQ=255)
                                             ggCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGcat                   >  1:2360725/1‑71 (MQ=255)
                                               cgcATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGcatca                 <  1:1848317/71‑1 (MQ=255)
                                                     aCCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCt           >  1:485875/1‑71 (MQ=255)
                                                     aCCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCt           >  1:58817/1‑71 (MQ=255)
                                                     aCCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCt           >  1:1170282/1‑71 (MQ=255)
                                                        gCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGcat                   <  1:660887/60‑1 (MQ=255)
                                                         cTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCg        >  1:577512/1‑70 (MQ=255)
                                                         cTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCg        >  1:1135903/1‑70 (MQ=255)
                                                         cTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCg        >  1:2516164/1‑70 (MQ=255)
                                                           gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                     <  1:462952/39‑1 (MQ=255)
                                                           gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                     <  1:3322816/55‑1 (MQ=255)
                                                           gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                     <  1:3651249/55‑1 (MQ=255)
                                                           gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                     <  1:3211982/55‑1 (MQ=255)
                                                           gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                     <  1:1370783/55‑1 (MQ=255)
                                                           gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                     <  1:794423/55‑1 (MQ=255)
                                                            tAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                     <  1:997286/38‑1 (MQ=255)
                                                            tAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                     <  1:1323837/38‑1 (MQ=255)
                                                            tAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                     <  1:2330246/38‑1 (MQ=255)
                                                              aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGcatca                 >  1:2930891/1‑56 (MQ=255)
                                                              aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa  >  1:1229764/1‑71 (MQ=255)
                                                              aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa  >  1:1609917/1‑71 (MQ=255)
                                                              aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa  >  1:1883955/1‑71 (MQ=255)
                                                              aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa  >  1:2972429/1‑71 (MQ=255)
                                                              |                                                                      
TGCAGATTTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGA  >  W3110S.gb/670214‑670346

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: