Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 601,357 C→T 10.5% P60S (CCT→TCT)  btuF → Vitamin B12‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440601,3570CT10.5% 111.9 / 4.5 38P60S (CCT→TCT) btuFVitamin B12‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (18/16);  new base T (2/2);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGGCTGCAGAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGCCCAGCGCGACCAACTCAAGCGCCTGGGTATCGCC  >  KT2440/601226‑601485
                                                                                                                                   |                                                                                                                                
tGGCTGCAGAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCtt                                                                                                                               >  8:101539/1‑135 (MQ=255)
tGGCTGCAGAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCtt                                                                                                                               >  4:181394/1‑135 (MQ=255)
tGGCTGCAGAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCtt                                                                                                                               >  7:148751/1‑135 (MQ=255)
tGGCTGCAGAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCtt                                                                                                                               >  8:65097/1‑135 (MQ=255)
        gAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCggg                                                                                                                       >  2:34583/1‑135 (MQ=255)
        gAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCggg                                                                                                                       >  4:117740/1‑135 (MQ=255)
        gAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCggg                                                                                                                       >  1:210073/1‑135 (MQ=255)
          aCCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGgcg                                                                                                                     >  6:171603/1‑135 (MQ=255)
                 cGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGCCCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAgg                                                                                                               >  8:89978/1‑134 (MQ=255)
                  gtgtCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGc                                                                                                             <  2:166841/135‑1 (MQ=255)
                  gtgtCGGCAGGCGGGCCCCATCTATGAGCGAGATCATGCTTGAGCTGAAGGCCGACGACCTGCTGGTGGGCGGGCTCTATGGCGGAGAACGGCCGGCTGCACTGCGCCACTTGTCTTCGGTCGGGCGCCAAGGGc                                                                                                             <  8:246939/135‑1 (MQ=255)
                           gTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGTCTTCGGTCGGGCGCCAAGGGc                                                                                                             <  8:242896/126‑1 (MQ=255)
                           gTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGTCTTCGGTCGGGCGCCAAGGGc                                                                                                             >  7:242896/1‑126 (MQ=255)
                           gTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGTCTTCGGTCGGGCGCCAAGGGCAACTGGACa                                                                                                    >  7:246939/1‑135 (MQ=255)
                                 cccccTCGATGAGGGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGa                                                                                                <  4:263947/133‑1 (MQ=255)
                                 cccccTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGa                                                                                                <  1:34583/133‑1 (MQ=255)
                                    ccTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGcc                                                                                           >  1:213081/1‑135 (MQ=255)
                                          tGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCa                                                                                     <  2:187247/135‑1 (MQ=255)
                                                    cATGCTGGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGccc                                                                           <  4:155947/135‑1 (MQ=255)
                                                      tGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGcccc                                                                          <  7:101539/134‑1 (MQ=255)
                                                      tGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCg                                                                         <  6:288010/135‑1 (MQ=255)
                                                               tGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTtgc                                                                >  2:134315/1‑135 (MQ=255)
                                                                 gcggCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTtgctg                                                              <  6:57167/133‑1 (MQ=255)
                                                                    gCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGctgctg                                                           >  1:273823/1‑135 (MQ=255)
                                                                         cgacCTGCTGGTGGGCGTGCGCGAGGGCGGCGAACGGCCGGCGGCACTGCGCCACTTGCCTTCGCTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCGGCGCCCCGACCTGTTGCTGCTGTGGcc                                                      <  4:9008/135‑1 (MQ=255)
                                                                            ccTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAg                                                    >  8:152623/1‑134 (MQ=255)
                                                                                 tggtggGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAgtg                                              >  7:269105/1‑135 (MQ=255)
                                                                                           gCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGc                                    >  6:83627/1‑135 (MQ=255)
                                                                                             tCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCACCTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCGGTGGCCAGGCAGTGTGCCGCCGGc                                    >  7:259090/1‑133 (MQ=255)
                                                                                             gcgctGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGGCGGGCGCCAATGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCGGTGGCCAGGCAGTGTGCCGCCGGccc                                  <  3:211833/131‑1 (MQ=255)
                                                                                             gcGATGGCGGCGAACGGCCGGCGGCACTGCGCCACTTGCCTGCGGGCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGccc                                  <  7:291551/134‑1 (MQ=255)
                                                                                                  ggcggcGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGCCCAgcgc                             <  2:210073/135‑1 (MQ=255)
                                                                                                  ggcggcGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGGCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGCCCAgcgc                             <  3:117740/135‑1 (MQ=255)
                                                                                                        gAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGCCCAGCGCGACCaa                       <  2:176169/135‑1 (MQ=255)
                                                                                                         aaCGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGCCCAGCGCGACCAAc                      >  1:290398/1‑135 (MQ=255)
                                                                                                                     cACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGACAGTGTGCCGCCGGCCCAGCGCGACCAACTCAAGCGCCTgg          <  2:213081/135‑1 (MQ=255)
                                                                                                                          cgcCACTGGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTGGCTGCTGTGGCCAGGCGGTGTGCCGCCGGCCCAGCGCGACCAACTCAAGCGCCTGGGTATc     <  7:27710/135‑1 (MQ=255)
                                                                                                                             aaCTTGCCTTCGGTAGGGCGCCAAGGGCAACAGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCAGCCGGCCCAGAGCGACCAACTCAAGCGCCTGGGTATCGcc  >  4:246834/2‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                
TGGCTGCAGAACCGCTACGTGTCGTCAGTCTGGCCCCCTCGATGAGCGAGATCATGCTTGAGCTGCAGGCCGACGACCTGCTGGTGGGCGTGCTCGATGGCGGCGAACGGCCGGCTGCACTGCGCCACTTGCCTTCGGTCGGGCGCCAAGGGCAACTGGACATGGAGCGCCTGCTCAGCCTGCGCCCCGACCTGTTGCTGCTGTGGCCAGGCAGTGTGCCGCCGGCCCAGCGCGACCAACTCAAGCGCCTGGGTATCGCC  >  KT2440/601226‑601485

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