Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,500,241 A→C 32.0% L52V (TTG→GTG)  cysG_1 ← Siroheme synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,500,2410AC32.0% 49.0 / 3.2 51L52V (TTG→GTG) cysG_1Siroheme synthase
Reads supporting (aligned to +/- strand):  ref base A (26/7);  new base C (10/8);  total (36/15)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTG  >  ACB122/4500117‑4500373
                                                                                                                            |                                                                                                                                     
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)
                gggCATTCCGCCCTGGGTCGTCAGGCGCCGCGAACACCCGCCGCCACCCGACCAGGTCGGCTGCCTTATAGCCACCCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCg                                                                                                            <  1:277673/135‑1 (MQ=255)
                   cATTCAGCCCAGGGTCGTCAGTCGCCGCGCTCCCCAGCCGGCCACCGCCCCGGTCGGCCGCCTGCTAGCCCCGCCCCCGCCCCTCACCGCCCCCGTCCCGGGCCCACGCCGCCCACAGGCC‑GGCGCCGGCCGgcg                                                                                                         >  2:383367/1‑135 (MQ=255)
                    aTTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCGGATAGCCACGCCCCAGCACCTCCCCGCCCCCTTCCCGGGCCACCGCGGCCCACTGGCC‑GTCGACGTCCGgcgc                                                                                                        >  7:9001/1‑135 (MQ=255)
                         gCCCAGGGTCGTCAGTCGCCGCGCTCACCAGCCGGCAACCGACCAGGTCGGCGGCCTGATAGCCACGCCCCAGCCCCTCACCGCCACCTTGCCGGGCCACAGCGGCCAACTGGCC‑GTCGACGTCCGGCGccac                                                                                                     >  8:227983/1‑133 (MQ=255)
                         gCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCCACCGACCAGGTCGGCTGCCTGCTAGCCCCGCACCAGCCCCTCCCCGCCACCTTCCCGGGCCACCCCGGCCAACTGGGC‑GTCGCCGTCCGGCGccccca                                                                                                   >  7:116916/1‑135 (MQ=255)
                                gggggcAGGCGCCGCTAACCCCAGCCGGCACCCGCCCCGGGCGGGGGCCTGATCCCCACCCACCCGCCCCACCCCGCCACCTCCCCGGGCCAACGCGGCCAACTGGCC‑GGCGACGTCCGGCGCCACCACGCgcag                                                                                            <  2:243545/130‑1 (MQ=255)
                                       gTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCAcgccc                                                                                     >  7:379180/1‑134 (MQ=255)
                                       gTCGCCGCCATCACCAGCCGCCAACCCACCAGGTCGGCTGGCTCATACCCCCGCACCACCACCTCCCCCGAACCTTCCCCGGCCACCGCGGCCACCTTGCC‑GTCCACGTCCGGCGCCACCACGcgccgccccccc                                                                                     >  2:341347/1‑129 (MQ=255)
                                          gccgcGACCCCCCGCCGCCACCCGACCCGGGCGGCGGCCTGATCGCCACCCCCCCGCACCCCCCCCCCACCTCCCCGGCCCACCGCGCCCAACTGGCC‑GGCGACGTCCGGCGCCACCACGCGCAGCACGccgccg                                                                                  <  1:383367/135‑1 (MQ=255)
                                           ccgcGATCACCAGCCGGCAACCGACCCGGTCGGCTGCCTGATAGCCACCCACCAGCACCTCACCGCCACCTACCCGGGCCACCGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCgg                                                                                 <  8:394694/135‑1 (MQ=255)
                                             gcgATCACCAGCCGGCAACCGACCTGGGCGGCTGCGTGATAGCCACGCACCACCACCGCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GGCGACGGCCGGCGCCACCACGCGCAGCACGCCGCCGGCa                                                                               <  7:317500/135‑1 (MQ=255)
                                              cgATCACCAGCCGGTAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCAt                                                                              >  4:340182/1‑135 (MQ=255)
                                                  caccccccGCCACCCGCCCAGGGCGGCTGCCTGATCCCCCCCCACCCCCACCTCCCCCCCACCTCCCCGGGCCACCGCGGCCAACTGGCC‑GGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGc                                                                          <  2:82169/129‑1 (MQ=255)
                                                  cACCAGCCGGCAACCGACCAGGGCGGCGGCCTGATAGCCACGCACCCGCACCGCACCCCCACCTACCCGGCCCACCGCGGCCAAGTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGc                                                                          <  8:77041/135‑1 (MQ=255)
                                                                 gACCCGGTCGGCTGCCTGATAGCCCCGCACCACCACCTCCCCCCCACCTCCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCtt                                                           <  7:453743/135‑1 (MQ=255)
                                                                  aCCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCCCCTTCCCGGGCCACCGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCCCCACCCGCCGCCCGCCCCCGGCATCCGCCCGCCGGCGCGccggg                                                          >  1:403740/1‑132 (MQ=255)
                                                                     aGGTCGGCTGCCTGATAGCCACGCACCAGCCCCTCACCGCCCCCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCGTgcg                                                        >  4:411594/1‑134 (MQ=255)
                                                                     aGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCGCCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTgcg                                                        >  8:334981/1‑134 (MQ=255)
                                                                        gcGGCTGCCGCATAGCCCCGCCCCAGCACCTCACCGCCACCTGCCCGGGCCACCGCGGCCAACTGGCC‑GGCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAg                                                    <  5:54356/134‑1 (MQ=255)
                                                                                    tAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCATGCGCAAGGCGCTCTCACCg                                        >  8:263381/1‑135 (MQ=255)
                                                                                    tAGCCACGCACCAGCACCTCACCGCCACCTTCACGGGCCAAAGCGGCCAACTCGCC‑GTCGACCTCCGGCGCCCCCCCGCGCAGCACGCCGCACGCATCGGCCAGAAGGCCCCCCTAGCGCCACGCGAtcaacccc                                        >  2:420172/1‑129 (MQ=255)
                                                                                    tAGCCACGCACCAGCACCTCACCGCCACCTACCCGGGCCAACGCGGCCAACTGGCC‑GTCGACCTCCGGCGCCCCCACGCGCAGCACGCCGCCGGCCTCGGCCAGCAGGCCCGCCTTGCGCCACGCGCTCACCCCg                                        >  2:420169/1‑135 (MQ=255)
                                                                                            cACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCGTGCGCAAGGCGCTCTCACCGCccccccc                                >  1:414501/1‑132 (MQ=255)
                                                                                                      tccccGCCACCTTCCCGGGCCAACGGGGCAAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCAcg                      <  7:348533/132‑1 (MQ=255)
                                                                                                      gccccGCCACCTTCCCGGGCCACCGCGGCCAACTGGCC‑GTCGACGGCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCAcg                      <  4:285684/132‑1 (MQ=255)
                                                                                                        aCCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCCCCCCCGACCAGCAcccg                    >  3:179748/1‑135 (MQ=255)
                                                                                                          cGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCCCCGACGACCAGCACGCGgc                  >  3:269537/1‑135 (MQ=255)
                                                                                                          cGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGgc                  >  2:88677/1‑135 (MQ=255)
                                                                                                              cccctCCCGGGCCCACGGGGCCCACGGGCCGGGCGACGCCCGGCGCCACCCCGCGCAGCCCGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGgccgc               <  5:463405/131‑1 (MQ=255)
                                                                                                              aCCTTCCCGGGCCACCGCGGCAAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGgccgcc              <  8:23591/135‑1 (MQ=255)
                                                                                                                  tCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCCCCGACGACCCGCCCGCGGCCGCCCt            >  3:393879/1‑133 (MQ=255)
                                                                                                                    ccGGGCAACCGCGGCCAACTGGCG‑GGCGACGGCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGc          <  8:113871/133‑1 (MQ=255)
                                                                                                                     cGGGCCAACGGGGCCAAGTGGCG‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGc       <  4:215495/135‑1 (MQ=255)
                                                                                                                     cGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGCTCTCACCGCCCCCGACGACCAGCACGCGGCCGCCCTGCAGc       >  6:465225/1‑135 (MQ=255)
                                                                                                                          cAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGCTCTCACCGCCCCCGCCGCCCAGCACGCGGCCGCCCTGCAGCGtgtg  >  4:219390/1‑135 (MQ=255)
                                                                                                                            |                                                                                                                                     
GCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTG  >  ACB122/4500117‑4500373

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: ‑