Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,296,418 T→C T203A (ACC→GCC)  ccmB ← heme trafficking system membrane protein CcmB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,296,4180TC97.1% 115.6 / ‑4.7 34T203A (ACC→GCC) ccmBheme trafficking system membrane protein CcmB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (12/21);  total (12/22)
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

GCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAG  >  NC_000913/2296281‑2296544
                                                                                                                                         |                                                                                                                              
gCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGgcgc                                                                                                                              >  2:862157/1‑140 (MQ=255)
 ccGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGcgcc                                                                                                                             >  2:657209/1‑140 (MQ=255)
 ccGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGcgcc                                                                                                                             >  1:463516/1‑140 (MQ=255)
     ggtggGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGcgcctgcc                                                                                                                         >  2:891550/1‑140 (MQ=255)
           aTCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAAc                                                                                                                   >  2:768462/1‑140 (MQ=255)
            tCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAAcg                                                                                                                  <  1:768462/140‑1 (MQ=255)
                  gTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTa                                                                                                            <  1:809823/140‑1 (MQ=255)
                          aGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGcc                                                                                                    <  2:828777/140‑1 (MQ=255)
                                 tCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATAcc                                                                                             >  1:117626/1‑140 (MQ=255)
                                       tAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCaa                                                                                       <  1:509160/140‑1 (MQ=255)
                                          tttCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACgg                                                                                    >  2:64347/1‑140 (MQ=255)
                                            tCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGc                                                                                  <  2:716690/140‑1 (MQ=255)
                                             cGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCa                                                                                 <  1:362212/140‑1 (MQ=255)
                                             cGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCa                                                                                 <  1:906864/140‑1 (MQ=255)
                                               aTACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCaaa                                                                               <  2:463516/140‑1 (MQ=255)
                                                   cAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCa                                                                           >  2:32212/1‑140 (MQ=255)
                                                      aCTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAg                                                                        >  1:320963/1‑140 (MQ=255)
                                                             cAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCcgc                                                                 >  2:1289183/1‑140 (MQ=255)
                                                                        aaTCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCgcgg                                                      <  1:113591/140‑1 (MQ=255)
                                                                         aTCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCgcggc                                                     <  2:117626/140‑1 (MQ=255)
                                                                            cAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGt                                                  >  2:517410/1‑140 (MQ=255)
                                                                             aGCGTTATTGAATGCAGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTg                                                 <  2:1183551/140‑1 (MQ=255)
                                                                              gCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAGCGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTgg                                                <  2:1070753/140‑1 (MQ=255)
                                                                              gCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCgcggcg                                                    >  2:138455/1‑136 (MQ=255)
                                                                              gCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCgcggcg                                                    <  1:138455/136‑1 (MQ=255)
                                                                               cGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGc                                               >  2:26542/1‑140 (MQ=255)
                                                                               cGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGc                                               >  2:792534/1‑140 (MQ=255)
                                                                                gTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCa                                              <  1:26542/140‑1 (MQ=255)
                                                                                 ttATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCaa                                             >  2:475366/1‑140 (MQ=255)
                                                                                   aTTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAg                                           <  1:64347/140‑1 (MQ=255)
                                                                                      gAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATg                                        <  2:320963/140‑1 (MQ=255)
                                                                                               ttCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATggg                               <  1:475366/140‑1 (MQ=255)
                                                                                                        cTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGc                      <  1:1289183/140‑1 (MQ=255)
                                                                                                          gccgccGTCGCAAAAGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCgg                    >  2:478976/1‑140 (MQ=255)
                                                                                                                       aaGGACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATg       <  1:478976/140‑1 (MQ=255)
                                                                                                                         ggACTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCt     <  2:174901/140‑1 (MQ=255)
                                                                                                                            cTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGag  <  1:862157/140‑1 (MQ=255)
                                                                                                                            cTTAATGTCGCGGCGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGag  <  1:891550/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                              
GCCGTGGTGGGATCGCCAGTTGATGCAGTGTTTTCCACATAGTTTCGATACCAGACTCGAACAAAAATCAGTAATCCAGCGTTATTGAATGCTGATTCGTAACGCTGCCGCCGTCGCAAAAGGACTTAATGTCGCGGTGCCTGCCAGCAACGCGCCTAAAATTGCCAGATACCCGTCAACGGGCAAATGCATAGAAGCCGCGTCCATCGCGGCGGTGGCAAAGATGAGTAATGGGATAGTCAGCGGTAACACCAGTATGCTGAG  >  NC_000913/2296281‑2296544

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