Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,295,693 C→T 100% 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,6930CT80.0% 82.4 / 17.7 40A18T (GCG→ACG) ccmHheme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (7/1);  new base T (16/16);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCT  >  W3110S.gb/2295626‑2295757
                                                                   |                                                                
tcttcACTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAg                                                               >  1:1623/1‑71 (MQ=255)
 cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAgc                                                              >  1:301191/1‑71 (MQ=255)
                ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGat                                                >  1:2618876/1‑70 (MQ=255)
                      gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcag                                          >  1:1187688/1‑70 (MQ=255)
                       ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca                                        >  1:1989869/1‑71 (MQ=255)
                       ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGCAGATcatcagc                                         <  1:2521667/70‑1 (MQ=255)
                        tgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcat                                       >  1:2058547/1‑71 (MQ=255)
                         gttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcat                                       >  1:597130/1‑70 (MQ=255)
                          ttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcat                                       <  1:366564/69‑1 (MQ=255)
                          ttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcag                                          <  1:2554076/66‑1 (MQ=255)
                            ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                    <  1:353289/70‑1 (MQ=255)
                            ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                    <  1:1723589/70‑1 (MQ=255)
                            ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                    <  1:2002391/70‑1 (MQ=255)
                                tgCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCgg                                                    >  1:1872531/1‑50 (MQ=255)
                                     cATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                    <  1:1741980/61‑1 (MQ=255)
                                     cATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                    <  1:1975183/61‑1 (MQ=255)
                                      aTCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAAt                          <  1:2634964/70‑1 (MQ=255)
                                       tCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                    <  1:138790/59‑1 (MQ=255)
                                       tCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcatcag                                    <  1:1306867/59‑1 (MQ=255)
                                       tCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCC‑‑AGATcatcagcatcag                                    <  1:1644566/57‑1 (MQ=255)
                                         tttAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATCAtcagcatcagca                                  >  1:1281274/1‑59 (MQ=255)
                                           tAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca                                        >  1:512596/1‑51 (MQ=255)
                                           tAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca                                        >  1:75319/1‑51 (MQ=255)
                                           tAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaa                      <  1:2311679/69‑1 (MQ=255)
                                           tAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                       <  1:1554762/68‑1 (MQ=255)
                                               cTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcagca                                        >  1:887122/1‑47 (MQ=255)
                                                tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagca                                        >  1:1320303/1‑46 (MQ=255)
                                                tGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagca                                        >  1:2128191/1‑46 (MQ=255)
                                                      acaTCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCAt               >  1:1942490/1‑65 (MQ=255)
                                                        aTCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                       <  1:2137599/55‑1 (MQ=255)
                                                         tCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCaa                           >  1:506275/1‑50 (MQ=255)
                                                           gATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                     >  1:733793/1‑54 (MQ=255)
                                                           gATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc    <  1:928281/71‑1 (MQ=255)
                                                           gATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc    >  1:2561591/1‑71 (MQ=255)
                                                           gATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc    >  1:2521137/1‑71 (MQ=255)
                                                           gATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc    >  1:1458848/1‑71 (MQ=255)
                                                           gATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc    <  1:1156877/71‑1 (MQ=255)
                                                           gATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc    >  1:2449562/1‑71 (MQ=255)
                                                           gATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc    >  1:2198453/1‑71 (MQ=255)
                                                              ggTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctcct  >  1:110042/1‑70 (MQ=255)
                                                                   |                                                                
TCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCT  >  W3110S.gb/2295626‑2295757

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