Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 670,287 T→C 10.5% I15V (ATC→GTC)  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,2870TC10.5% 147.3 / 8.0 57I15V (ATC→GTC) cobCpredicted alpha‑ribazole‑5'‑P phosphatase
Reads supporting (aligned to +/- strand):  ref base T (29/22);  new base C (4/2);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTAT  >  W3110S.gb/670221‑670353
                                                                  |                                                                  
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)
                            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:1170282/1‑71 (MQ=255)
                                              aCCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCt                  >  1:485875/1‑71 (MQ=255)
                                              aCCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCt                  >  1:58817/1‑71 (MQ=255)
                                                 gCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGcat                          <  1:660887/60‑1 (MQ=255)
                                                  cTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCg               >  1:2516164/1‑70 (MQ=255)
                                                  cTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCg               >  1:577512/1‑70 (MQ=255)
                                                  cTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCg               >  1:1135903/1‑70 (MQ=255)
                                                    gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                            <  1:462952/39‑1 (MQ=255)
                                                    gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                            <  1:3211982/55‑1 (MQ=255)
                                                    gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                            <  1:794423/55‑1 (MQ=255)
                                                    gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                            <  1:3322816/55‑1 (MQ=255)
                                                    gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                            <  1:1370783/55‑1 (MQ=255)
                                                    gTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGc                            <  1:3651249/55‑1 (MQ=255)
                                                     tAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                            <  1:997286/38‑1 (MQ=255)
                                                     tAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                            <  1:2330246/38‑1 (MQ=255)
                                                     tAAAGACCATCGATATTCGCTTGCGTTTCACCATGACg                                            <  1:1323837/38‑1 (MQ=255)
                                                       aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGcatca                        >  1:2930891/1‑56 (MQ=255)
                                                       aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa         >  1:1609917/1‑71 (MQ=255)
                                                       aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa         >  1:1229764/1‑71 (MQ=255)
                                                       aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa         >  1:2972429/1‑71 (MQ=255)
                                                       aaGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGa         >  1:1883955/1‑71 (MQ=255)
                                                          aCCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAg      <  1:2407914/71‑1 (MQ=255)
                                                           ccATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAgg     >  1:1734978/1‑71 (MQ=255)
                                                           ccATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAgg     >  1:645654/1‑71 (MQ=255)
                                                            cATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAgg     >  1:2245093/1‑70 (MQ=255)
                                                            cATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGt    <  1:3758292/71‑1 (MQ=255)
                                                               cGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGcat                          >  1:3296631/1‑46 (MQ=255)
                                                                gATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTAt  >  1:1689415/1‑69 (MQ=255)
                                                                 atatTCGCTTGCGTTTCACCATGACGAATTAACCACAg                                >  1:3369576/1‑38 (MQ=255)
                                                                 atatTCGCTTGCGTTTCACCATGACGAATTAACCACAg                                >  1:1001339/1‑38 (MQ=255)
                                                                  |                                                                  
TTTGCGCTTGCTCAATACCGCGCGCGGTCAGGGGGGTGGGCGCATGACCGCTGTAAAGACCATCGATATTCGCTTGCGTTTCACCATGACGAATTAACCACAGTCGCATCATGCCCTCCGTAACGACAGGTAT  >  W3110S.gb/670221‑670353

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