Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,295,693 C→T 22.2% A18T (GCG→ACG)  ccmH ← heme lyase, CcmH subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,295,6930CT22.2% 74.3 / 14.6 36A18T (GCG→ACG) ccmHheme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (16/12);  new base T (4/4);  total (20/16)
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

TTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTT  >  W3110S.gb/2295625‑2295758
                                                                    |                                                                 
ttcttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGc                                                                   <  1:318095/69‑1 (MQ=255)
ttcttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGc                                                                   <  1:2770343/69‑1 (MQ=255)
  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCa                                                                 >  1:2648637/1‑69 (MQ=255)
     cAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcg                                                              >  1:3163609/1‑69 (MQ=255)
         gAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGa                                                          <  1:1601982/69‑1 (MQ=255)
                    aCTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatca                                               <  1:950988/69‑1 (MQ=255)
                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatcatc                                             >  1:1019865/1‑68 (MQ=255)
                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatcatc                                             >  1:1715012/1‑68 (MQ=255)
                          gttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcag                                           <  1:2689098/67‑1 (MQ=255)
                          gttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcag                                           <  1:2689096/67‑1 (MQ=255)
                          gttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcag                                           <  1:2180984/67‑1 (MQ=255)
                          gttCCTGTGCTTCATCTTTAAACCGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcag                                           <  1:2700912/67‑1 (MQ=255)
                             ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatc                                       <  1:1009444/68‑1 (MQ=255)
                                 tgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCAc                                  >  1:1493497/1‑69 (MQ=255)
                                 tgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCAc                                  >  1:2541602/1‑69 (MQ=255)
                                 tgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCAc                                  >  1:1387904/1‑69 (MQ=255)
                                          tttAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATa                          >  1:2062633/1‑68 (MQ=255)
                                          tttAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagc                                          >  1:813929/1‑52 (MQ=255)
                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGc                                                        <  1:2933604/36‑1 (MQ=255)
                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGc                                                        <  1:2270221/36‑1 (MQ=255)
                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGc                                                        <  1:2265881/36‑1 (MQ=255)
                                                cTGCAAGACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCt                   >  1:8348/1‑69 (MQ=255)
                                                cTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCt                   <  1:470861/69‑1 (MQ=255)
                                                cTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAAc                     <  1:2617415/67‑1 (MQ=255)
                                                cTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCt                   >  1:130729/1‑69 (MQ=255)
                                                cTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCt                   >  1:2193733/1‑69 (MQ=255)
                                                cTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCt                   >  1:53923/1‑69 (MQ=255)
                                                cTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCt                   >  1:1416440/1‑69 (MQ=255)
                                                 tGCAACACATCGATGGTCGCCAGCGCTGAGGCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTc                  >  1:2230463/1‑69 (MQ=255)
                                                     acacATCGATGGTCGCCAGCGCTGAGCCGGAGatca                                               <  1:3275142/36‑1 (MQ=255)
                                                         aTCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCg          <  1:3461729/69‑1 (MQ=255)
                                                          tCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                     >  1:3581453/1‑41 (MQ=255)
                                                            gATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatcat                                     <  1:2077893/39‑2 (MQ=255)
                                                            gATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatcag                                     <  1:1200925/39‑1 (MQ=255)
                                                            gATGGTCGCCAGCGCTGAGCCGGAGATcatcaacatcag                                     >  1:3429258/1‑39 (MQ=255)
                                                              tGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc     >  1:286797/1‑69 (MQ=255)
                                                              tGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc     >  1:1661827/1‑69 (MQ=255)
                                                                 tCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCtt  >  1:370674/1‑69 (MQ=255)
                                                                    |                                                                 
TTCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTT  >  W3110S.gb/2295625‑2295758

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