Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,293,137 C→T 9.7% A590T (GCC→ACC)  ccmF ← heme lyase, CcmF subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,293,1370CT 9.7% 117.4 / 3.3 41A590T (GCC→ACC) ccmFheme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (19/18);  new base T (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

TCACGCGCTTACGATAGCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGCCGATAGTCGCCACACCGC  >  NC_000913/2292993‑2293277
                                                                                                                                                |                                                                                                                                            
tCACGCGCTTACGATAGCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCcgcg                                                                                                                                          >  4:54629/1‑149 (MQ=255)
  aCGCGCTTACGATAGCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGt                                                                                                                                         <  8:17097/148‑1 (MQ=255)
  aCGCGCTTACGATAGCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTa                                                                                                                                        >  4:190567/1‑149 (MQ=255)
  aCGCGCTTACGATAGCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTa                                                                                                                                        >  8:219272/1‑149 (MQ=255)
    gcgcTTACGATAGCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCCGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCTCGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACa                                                                                                                                      <  6:95802/149‑1 (MQ=255)
                gCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCCCGCGTGa                                                                                                                          >  5:152825/1‑149 (MQ=255)
                         aaaCAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGt                                                                                                                 <  7:219272/149‑1 (MQ=255)
                            cAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCaa                                                                                                              <  6:341473/149‑1 (MQ=255)
                             aGACACAGTAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCAACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAAt                                                                                                             >  2:97924/1‑149 (MQ=255)
                               acacAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTg                                                                                                           >  8:321931/1‑149 (MQ=255)
                                          cGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGtc                                                                                                <  7:241586/149‑1 (MQ=255)
                                                              ccgcccgcccAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGGGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGtgtt                                                                            <  8:334966/149‑1 (MQ=255)
                                                                        aaaTCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGt                                                                  >  7:120246/1‑149 (MQ=255)
                                                                          aTCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGttt                                                                <  3:190567/149‑1 (MQ=255)
                                                                             cAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTc                                                             <  3:100518/149‑1 (MQ=255)
                                                                               gCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAg                                                                                >  1:371625/1‑128 (MQ=255)
                                                                               gCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAg                                                                                <  2:371625/128‑1 (MQ=255)
                                                                               gCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCg                                                           >  3:225000/1‑149 (MQ=255)
                                                                                   aCAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGATGGCCGCGTACAGGTCACGAGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATa                                                       <  8:85634/149‑1 (MQ=255)
                                                                                   aCAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATa                                                       >  2:35868/1‑149 (MQ=255)
                                                                                      aaTGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAg                                                    >  3:14193/1‑149 (MQ=255)
                                                                                       aTGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGc                                                   >  8:13534/1‑149 (MQ=255)
                                                                                                   aaaGCCGCCCGGCCCCCGCGCCGTTTTCCAGCTCTGCACCGCGGGCCGCGTACAGGTCCCGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATCACGTTTTTCCGCATACAGCACCGTTTCCGGc                                       <  4:68307/149‑1 (MQ=255)
                                                                                                   aaaGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGc                                       >  3:305271/1‑149 (MQ=255)
                                                                                                    aaGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTcatcat                                                                                           <  2:207309/96‑1 (MQ=255)
                                                                                                    aaGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTcatcat                                                                                           >  1:207309/1‑96 (MQ=255)
                                                                                                      gACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCg                                         >  7:203873/1‑144 (MQ=255)
                                                                                                      gACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCg                                         <  8:203873/144‑1 (MQ=255)
                                                                                                               cccACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATcgcgc                           <  4:14193/149‑1 (MQ=255)
                                                                                                                   cgcgCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTg                                                                                                           >  2:369021/1‑65 (MQ=255)
                                                                                                                   cgcgCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTg                                                                                                           <  1:369021/65‑1 (MQ=255)
                                                                                                                      gcCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGc                    >  3:199338/1‑149 (MQ=255)
                                                                                                                       ccGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGcc                   <  1:97924/149‑1 (MQ=255)
                                                                                                                        cGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGcc                   >  5:299825/1‑148 (MQ=255)
                                                                                                                                aGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGCCGATAGTCGc          <  3:160422/149‑1 (MQ=255)
                                                                                                                                aGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGCCGATAGTCGc          <  1:35868/149‑1 (MQ=255)
                                                                                                                                    ctTCACCGAGGGTCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGc                                                          >  4:230306/1‑97 (MQ=255)
                                                                                                                                    ctTCACCGAGGGTCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGc                                                          >  4:225185/1‑97 (MQ=255)
                                                                                                                                    ctTCACCGAGGGTCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGc                                                          <  3:230306/97‑1 (MQ=255)
                                                                                                                                    ctTCACCGAGGGTCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGc                                                          <  3:225185/97‑1 (MQ=255)
                                                                                                                                        aCCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGCCGATAGTCGCCACACcgc  >  3:235579/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                            
TCACGCGCTTACGATAGCGAGGATCAAACAGACACAGCAGTCCGCCCAACGCCATCATCAGCCCGCCCGCCCAAATCCAGCGAACAAATGGTTTGTAGTAAAGACGCACGGCCCACGCGCCGTTTTCCAGCTCTTCACCGAGGGCCGCGTACAGGTCACGCGTGATGCCGCCGTCAATTGCCGCTTCGGTCATCATCGACCCGGCAGTGTTGTAATAACGTTTTTCCGCATACAGCACCGTTTCCGGCTTGCCATCGCGCGTTACGCCGATAGTCGCCACACCGC  >  NC_000913/2292993‑2293277

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