Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,911,836 Δ1 bp 12.2% coding (390/477 nt) ccmH ← cytochrome c heme lyase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,911,8360A.12.2% 116.3 / 11.8 41coding (390/477 nt)ccmHcytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base . (2/3);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

GGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGAT  >  NC_002947/4911687‑4911982
                                                                                                                                                      |                                                                                                                                                  
ggATCTGGAGCGGCAAAAAGCTGAGGGCAGCGAGCAGAAGAAGGCCCGCGCTAAGCGAGATTTCAATCAGGGGTCTGTTCTTTTTCCAGCTGTTGGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTa                                                                                                                                                    <  2:269265/150‑1 (MQ=255)
    cGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCAc                                                                                                                                                <  1:430973/150‑1 (MQ=255)
         gcagcaAAAAGCTGAGGGCAGCGAGCTGCAGGAGGCCTTCGCTAAGCCAGAATTCAATCATGGGTTCATAATCGCTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGgcgc                                                                                                                                           <  1:1040638/150‑1 (MQ=255)
         gcagcaAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGgcgc                                                                                                                                           >  1:905067/1‑150 (MQ=255)
              aaaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGCCgatg                                                                                                                                      >  2:1109405/1‑147 (MQ=255)
               aaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGCCgatga                                                                                                                                     <  1:1109405/150‑5 (MQ=255)
                 aaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACATATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTAc                                                                                                                                   >  1:289659/1‑150 (MQ=255)
                    cTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGAt                                                                                                                                <  1:771936/150‑1 (MQ=255)
                                       gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCAc                                                                                                                                                >  2:420008/1‑115 (MQ=255)
                                       gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCAc                                                                                                                                                <  1:420008/115‑1 (MQ=255)
                                                        cAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGAt                                                                                            >  2:890048/1‑150 (MQ=255)
                                                             ttCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCAAAGCTCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                                       <  1:850972/150‑1 (MQ=255)
                                                                aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                                    >  2:890616/1‑150 (MQ=255)
                                                                aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                                    <  1:696581/150‑1 (MQ=255)
                                                                         tCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcaca                                                                           >  2:900458/1‑150 (MQ=255)
                                                                             ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagcag                                                                       >  1:934882/1‑150 (MQ=255)
                                                                                 ttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAg                                                                   <  1:890616/150‑1 (MQ=255)
                                                                                 ttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAg                                                                   <  1:1059301/150‑1 (MQ=255)
                                                                                       agcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGc                                                             >  2:231989/1‑150 (MQ=255)
                                                                                                 cGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAgcg                                                   >  1:546204/1‑150 (MQ=255)
                                                                                                  gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGc                                                                                                                      >  1:804340/1‑82 (MQ=255)
                                                                                                  gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGc                                                                                                                      <  2:804340/82‑1 (MQ=255)
                                                                                                   agaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGTgcggt                                                 <  2:546204/150‑3 (MQ=255)
                                                                                                    gaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCg                                                >  1:962995/1‑150 (MQ=255)
                                                                                                           gcgTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                                         >  1:663461/1‑150 (MQ=255)
                                                                                                                    ggCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                                >  2:320303/1‑150 (MQ=255)
                                                                                                                     gCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAg                               >  2:178927/1‑150 (MQ=255)
                                                                                                                      ctgACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGt                              <  1:308053/148‑1 (MQ=255)
                                                                                                                         aCAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTCGCCCACTaatgcag                           >  2:308053/1‑145 (MQ=255)
                                                                                                                           aaCTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                                         >  1:458744/1‑134 (MQ=255)
                                                                                                                           aaCTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                                         <  2:458744/134‑1 (MQ=255)
                                                                                                                             cTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                                       <  1:753877/86‑1 (MQ=255)
                                                                                                                             cTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                                       >  2:753877/1‑86 (MQ=255)
                                                                                                                              tCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg                      >  1:338553/1‑150 (MQ=255)
                                                                                                                                  cGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt                  <  1:681888/150‑1 (MQ=255)
                                                                                                                                   ggCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTc                 >  1:605393/1‑150 (MQ=255)
                                                                                                                                    gcgcCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCg                <  2:962995/150‑1 (MQ=255)
                                                                                                                                       ccgtgctaCGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCg                >  2:597737/8‑147 (MQ=255)
                                                                                                                                       ccgtgctaCGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCg                <  1:597737/140‑1 (MQ=255)
                                                                                                                                        ccTGCACGGCCGTAGCACGGCGCCGGCGTACGAACACCGCCGGCACCACGAAGCCACCGGCCAGCAAGATCCGCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAccc            >  1:933388/1‑149 (MQ=255)
                                                                                                                                                   cGTAGCCCGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCgat  <  1:490197/150‑1 (MQ=255)
                                                                                                                                                   cGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGTCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCgat  <  2:605393/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GGATCGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGAT  >  NC_002947/4911687‑4911982

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: