Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 301,283 A→C 10.9% V172G (GTG→GGG)  yagT ← predicted xanthine dehydrogenase, 2Fe‑2S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb301,2830AC10.9% 117.5 / 11.0 55V172G (GTG→GGG) yagTpredicted xanthine dehydrogenase, 2Fe‑2S subunit
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base C (3/3);  total (25/30)
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

GTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATA  >  W3110S.gb/301217‑301348
                                                                  |                                                                 
gttTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGc                                                               >  1:36216/1‑71 (MQ=255)
gttTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGc                                                               >  1:3368661/1‑71 (MQ=255)
   tCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCAc                                                                  >  1:1440820/1‑65 (MQ=255)
   tCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTAc                                                            >  1:3732658/1‑71 (MQ=255)
       gAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACt                                                           <  1:3635589/68‑1 (MQ=255)
       gAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACt                                                           <  1:3009258/68‑1 (MQ=255)
          cGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGc                                                               <  1:2309910/61‑1 (MQ=255)
          cGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGc                                                               <  1:153853/61‑1 (MQ=255)
             aaaCCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCa                                                                   <  1:3664576/54‑1 (MQ=255)
              aaCCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCCCCGCTACTGATGAGCaaa                                                 >  1:285191/1‑71 (MQ=255)
              aaCCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCaaa                                                 <  1:1112477/71‑1 (MQ=255)
              aaCCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCaaa                                                 <  1:705944/71‑1 (MQ=255)
               aCCAAATCGACCGTGACGTGACTGGGAATGCCCTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAAt                                                <  1:2112096/71‑1 (MQ=255)
                ccAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCaaa                                                 <  1:3291747/69‑1 (MQ=255)
                 cAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCa                                                                   >  1:4045708/1‑50 (MQ=255)
                   aaTCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGc                                            <  1:2475283/71‑1 (MQ=255)
                    aTCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTa                                                             >  1:128706/1‑53 (MQ=255)
                    aTCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTa                                                             >  1:489764/1‑53 (MQ=255)
                       gACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCgg                                        <  1:2345008/71‑1 (MQ=255)
                          cctgacgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGcccc                                                                 <  1:527671/41‑1 (MQ=38)
                            tgacgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTg                                   <  1:820650/71‑1 (MQ=255)
                             gacgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCCCCGCTACtgatga                                                      >  1:2448013/1‑51 (MQ=255)
                             gacgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTg                                   >  1:4028378/1‑70 (MQ=255)
                              acgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGc                                                               <  1:2032751/41‑1 (MQ=255)
                              acgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGc                                  <  1:1748404/70‑1 (MQ=255)
                              acgtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCa                                 <  1:865858/71‑1 (MQ=255)
                                gtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCAc                                                                  >  1:2534917/1‑36 (MQ=255)
                                gtgacTGGGAATGCCGTCCTGAATCTCTTTTAGCAc                                                                  >  1:1351747/1‑36 (MQ=255)
                                      gggAATGCCGTCCTGAATCTCTTTTAGCCCCGCTACTGATGAGCAAATTTGcc                                           <  1:2034872/53‑1 (MQ=255)
                                       ggAATGCCGTCCTGAATCTCTTTTAGCCCCGCTACTGATGAGCAAATTTGcc                                           <  1:3375388/52‑1 (MQ=255)
                                       ggAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTTCAGTAGCCg                          <  1:4607378/69‑1 (MQ=255)
                                       ggAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCg                          <  1:3830378/69‑1 (MQ=255)
                                          aTGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAg                                >  1:1873791/1‑60 (MQ=255)
                                          aTGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGt                               >  1:3211641/1‑61 (MQ=255)
                                          aTGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGt                               >  1:73781/1‑61 (MQ=255)
                                           tGCCGTCCTGAATCTCTTTTAGCACCGCTACtgatg                                                       <  1:2540998/36‑1 (MQ=255)
                                             ccGTCCTGAATCTCTTTTAGCCCCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTgg                    >  1:1149696/1‑69 (MQ=255)
                                             ccGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCAc                       >  1:4542570/1‑66 (MQ=255)
                                             ccGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTgg                    >  1:748954/1‑69 (MQ=255)
                                              cGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAg                 <  1:2994137/71‑1 (MQ=255)
                                                tCCTGAATCTCTTTTAGCCCCGCTACTGATGAGCAAATTTGcc                                           <  1:1928571/43‑1 (MQ=255)
                                                  cTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTa                              <  1:1389519/54‑1 (MQ=255)
                                                  cTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCcat             >  1:2143418/1‑71 (MQ=255)
                                                   tGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCcat             >  1:439872/1‑70 (MQ=255)
                                                   tGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCcat             >  1:2257226/1‑70 (MQ=255)
                                                       tctcTTTTAGCACCGCTACTGATGAGCAAATTTGCCCgg                                        <  1:866826/39‑1 (MQ=255)
                                                       tctcTTTTAGCACCGCTACTGATGAGCAAATTTGCCCgg                                        <  1:2021903/39‑1 (MQ=255)
                                                       tctcTTTTAGCACCGCTACTGATGAGCAAATTTGCCCg                                         <  1:1073862/38‑1 (MQ=255)
                                                       tctcTTTTAGCACCGCTACTGATGAGCAAATTTGCCCg                                         <  1:189348/38‑1 (MQ=255)
                                                       tctcTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGc                                  <  1:4634602/45‑1 (MQ=255)
                                                       tctcTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCcatca           >  1:1641291/1‑68 (MQ=255)
                                                        ctctTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCt       >  1:1770846/1‑71 (MQ=255)
                                                         tctTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCAc                       >  1:2332513/1‑54 (MQ=255)
                                                           ttttAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCcat             <  1:3063432/62‑1 (MQ=255)
                                                             ttAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCcatca           >  1:3729322/1‑62 (MQ=255)
                                                             ttAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCcatca           >  1:4468168/1‑62 (MQ=255)
                                                              tAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATa  <  1:411538/70‑1 (MQ=255)
                                                              tAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATa  <  1:2665450/70‑1 (MQ=255)
                                                                  |                                                                 
GTTTCTGGAGCGGAAACCAAATCGACCGTGACGTGACTGGGAATGCCGTCCTGAATCTCTTTTAGCACCGCTACTGATGAGCAAATTTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATA  >  W3110S.gb/301217‑301348

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: