Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,500,221 A→C 33.9% G58G (GGT→GGG)  cysG_1 ← Siroheme synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,500,2210AC33.9% 29.9 / 3.3 44G58G (GGT→GGG) cysG_1Siroheme synthase
Reads supporting (aligned to +/- strand):  ref base A (19/8);  new base C (12/4);  total (32/12)
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

ACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCG  >  ACB122/4500090‑4500355
                                                                                                                                   |                                                                                                                                      
acaTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCg                                                                                                                                     <  4:453935/135‑1 (MQ=255)
 caTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGc                                                                                                                                    <  8:464403/135‑1 (MQ=255)
  aTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCCACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCCCCTCCCCGcc                                                                                                                                   >  8:281792/1‑135 (MQ=255)
         ggCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTGAGTCGCCGCGATCACCAGCCGGCAACCGCCCCGGTCGGCTGCCTGATTGCCCCGCCCCAGCACCTCCCCGCCCCCGTcc                                                                                                                            >  5:94219/1‑135 (MQ=255)
               cTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGcc                                                                                                                      <  3:221089/135‑1 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCCCCAGCCGGCAACCGACCCGGTCGGCTGCCGGCTAGCCCCGCCCCCGCACCGCCCCGCCCCCGGCCCGGGCCCCCGCGGccccc                                                                                                          >  1:166759/1‑132 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCAACCGACCCGGTCGGCTGCCGGCTAGCCACGCACCCCCACCACACCGCCCCCGTCCCGGGCCCACCCCGCCCAc                                                                                                          >  2:277673/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCCCGCACCAGCACCTCACCGCCCCCGTCCCGGGCCCACCCGGCCCAc                                                                                                          >  1:87601/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCACCACCCGGCAACCGACCACGTCGGCAGCCAGATAGCCACGCCCCCGCCCCCCCCCGCCCCCCTCCCCGGCCCACGCGGccccc                                                                                                          >  6:358172/1‑132 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGCCCACGTCGGCTGCCGGATAGCCACCCACCCGCACCGCCCCCCCCCCGGCCCCGGCCACCCCGGccc                                                                                                            >  7:213815/1‑132 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCACCCGCCCAGGTCGGCTGCCGGATAGCCCCGCGCCAGCCCCTCACCGCCCCCTTCCCGGGCCCCCGCGGCCAAc                                                                                                          >  8:220819/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCACACGACCAGGTCGGCTTCCAGATATACCCGCACCCGGCCCGCCCCCCCACCCTCCCGGGCCCACCCGGCCCAc                                                                                                          >  4:397989/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGCCTGATAGCCCCGCACCAGCACCTCACCGCCACCGTCCCCGGCCAACGCGGCCAAc                                                                                                          >  8:317500/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCCCCCGCCCCTCCCCGCCCCCGTCCCGGGCCCACCCGGCCCAc                                                                                                          >  1:82169/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCCCCGCCACCTTCCCGGGCCCACGCGGCCCAc                                                                                                          >  8:352534/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAAc                                                                                                          >  6:112945/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCACGTCGGCTGCCAGATACCCACGCACCAGCACCACACCGCACCCATCCCCGGCCAACCCGGCCCAc                                                                                                          >  6:127414/1‑135 (MQ=255)
                           gcATCGGCCGACACCTGGGCATTCAGCCAAGGGTCGTCAGTCGCCGAGATAACCAGCCGGCCCACGACCAGGTCGGCCGCGTGATCGCCACCCCCCAGCCCCTCACCGACACCTTCCCGGCCCCACTCGGCCAAc                                                                                                          >  8:125396/1‑135 (MQ=255)
                            cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGCCCAGGTCGGCTGCCGGATTGCCCCCCCCCCGCACCTCACCCCCCCCCTCCCGGGCCCCCCCGGCCCACt                                                                                                         >  5:342033/1‑135 (MQ=255)
                            cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGCCCAGGTCGGCTGCCGGATTGCCACGCCCCAGCACCTCCCCGCCCCCTTCCCGGGCCACCGCGGCCAACt                                                                                                         >  7:437109/1‑135 (MQ=255)
                             aTCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCCCCTTCCCGGGCCCACGCGGCCAACTg                                                                                                        >  3:215495/1‑135 (MQ=255)
                                           gggCATTCCGCCCTGGGTCGTCAGGCGCCGCGAACACCCGCCGCCACCCGACCAGGTCGGCTGCCTTATAGCCACCCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCg                                                                                          <  1:277673/135‑1 (MQ=255)
                                              cATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCCCCAGCCGGCCACCGCCCCGGTCGGCCGCCTGCTAGCCCCGCCCCCGCCCCTCACCGCCCCCGTCCCGGGCCCACGCCGCCCACAGGCCGGCGCCGGCCGgcg                                                                                       >  2:383367/1‑135 (MQ=255)
                                               aTTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCGGATAGCCACGCCCCAGCACCTCCCCGCCCCCTTCCCGGGCCACCGCGGCCCACTGGCCGTCGACGTCCGgcgc                                                                                      >  7:9001/1‑135 (MQ=255)
                                                    gCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCAACCGACCAGGTCGGCGGCCTGATAGCCACGCCCCAGCCCCTCACCGCCACCTTGCCGGGCCACAGCGGCCAACTGGCCGTCGACGTCCGGCGccac                                                                                   >  8:227983/1‑133 (MQ=255)
                                                    gCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGACCAGGTCGGCTGCCTGCTAGCCCCGCACCAGCCCCTCCCCGCCACCTTCCCGGGCCACCCCGGCCAACTGGGCGTCGCCGTCCGGCGccccca                                                                                 >  7:116916/1‑135 (MQ=255)
                                                           gggggcAGGCGCCGCTAACCCCAGCCGGCACCCGCCCCGGGCGGGGGCCTGATCCCCACCCACCCGCCCCACCCCGCCACCTCCCCGGGCCAACGCGGCCAACTGGCCGGCGACGTCCGGCGCCACCACGCgcag                                                                          <  2:243545/130‑1 (MQ=255)
                                                                  gTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCAcgccc                                                                   >  7:379180/1‑134 (MQ=255)
                                                                  gTCGCCGCCATCACCAGCCGCCAACCCACCAGGTCGGCTGGCTCATACCCCCGCACCACCACCTCCCCCGAACCTTCCCCGGCCACCGCGGCCACCTTGCCGTCCACGTCCGGCGCCACCACGcgccgccccccc                                                                   >  2:341347/1‑129 (MQ=255)
                                                                     gccgcGACCCCCCGCCGCCACCCGACCCGGGCGGCGGCCTGATCGCCACCCCCCCGCACCCCCCCCCCACCTCCCCGGCCCACCGCGCCCAACTGGCCGGCGACGTCCGGCGCCACCACGCGCAGCACGccgccg                                                                <  1:383367/135‑1 (MQ=255)
                                                                      ccgcGATCACCAGCCGGCAACCGACCCGGTCGGCTGCCTGATAGCCACCCACCAGCACCTCACCGCCACCTACCCGGGCCACCGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCgg                                                               <  8:394694/135‑1 (MQ=255)
                                                                        gcgATCACCAGCCGGCAACCGACCTGGGCGGCTGCGTGATAGCCACGCACCACCACCGCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGGCGACGGCCGGCGCCACCACGCGCAGCACGCCGCCGGCa                                                             <  7:317500/135‑1 (MQ=255)
                                                                         cgATCACCAGCCGGTAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCAt                                                            >  4:340182/1‑135 (MQ=255)
                                                                             caccccccGCCACCCGCCCAGGGCGGCTGCCTGATCCCCCCCCACCCCCACCTCCCCCCCACCTCCCCGGGCCACCGCGGCCAACTGGCCGGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGc                                                        <  2:82169/129‑1 (MQ=255)
                                                                             cACCAGCCGGCAACCGACCAGGGCGGCGGCCTGATAGCCACGCACCCGCACCGCACCCCCACCTACCCGGCCCACCGCGGCCAAGTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGc                                                        <  8:77041/135‑1 (MQ=255)
                                                                                            gACCCGGTCGGCTGCCTGATAGCCCCGCACCACCACCTCCCCCCCACCTCCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCtt                                         <  7:453743/135‑1 (MQ=255)
                                                                                             aCCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCCCCTTCCCGGGCCACCGCGGCCAACTGGCCGTCGACGTCCGGCGCCCCCACCCGCCGCCCGCCCCCGGCATCCGCCCGCCGGCGCGccggg                                        >  1:403740/1‑132 (MQ=255)
                                                                                                aGGTCGGCTGCCTGATAGCCACGCACCAGCCCCTCACCGCCCCCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCGTgcg                                      >  4:411594/1‑134 (MQ=255)
                                                                                                aGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCGCCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTgcg                                      >  8:334981/1‑134 (MQ=255)
                                                                                                   gcGGCTGCCGCATAGCCCCGCCCCAGCACCTCACCGCCACCTGCCCGGGCCACCGCGGCCAACTGGCCGGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAg                                  <  5:54356/134‑1 (MQ=255)
                                                                                                               tAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCATGCGCAAGGCGCTCTCACCg                      >  8:263381/1‑135 (MQ=255)
                                                                                                               tAGCCACGCACCAGCACCTCACCGCCACCTTCACGGGCCAAAGCGGCCAACTCGCCGTCGACCTCCGGCGCCCCCCCGCGCAGCACGCCGCACGCATCGGCCAGAAGGCCCCCCTAGCGCCACGCGAtcaacccc                      >  2:420172/1‑129 (MQ=255)
                                                                                                               tAGCCACGCACCAGCACCTCACCGCCACCTACCCGGGCCAACGCGGCCAACTGGCCGTCGACCTCCGGCGCCCCCACGCGCAGCACGCCGCCGGCCTCGGCCAGCAGGCCCGCCTTGCGCCACGCGCTCACCCCg                      >  2:420169/1‑135 (MQ=255)
                                                                                                                       cACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCGTGCGCAAGGCGCTCTCACCGCccccccc              >  1:414501/1‑132 (MQ=255)
                                                                                                                                   aCCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCCCCCCCGACCAGCAcccg  >  3:179748/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                      
ACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCG  >  ACB122/4500090‑4500355

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