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

AACATGGTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCATCGTCGCCCG  >  ACB122/1472024‑1472286
                                                                                                                                  |                                                                                                                                      
aaCATGGTCGACGGCACTGAAAAGCCCGTGAACGATCGCGTGGGTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTCCAGATGATCGTCGACGGCGAACACCCCTAGGGTGACGTGTTCGCCGTGgcgc                                                                                                                                   <  8:290433/136‑1 (MQ=255)
  cATGGTCGACGTCACTGAAACGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGgcgcg                                                                                                                                  <  6:313115/135‑1 (MQ=255)
    tGGTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGGGCACACCCCAAGCGGGTGGTGGTGGCCGGGgcggggc                                                                                                                                >  2:124831/1‑131 (MQ=255)
     ggTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CAt                                                                                                                              <  1:370137/135‑1 (MQ=255)
       tCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CAt                                                                                                                              <  3:5538/133‑1 (MQ=255)
                       ggCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAg‑gcggc                                                                                                           <  6:255855/135‑1 (MQ=255)
                           gTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGAGGGCGCCCACCCCAAGGGGGACGTGTTGGCCGGGGCGGG‑CCTTGGCGGGCTACAG‑gcggcgca                                                                                                        >  1:231369/1‑132 (MQ=255)
                           gTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGCACACCCCCAGGGTGTCGTGTTCGCCGTGGCGCG‑CAATGCCGGGATACAG‑GCGGCGaa                                                                                                        >  6:210455/1‑134 (MQ=255)
                             gACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGACCACCCCAAGGGTGTCGGGTTCGCCGGGGCGCG‑CCGTGGCCGGGTAAAG‑GCGGCGcaaa                                                                                                      >  3:277868/1‑134 (MQ=255)
                              aCCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGCACCCCCCCAGGGTGCCGGGGTGGCCGGGGGGCG‑GCATGGCGGGGTTCAG‑GCGGGGcaaaaa                                                                                                    >  4:202674/1‑135 (MQ=255)
                                cGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGCAAACCCCAAGGGTGCCGTGTTAGCCGTGGCGCG‑CATTGCCGGGATACCG‑GCGGCGAAAAAAAcc                                                                                                 >  7:130632/1‑136 (MQ=255)
                                   gcgcGAGGCTACCGCCGAGGCGGGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCTAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GGGGCGAAAAAAACCAGc                                                                                              <  7:197461/136‑1 (MQ=255)
                                     gcgAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACCGCGAAAACCCCCACGGGGGCGGGGTAGCCCGGGCCCGCCATTG‑CGGGGTGCAGCGCGGCGAAAAAAccccccc                                                                                             >  7:12582/1‑132 (MQ=255)
                                        aGGCTACCGCCGAGGCGGGGGGGCGCATGTGGCCGCAGACCTTGCAGTTGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCGGTGGCGCG‑AATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGAcc                                                                                          <  4:123225/135‑1 (MQ=255)
                                              ccgccgAGGCGCGGGTGCGCATGTTGCCGCTGACCTTGCAGAGGATCGGCGACGGCTAACACCCCAAGGGTGACGGTTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCCGCGACCTGAtt                                                                                     <  8:4545/134‑1 (MQ=255)
                                              ccgccgAGGCGCGGGTGCGCATGTGGCCGCAGACCTTGCAGAGGAGCGTCGAGGGCCAACACCCCAAGGGGGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCTAAAAAAACCAGCGACCTGAtt                                                                                     <  7:162594/134‑1 (MQ=255)
                                                    aGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCtgtg                                                                             <  7:266753/136‑1 (MQ=255)
                                                               cgcAGGTGGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGGGAGGTGTGCGCGGTGGCGCG‑CTTTGCCGGGATACCG‑GCGGCTAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTg                                                                   <  2:215347/135‑1 (MQ=255)
                                                                     tggCCGCAGACCTTGCAGATGATCGGCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGGGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTc                                                             <  6:294298/133‑1 (MQ=255)
                                                                       gCCGCAGACCTTGCAGATGATCGTCGACGGCGCACACCCCAAGGGTGGCGTGTTAGCCGGGGCGCG‑CCTTTGCGGGATACAG‑GCGGCGCAAAAAACCACCGACCTGATTACGCTGTGCCCCCCCCGGCTGGTCAc                                                           >  6:164149/1‑135 (MQ=255)
                                                                       gCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CATGGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCAc                                                           <  8:56940/135‑1 (MQ=255)
                                                                          gCAGACCTTGCATATGCGCGTCGCGGGCGAACACCCAAGGGGTGAGGGGTTCGCGGGGGCGCG‑CATTGCCGGGAACCTG‑GCGGCTAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAg                                                        <  6:273531/135‑1 (MQ=255)
                                                                                     aGATGATCGTCGACGGCGAACACCCCCACGGTGTCGTGTTAGGCGGGGCGCG‑CCGTGCCGGGCTACCG‑GCGGCGCAAAAAAACAGCGACCTGATTCAGCTGTGCCCCCCCCGGATGCGCACCCGCGGCAAGGTCg                                             >  2:134941/1‑135 (MQ=255)
                                                                                           tcggcgCCGGCGAACCCCCCCCGGGTGCGGGGGGCGCCGGGGCGCG‑CTTGGCCGGGATACCG‑GCGGCGAAAAACCCCCGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCa                                       <  3:104878/135‑1 (MQ=255)
                                                                                             gtcgCCGGCTAACACCCCAAGGGTGAGGTGTGCGCCGGGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAACCCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAgcg                                    <  8:201705/136‑1 (MQ=255)
                                                                                              cccccGGCCAACACCCCAGGGGGGCCGTGTCCGCGGTGGCGCG‑CCTTGCCGGCAAACAC‑GCGCCTAAAAAACCCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAgcg                                    <  1:358937/131‑1 (MQ=255)
                                                                                                       aaCCCCCCACGGGTGACGTGTGCGCCGGGGCGCG‑CATGGCCGGGACACCG‑GCGGCTAAAAAAACCCGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGc                           <  3:122551/135‑1 (MQ=255)
                                                                                                       aaCACCCCATGGGTGACGGGTGCGCCGTGGCGCG‑CATTGCCGGCATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGc                           <  3:122550/135‑1 (MQ=255)
                                                                                                         cacCCCACGGGTGACGGGGTCGCCGTGGCGCG‑CATTGCCGGGATACCG‑GCGGCGAAAAAACCCCGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAg                         <  2:231369/135‑1 (MQ=255)
                                                                                                               aaGGGTGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGcc                   >  4:42315/1‑135 (MQ=255)
                                                                                                                 gggTGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGt                 >  6:1778/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgt        >  7:369931/1‑136 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgt        >  7:125931/1‑136 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  3:106696/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  2:178720/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  3:189256/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  5:128698/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcg         >  1:188699/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCG‑AAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgt        >  6:344919/1‑135 (MQ=255)
                                                                                                                         tgtTCGCCGTGGCGCG‑CATTGCCGGGATAAAG‑GCGGCTAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCGCCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtctt        >  7:354937/1‑134 (MQ=255)
                                                                                                                          gtgcGCCGGGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgtc       <  8:189235/133‑1 (MQ=255)
                                                                                                                           ttCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgtc       <  8:369931/135‑1 (MQ=255)
                                                                                                                           gcggCCGGGGCGCG‑CTTGGCCGGCAACCCG‑GGGGCTAAAAACCCCCGCGACCTGATTCCGCTGGGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCAtcgtc       <  5:379024/132‑1 (MQ=255)
                                                                                                                              gCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGACGAGGGGCAGGATGCCGTGCGCATCGTCGccc   >  8:148918/1‑136 (MQ=255)
                                                                                                                                cGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCATCGTCGCCCg  >  5:173197/1‑135 (MQ=255)
                                                                                                                                  |                                                                                                                                      
AACATGGTCGACGTCACTGAAAAGGCCGTGACCGAGCGCGAGGCTACCGCCGAGGCGCGGGTGCGCATGTTGCCGCAGACCTTGCAGATGATCGTCGACGGCGAACACCCCAAGGGTGACGTGTTCGCCGTGGCGCG‑CATTGCCGGGATACAG‑GCGGCGAAAAAAACCAGCGACCTGATTCCGCTGTGCCACCCGCTGATGCTCACCAGCGTCAAGGTCGAGCTCAGCGCCGAGGGGCAGGATGCCGTGCGCATCGTCGCCCG  >  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: ‑