Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,502,376 A→C 40.5% V17G (GTC→GGC)  cysG_1 ← Siroheme synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,502,3760AC40.5% 30.2 / 4.7 37V17G (GTC→GGC) cysG_1Siroheme synthase
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base C (14/2);  total (25/12)
Fisher's exact test for biased strand distribution p-value = 3.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGCACTAGGATCGAACCGTCGGCCTGCTGGTAGTT  >  KT2440/4502245‑4502506
                                                                                                                                     |                                                                                                                                  
cccccgcACCCCCACCTCCCCGGCCCACCGGGCCCACGGGGCG‑GGCGACGGCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGGGATCTCACCGCCACCGACGACCa                                                                                                                                 <  3:151834/129‑1 (MQ=255)
   ccTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCCAGGCGCTCTCACCGCCACCGACGACCAGCa                                                                                                                              >  7:324538/1‑135 (MQ=255)
   ccTCACCCCCACATCCCCGGG‑CAACGGGGCCAAAGGGCCGGGCGACGGCCGGCGCCACCACGCGAAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCa                                                                                                                              <  4:134347/135‑1 (MQ=255)
       aCCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCCGCC‑CGCCGCCGGCATCGGCCCGCAGGCGCGCCTTGCGCAAGGCGCTCTCCCCGCCACCGCCGCCCCGCCcgcg                                                                                                                          >  1:291728/1‑135 (MQ=255)
         cGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCCAGGCGCTCTCACCGCCACCGACGCCCCGCACGCGgc                                                                                                                        >  6:80825/1‑135 (MQ=255)
         cGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGCCGCCCAGCCCGCGgc                                                                                                                        >  5:373926/1‑135 (MQ=255)
         cGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCATGCGACACAACGCCCAGCA‑CGACGCCGGCATCGGCCAGCACGCGCGCCTTGCGCAAGGCGATCTAACAGCCACCGACGCCCAGCACGCGgc                                                                                                                        >  6:309950/1‑135 (MQ=255)
         cGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACCACCCGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCACGCGCGCCGTGCGCAAGGCGAACTCACCCCCACCGACCCCcaccacccccc                                                                                                                        >  2:130296/1‑130 (MQ=255)
         cGCCACCTTCCCGGGCCAACGCGACCAACTGGCC‑GTCGACGTCAGGCGCCACCACGCGCAGCA‑CGCCGCCGGCAACGGCCAGCAGGCGCGCCTTGCGCCAGGCGATCTCACCGCCACCGACGACCAGCACGCGgc                                                                                                                        >  4:55969/1‑135 (MQ=255)
             aCCTCCCCGGCCCAACGCGGCCAACTGGCC‑GGCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGgccgcc                                                                                                                    <  7:34017/135‑1 (MQ=255)
              ccTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCCCCACGCGCCGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCCCGGCGAGCTCACCGCCCCCGCCGCCCCGCACGCGGCCGccc                                                                                                                   >  1:99012/1‑135 (MQ=255)
                 tCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGc                                                                                                                >  2:272692/1‑135 (MQ=255)
                         cAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTtgtg                                                                                                        >  8:272514/1‑135 (MQ=255)
                             gcgGCCAACTGGCC‑GTCGCCGTCCGGCGCCACCACGCGCCGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCGTGCGCAAGGCGCTCTCACCGCCCCCGCCGCCCAGCACGCGGCCGCCCGGCAGCTGGTGGa                                                                                                      >  7:62163/1‑133 (MQ=255)
                               ggCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAACGCCATCTAACCGCCCCCCACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGaaacg                                                                                                  >  6:154476/1‑132 (MQ=255)
                                    aCTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGCTCTCACCGCCCCCGACGCCCAGCCCGCGGCCGCCCGGCAGCTTGTGGCACAGCGGcc                                                                                             >  7:260733/1‑134 (MQ=255)
                                    aCTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCGTGCGCCAGGCGCTCTCCCCGCCACCCCCGCCCAGCCCCCGGCCGCCCGGCCCCGGGGGGGACAGCGGcc                                                                                             >  1:210193/1‑134 (MQ=255)
                                                 cGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGCTCTCACCGCCCCCGCCGCCCCGCCCGCGGCCCCCCGGCAGCGTGTGGAACCGCGGCCGGGAATCCAttt                                                                                 >  4:303825/1‑135 (MQ=255)
                                                 cGTCCGGCGCCACCAAGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCCCCGACGACCAGCACGCCGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCaa                                                                                   >  8:125356/1‑132 (MQ=255)
                                                   gcgggCGCCACCCGGCGC‑GCACCGCGGCGGGCAGGGGCCAGCCGGGGCGCGTGGCGCACGGGGATCTCACCGCCACCGACCACCCGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAg                                                                               <  2:210193/132‑1 (MQ=255)
                                                     gggCGCCCCCCCGCGCCCCC‑CGCCGCCGGCTGCGGCCCGCCGGCGCGCCTTGCGCACGGGGAGCACACCGCCCCCGCCGACACCCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCg                                                                             <  1:338029/134‑1 (MQ=255)
                                                                cgcgCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGCCGACCCGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTTATCCATTTAGCGGATGACCTCgc                                                                  >  1:102656/1‑134 (MQ=255)
                                                                   gcagca‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGAcgc                                                               <  5:142659/135‑1 (MQ=255)
                                                                    ctgca‑CGCCGCCGGCTTCGGCCAGCAGGGGCGCGTTGCGCAAGGCGATCTCACCGCCACCGACGACCACCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGAcgcc                                                              <  4:146549/133‑1 (MQ=255)
                                                                        a‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCCCAGCCACCGCCGACCAGCACGCGGCCGCCCTGCAGCGTGGGGTACAGCGGCACGTAATCCATTTAGCGGATGCCCTCgccgccgcc                                                           >  4:9198/1‑134 (MQ=255)
                                                                          cgccgcCGGCATCGGCAAGCAGGCGCGCCTTGCGCAAGGCGAGCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCa                                                         <  7:46670/135‑1 (MQ=255)
                                                                           gccgccGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGCTCCCCCCGCCCCCGCCGCCCAGCCCGCGGCCGCCCGGCCGCGTGTGGCACAGCGGCAGGTAATCCATTTAGCGGATGGCCTCGCCGCCGCCCag                                                        >  3:318248/1‑134 (MQ=255)
                                                                                cgGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCCCCGACCCCCCGCACGCGGCCGCCCTGCAGCTTGTGGCACCGCGGCAGGTAATCCATTTAGCGGATGACCTCGCCGCCGCCCCTGTAAg                                                   >  1:21763/1‑135 (MQ=255)
                                                                                cgGCATCGGCCAGCAGGCGCGCCGTGCGCAAGGCGATCTCACCGCCCCCGACGCCCAGCCCGCGGCCGCCCTGCAGCTTGTGGCACAGCGGCAGGTAATCCATTTAGCGGATGCCCTCGCCGCCGCCCATGTAAg                                                   >  2:261988/1‑135 (MQ=255)
                                                                                          cagcagGCGCGCCTTGCGCAAGGCGATCTCCCCGCCCCCGACGCCCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCGTCAGCAc                                         >  5:218671/1‑135 (MQ=255)
                                                                                                gcgcgcCTTGCGCAAGGCGATCTCACCGCCCCCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCgg                                   >  5:267236/1‑135 (MQ=255)
                                                                                                gcgcgcCTTGCGCAAGGCGATCTCACCGCCACCGACCACCAGCCCGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCgg                                   <  2:99012/135‑1 (MQ=255)
                                                                                                 cgcgcCTTGCGCAAGGCGATCTCACCGCCACCGCCGCCCCGCCCGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGc                                  <  5:80825/135‑1 (MQ=255)
                                                                                                                     ctcACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGGGGAACAGCGGAAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGCACTAGGATCGAACCGTc                 <  6:161339/132‑1 (MQ=255)
                                                                                                                       cACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGCACTAGGATCGAACCGTCGGCct            >  5:363957/1‑135 (MQ=255)
                                                                                                                       cACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGCACTAGGATCGAACCGTCGGCct            >  1:135308/1‑135 (MQ=255)
                                                                                                                          cGCCACCGACGCCCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGCACTAGGATCGAACCGTCGGCctgct         <  2:249115/135‑1 (MQ=255)
                                                                                                                                 gacgacCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGCACTAGGATCGAACCGTCGGCCTGCTGGTAGtt  >  4:228123/1‑135 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCC‑GTCGACGTCCGGCGCCACCACGCGCAGCA‑CGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCACCGACGACCAGCACGCGGCCGCCCTGCAGCTTGTGGAACAGCGGCAGGTAATCCATTTAGCGGATGACCTCGACGCCGCCCATGTAAGGCTTCAGCACTTCCGGCACTAGGATCGAACCGTCGGCCTGCTGGTAGTT  >  KT2440/4502245‑4502506

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