Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,472,154 T→G 27.9% V55G (GTG→GGG)  moaC → Cyclic pyranopterin monophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,472,1540TG27.9% 68.6 / 2.4 45V55G (GTG→GGG) moaCCyclic pyranopterin monophosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (15/16);  new base G (7/7);  total (22/23)
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

AACATGGTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCATCGTCGCCCG  >  ACB122/1472024‑1472286
                                                                                                                                  |                                                                                                                                      
aaCATGGTCGACGGCACTGAAAAGCCCGTGAACGATCGCGTGGGTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTCCAGATGATCGTCGACGGCGAACACCCCTAGGGTGACGTGTTCGCCGTGgcgc                                                                                                                                   <  8:290433/136‑1 (MQ=255)
  cATGGTCGACGTCACTGAAACGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGgcgcg                                                                                                                                  <  6:313115/135‑1 (MQ=255)
    tGGTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGGGCACACCCCAAGCGGGTGGTGGTGGCCGGGgcggggc                                                                                                                                >  2:124831/1‑131 (MQ=255)
     ggTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCAt                                                                                                                              <  1:370137/135‑1 (MQ=255)
       tCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCAt                                                                                                                              <  3:5538/133‑1 (MQ=255)
                       ggCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAggcggc                                                                                                           <  6:255855/135‑1 (MQ=255)
                           gTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGAGGGCGCCCACCCCAAGGGGGACGTGTTGGCCGGGGCGGGCCTTGGCGGGCTACAGgcggcgca                                                                                                        >  1:231369/1‑132 (MQ=255)
                           gTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGCACACCCCCAGGGTGTCGTGTTCGCCGTGGCGCGCAATGCCGGGATACAGGCGGCGaa                                                                                                        >  6:210455/1‑134 (MQ=255)
                             gACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGACCACCCCAAGGGTGTCGGGTTCGCCGGGGCGCGCCGTGGCCGGGTAAAGGCGGCGcaaa                                                                                                      >  3:277868/1‑134 (MQ=255)
                              aCCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGCACCCCCCCAGGGTGCCGGGGTGGCCGGGGGGCGGCATGGCGGGGTTCAGGCGGGGcaaaaa                                                                                                    >  4:202674/1‑135 (MQ=255)
                                cGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGCAAACCCCAAGGGTGCCGTGTTAGCCGTGGCGCGCATTGCCGGGATACCGGCGGCGAAAAAAAcc                                                                                                 >  7:130632/1‑136 (MQ=255)
                                   gcgcGAGGCTACCGCCGAGGCGGGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCTAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAGGGGGCGAAAAAAACCAGc                                                                                              <  7:197461/136‑1 (MQ=255)
                                     gcgAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACCGCGAAAACCCCCACGGGGGCGGGGTAGCCCGGGCCCGCCATTGCGGGGTGCAGCGCGGCGAAAAAAccccccc                                                                                             >  7:12582/1‑132 (MQ=255)
                                        aGGCTACCGCCGAGGCGGGGGGGCGCATGTGGCCGCAGACCTTGCAGTTGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCGGTGGCGCGAATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGAcc                                                                                          <  4:123225/135‑1 (MQ=255)
                                              ccgccgAGGCGCGGGTGCGCATGTTGCCGCTGACCTTGCAGAGGATCGGCGACGGCTAACACCCCAAGGGTGACGGTTTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCCGCGACCTGAtt                                                                                     <  8:4545/134‑1 (MQ=255)
                                              ccgccgAGGCGCGGGTGCGCATGTGGCCGCAGACCTTGCAGAGGAGCGTCGAGGGCCAACACCCCAAGGGGGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCTAAAAAAACCAGCGACCTGAtt                                                                                     <  7:162594/134‑1 (MQ=255)
                                                    aGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCtgtg                                                                             <  7:266753/136‑1 (MQ=255)
                                                               cgcAGGTGGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGGGAGGTGTGCGCGGTGGCGCGCTTTGCCGGGATACCGGCGGCTAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTg                                                                   <  2:215347/135‑1 (MQ=255)
                                                                     tggCCGCAGACCTTGCAGATGATCGGCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGGGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTc                                                             <  6:294298/133‑1 (MQ=255)
                                                                       gCCGCAGACCTTGCAGATGATCGTCGACGGCGCACACCCCAAGGGTGGCGTGTTAGCCGGGGCGCGCCTTTGCGGGATACAGGCGGCGCAAAAAACCACCGACCTGATTACGCTGTGCCCCCCCCGGCTGGTCAc                                                           >  6:164149/1‑135 (MQ=255)
                                                                       gCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCATGGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCAc                                                           <  8:56940/135‑1 (MQ=255)
                                                                          gCAGACCTTGCATATGCGCGTCGCGGGCGAACACCCAAGGGGTGAGGGGTTCGCGGGGGCGCGCATTGCCGGGAACCTGGCGGCTAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAg                                                        <  6:273531/135‑1 (MQ=255)
                                                                                     aGATGATCGTCGACGGCGAACACCCCCACGGTGTCGTGTTAGGCGGGGCGCGCCGTGCCGGGCTACCGGCGGCGCAAAAAAACAGCGACCTGATTCAGCTGTGCCCCCCCCGGATGCGCACCCGCGGCAAGGTCg                                             >  2:134941/1‑135 (MQ=255)
                                                                                           tcggcgCCGGCGAACCCCCCCCGGGTGCGGGGGGCGCCGGGGCGCGCTTGGCCGGGATACCGGCGGCGAAAAACCCCCGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCa                                       <  3:104878/135‑1 (MQ=255)
                                                                                             gtcgCCGGCTAACACCCCAAGGGTGAGGTGTGCGCCGGGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAACCCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAgcg                                    <  8:201705/136‑1 (MQ=255)
                                                                                              cccccGGCCAACACCCCAGGGGGGCCGTGTCCGCGGTGGCGCGCCTTGCCGGCAAACACGCGCCTAAAAAACCCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAgcg                                    <  1:358937/131‑1 (MQ=255)
                                                                                                       aaCCCCCCACGGGTGACGTGTGCGCCGGGGCGCGCATGGCCGGGACACCGGCGGCTAAAAAAACCCGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGc                           <  3:122551/135‑1 (MQ=255)
                                                                                                       aaCACCCCATGGGTGACGGGTGCGCCGTGGCGCGCATTGCCGGCATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGc                           <  3:122550/135‑1 (MQ=255)
                                                                                                         cacCCCACGGGTGACGGGGTCGCCGTGGCGCGCATTGCCGGGATACCGGCGGCGAAAAAACCCCGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAg                         <  2:231369/135‑1 (MQ=255)
                                                                                                               aaGGGTGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGcc                   >  4:42315/1‑135 (MQ=255)
                                                                                                                 gggTGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGt                 >  6:1778/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgt        >  7:369931/1‑136 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgt        >  7:125931/1‑136 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  3:106696/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  2:178720/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  3:189256/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  5:128698/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  1:188699/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgt        >  6:344919/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCGCATTGCCGGGATAAAGGCGGCTAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCGCCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtctt        >  7:354937/1‑134 (MQ=255)
                                                                                                                          gtgcGCCGGGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgtc       <  8:189235/133‑1 (MQ=255)
                                                                                                                           ttCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgtc       <  8:369931/135‑1 (MQ=255)
                                                                                                                           gcggCCGGGGCGCGCTTGGCCGGCAACCCGGGGGCTAAAAACCCCCGCGACCTGATTCCGCTGGGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgtc       <  5:379024/132‑1 (MQ=255)
                                                                                                                              gCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGACGAGGGGCAGGATGCCGTGCGCATCGTCGccc   >  8:148918/1‑136 (MQ=255)
                                                                                                                                cGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCATCGTCGCCCg  >  5:173197/1‑135 (MQ=255)
                                                                                                                                  |                                                                                                                                      
AACATGGTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCGCATTGCCGGGATACAGGCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCATCGTCGCCCG  >  ACB122/1472024‑1472286

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