Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,291,973 T→C G146G (GGA→GGG ccmH ← heme lyase, CcmH subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,291,9730TC96.7% 207.6 / ‑0.3 60coding (438/1053 nt)ccmHheme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (26/32);  minor base . (1/1);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.33e-01

CGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCACG  >  NC_000913/2291838‑2292108
                                                                                                                                       |                                                                                                                                       
cgcCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCg                                                                                                                                      <  2:22468/139‑1 (MQ=255)
       tCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCt                                                                                                                               <  1:777094/139‑1 (MQ=255)
                 ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGATGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCAGCACGCTTACCTTCCgg                                                                                                                     >  2:670281/1‑139 (MQ=255)
                  cGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                                    <  2:1364279/139‑1 (MQ=255)
                           tgCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGt                                                                                                           >  1:430894/1‑139 (MQ=255)
                           tgCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGt                                                                                                           >  1:1386509/1‑139 (MQ=255)
                            gCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGtt                                                                                                          <  2:1526816/139‑1 (MQ=255)
                                  gCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaa                                                                                                    >  2:136017/1‑139 (MQ=255)
                                  gCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCCCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaaa                                                                                                   >  1:697668/1‑139 (MQ=255)
                                     tgctgcCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACg                                                                                                 >  2:18205/1‑139 (MQ=255)
                                               ttttCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCtt                                                                                                                              <  2:765399/100‑1 (MQ=255)
                                               ttttCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCtt                                                                                                                              >  1:765399/1‑100 (MQ=255)
                                                            aTAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                                       >  2:702865/1‑126 (MQ=255)
                                                            aTAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                                       <  1:702865/126‑1 (MQ=255)
                                                                     aGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                 <  1:18205/139‑1 (MQ=255)
                                                                     aGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                 <  1:670281/139‑1 (MQ=255)
                                                                         tGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                             <  1:136017/139‑1 (MQ=255)
                                                                         tGGTAGTAGCTGACGCCAGCCACAATTAACGCCCCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGc                                                            <  2:697668/139‑1 (MQ=255)
                                                                            tagtagCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcaccac                                                                                    >  2:1475556/1‑113 (MQ=255)
                                                                            tagtagCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcaccac                                                                                    <  1:1475556/113‑1 (MQ=255)
                                                                                 gCTGACGCCAGCCACAATTAATGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                             >  2:1154604/1‑131 (MQ=255)
                                                                                 gCTGACGCCAGCCACAATTAATGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                             <  1:1154604/131‑1 (MQ=255)
                                                                                 gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTa                                                                                                                             >  1:747445/1‑67 (MQ=255)
                                                                                 gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTa                                                                                                                             <  2:747445/67‑1 (MQ=255)
                                                                                 gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                               >  1:1457649/1‑81 (MQ=255)
                                                                                 gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                               <  2:1457649/81‑1 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCt                                                                                                                               <  1:59046/61‑1 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCt                                                                                                                               >  2:59046/1‑61 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaaa                                                                                                   >  2:408016/1‑89 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaaa                                                                                                   <  1:408016/89‑1 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGca                                                                                        >  1:1125843/1‑100 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGca                                                                                        <  2:1125843/100‑1 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                 >  1:172135/1‑123 (MQ=255)
                                                                                     aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                 <  2:172135/123‑1 (MQ=255)
                                                                                            ccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCGGCACGCTTACCtt                                                                                                                         >  2:932470/1‑59 (MQ=255)
                                                                                            ccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCGGCACGCTTACCtt                                                                                                                         <  1:932470/59‑1 (MQ=255)
                                                                                                   tAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcca                                   <  2:1179135/139‑1 (MQ=255)
                                                                                                   tAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTCCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcca                                   <  2:925743/139‑1 (MQ=255)
                                                                                                      cGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATa                                >  1:923004/1‑139 (MQ=255)
                                                                                                           ccacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcca                                   <  1:223247/131‑1 (MQ=255)
                                                                                                           ccacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcca                                   >  2:223247/1‑131 (MQ=255)
                                                                                                                   cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGca                                                                                        <  2:566175/70‑1 (MQ=255)
                                                                                                                   cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGca                                                                                        >  1:566175/1‑70 (MQ=255)
                                                                                                                   cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGgca                   <  2:923004/139‑1 (MQ=255)
                                                                                                                   cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGgca                   <  2:430894/139‑1 (MQ=255)
                                                                                                                   cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGgca                   <  1:867142/139‑1 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                                    >  2:1274508/1‑39 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                                    <  1:1274508/39‑1 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCcgcg                                                                   <  2:214120/88‑1 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCcgcg                                                                   >  1:214120/1‑88 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGat                                                                 >  1:172435/1‑89 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                 <  2:172435/90‑1 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCgg                                                              <  2:751084/93‑1 (MQ=255)
                                                                                                                      ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCgg                                                              >  1:751084/1‑93 (MQ=255)
                                                                                                                       gCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCGGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                 <  1:1297112/89‑1 (MQ=255)
                                                                                                                       gCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCGGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                                 >  2:1297112/1‑89 (MQ=255)
                                                                                                                         aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                             <  1:126137/91‑1 (MQ=255)
                                                                                                                         aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                             >  2:126137/1‑91 (MQ=255)
                                                                                                                             aaaCAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCAcag         >  2:112974/1‑139 (MQ=255)
                                                                                                                                caacaTACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCAcagcag      <  2:786065/139‑1 (MQ=255)
                                                                                                                                    atacCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCACg  >  2:692109/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGCAGCACG  >  NC_000913/2291838‑2292108

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