Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,917,404 A→T 9.3% F25I (TTC→ATC)  ccmA ← cytochrome c biogenesis ABC transporter ATP‑binding protein CcmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,917,4040AT 9.3% 128.1 / 4.6 43F25I (TTC→ATC) ccmAcytochrome c biogenesis ABC transporter ATP‑binding protein CcmA
Reads supporting (aligned to +/- strand):  ref base A (17/22);  new base T (2/2);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAA  >  NC_002947/4917256‑4917542
                                                                                                                                                    |                                                                                                                                          
gCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGaaa                                                                                                                                           >  2:72999/1‑150 (MQ=255)
         tGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACATGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGATTTCGAGCTgg                                                                                                                                  <  2:688955/150‑2 (MQ=255)
         tGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGATTTCGAGCTgg                                                                                                                                  <  2:583390/150‑2 (MQ=255)
                ccAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAAca                                                                                                                           >  2:740434/1‑150 (MQ=255)
                        ttGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGgcgc                                                                                                                   <  1:389467/150‑1 (MQ=255)
                           ccACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAg                                                                                                                >  2:575196/1‑150 (MQ=255)
                           ccACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAg                                                                                                                <  2:427951/150‑1 (MQ=255)
                             aCCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTc                                                                                                              <  1:16005/150‑1 (MQ=255)
                                       aTCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAAt                                                                                                                                          >  2:768513/1‑112 (MQ=255)
                                       aTCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAAt                                                                                                                                          <  1:768513/112‑1 (MQ=255)
                                           gCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggcca                                                                                                <  2:517066/150‑1 (MQ=255)
                                           gCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATATGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggcca                                                                                                <  2:378338/150‑1 (MQ=255)
                                             ccGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTc                                                                                                              <  2:754184/134‑1 (MQ=255)
                                             ccGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTc                                                                                                              >  1:754184/1‑134 (MQ=255)
                                             ccGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggccagg                                                                                              >  1:541631/1‑150 (MQ=255)
                                                  cGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCcgc                                                                                         >  1:375719/1‑150 (MQ=255)
                                                   gTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCcgcg                                                                                        >  2:260924/1‑150 (MQ=255)
                                                               aGCCCGGCCAGCAGGCGGAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTc                                                                                                              >  2:408156/1‑116 (MQ=255)
                                                               aGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTc                                                                                                              <  1:408156/116‑1 (MQ=255)
                                                                  ccGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGATTGGAGGTGAAgg                                                                         <  1:72999/150‑1 (MQ=255)
                                                                        agcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACg                                                                   >  2:548056/1‑150 (MQ=255)
                                                                                cAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGcc                                                           <  2:541631/150‑1 (MQ=255)
                                                                                       cTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCtt                                                    <  1:323474/150‑1 (MQ=255)
                                                                                            cTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTa                                               >  2:482870/1‑150 (MQ=255)
                                                                                                    tGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTg                                       <  2:186998/150‑1 (MQ=255)
                                                                                                              gCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCaa                             <  2:497676/150‑1 (MQ=255)
                                                                                                                cACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAaca                           <  1:740434/150‑1 (MQ=255)
                                                                                                                 aCTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAgg                                                                         <  2:521955/103‑1 (MQ=255)
                                                                                                                 aCTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAgg                                                                         >  1:521955/1‑103 (MQ=255)
                                                                                                                   tGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCg                        >  1:453237/1‑150 (MQ=255)
                                                                                                                         gcagcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagc                                   <  2:109256/133‑1 (MQ=255)
                                                                                                                         gcagcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagc                                   >  1:109256/1‑133 (MQ=255)
                                                                                                                         gcagcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAg                  >  2:211332/1‑150 (MQ=255)
                                                                                                                          cagcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGa                 >  1:558725/1‑150 (MQ=255)
                                                                                                                            gcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagcagc                                >  2:537483/1‑133 (MQ=255)
                                                                                                                            gcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagcagc                                <  1:537483/133‑1 (MQ=255)
                                                                                                                                  cGCCGGCGCCCAGCTCGAAATCGGGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtgt         <  2:435354/150‑1 (MQ=255)
                                                                                                                                  cGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtgt         <  1:482870/150‑1 (MQ=255)
                                                                                                                                  cGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtgt         <  1:260924/150‑1 (MQ=255)
                                                                                                                                  cGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtgt         <  1:211332/150‑1 (MQ=255)
                                                                                                                                         gaacAGCTCGATTTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCaa  >  2:176328/4‑150 (MQ=255)
                                                                                                                                         gaacAGCTCGATTTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCaa  >  2:13138/4‑150 (MQ=255)
                                                                                                                                                tCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCAc                                                                    <  2:493168/77‑1 (MQ=255)
                                                                                                                                                tCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCAc                                                                    >  1:493168/1‑77 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAA  >  NC_002947/4917256‑4917542

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