Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 603,170 G→A 8.5% H214H (CAC→CAT)  btuB_1 ← Vitamin B12 transporter BtuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440603,1700GA 8.5% 145.9 / 3.8 47H214H (CAC→CAT) btuB_1Vitamin B12 transporter BtuB
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

CCTGAGGGCTTGCACCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGCCCCAGGCTGCGTTGGAATGTCTGGTTGCTG  >  KT2440/603037‑603305
                                                                                                                                     |                                                                                                                                       
ccTGAGGGCTTGCACCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCAGGCCTCGAAGCGCTCACTGAATCTGTGGCTCAGGCTCAGCTTCAACCCCTTGTTGCGGTACGcctcga                                                                                                                                        >  2:84659/1‑134 (MQ=255)
ccTGAGGGCTTGCACCGAGCGGCTCGTCAAATTCGCGGCGACCACGGCAGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGCCGTCGTg                                                                                                                                       >  3:106840/1‑136 (MQ=255)
    aGGGCTTGCACCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTc                                                                                                                                    >  1:199544/1‑135 (MQ=255)
       gCTTGCACCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCa                                                                                                                                <  4:27928/136‑1 (MQ=255)
            cACCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGCCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCt                                                                                                                           <  4:316167/136‑1 (MQ=255)
             aCCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCCTCGcgg                                                                                                                          >  3:23118/1‑134 (MQ=255)
             aCCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCt                                                                                                                           >  5:137436/1‑135 (MQ=255)
              ccGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTgg                                                                                                                         >  3:325404/1‑136 (MQ=255)
              ccGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTg                                                                                                                          <  8:330953/135‑1 (MQ=255)
                   cgGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGaa                                                                                                                     <  6:83650/135‑1 (MQ=255)
                   cgGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGaa                                                                                                                     <  8:143648/135‑1 (MQ=255)
                   cgGCTCGTCAAATTCGCGGCGACCACGGCTGGCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGaa                                                                                                                     <  5:348386/135‑1 (MQ=255)
                    ggCTCGTCAAATTCGCGGCGCCCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGCAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCGCCCTCGCGGGCGAAgg                                                                                                                   >  3:326891/1‑136 (MQ=255)
                    ggCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAGCGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAg                                                                                                                    >  5:45004/1‑135 (MQ=255)
                    ggCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGaa                                                                                                                     >  1:204040/1‑134 (MQ=255)
                    ggCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGaa                                                                                                                     >  7:298827/1‑134 (MQ=255)
                    ggCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGaa                                                                                                                     >  2:303772/1‑134 (MQ=255)
                    ggCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGaa                                                                                                                     >  8:203795/1‑134 (MQ=255)
                    ggCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAg                                                                                                                    >  1:106293/1‑135 (MQ=255)
                      cTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGAt                                                                                                                 >  4:205029/1‑136 (MQ=255)
                            aaaTTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGGGTGGCTCAGGCTCTGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCg                                                                                                           <  4:266119/136‑1 (MQ=255)
                              aTTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGACGCGCTCACTGAAGGGGTGGCTCAGGCTCAGGTTCAACGCCTGGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGCAGGAGGGCCCGgt                                                                                                         >  4:326358/1‑136 (MQ=255)
                                                     aGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCa                                                                                   >  8:219726/1‑135 (MQ=255)
                                                       cAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCATGGTCACCATCGCTg                                                                                                                          >  6:299892/1‑94 (MQ=255)
                                                       cAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCATGGTCACCATCGCTg                                                                                                                          <  5:299892/94‑1 (MQ=255)
                                                          gTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAg                                                                              <  5:196611/135‑1 (MQ=255)
                                                          gTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAg                                                                              <  2:291571/135‑1 (MQ=255)
                                                               aGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCCTCGCTGGCGAAGGATGGCCCGGGCGAGTCGGTTCCGGCAGTTTCATCCAGGCTgg                                                                         >  2:344574/1‑135 (MQ=255)
                                                                gCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGgc                                                                        <  1:271417/135‑1 (MQ=255)
                                                                        cGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGGTTCATCCAGGCTGGCGCCAAGGt                                                                >  5:180140/1‑135 (MQ=255)
                                                                                  aCTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTg                                                                                                                          >  6:303171/1‑67 (MQ=255)
                                                                                  aCTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTg                                                                                                                          <  5:303171/67‑1 (MQ=255)
                                                                                  aCTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGCAGCGGGt                                                      >  6:143132/1‑135 (MQ=255)
                                                                                         gtgtGGCTCAGGCGCAGGTGCACCGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGGCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGt                                              <  3:220656/136‑1 (MQ=255)
                                                                                               ctcaggctcaggTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCCCCATCGCTGGCGAAGGATGGCCCGGTCGCGTCGATAACGGCAGTTTAATCCATGCTGGCGCCAAGGTTGAAGCGGGTCGCGTCGTTgccg                                         >  2:228920/1‑135 (MQ=255)
                                                                                                   ggctcaggTGCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTgccgcc                                       <  1:303772/133‑1 (MQ=255)
                                                                                                     ctcaggTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGGTGAAGCGGGGCGCGTCGTTTCCCCCAgca                                   >  1:308410/1‑133 (MQ=255)
                                                                                                           ttCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGcccc                             <  8:298827/135‑1 (MQ=255)
                                                                                                           ttCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGcccc                             <  1:41344/135‑1 (MQ=255)
                                                                                                           ttCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGcccc                             <  2:106293/135‑1 (MQ=255)
                                                                                                           ttCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGcccc                             <  7:203795/135‑1 (MQ=255)
                                                                                                           ttCAACGCCTTGTTGCGGTAGGCGTCGGGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGcccc                             <  6:300343/135‑1 (MQ=255)
                                                                                                            tCAACGCCTTGTTGCGGTAGGCGTCATGGTCACCATCGCt                                                                                                                           >  2:178962/1‑40 (MQ=255)
                                                                                                            tCAACGCCTTGTTGCGGTAGGCGTCATGGTCACCATCGCt                                                                                                                           <  1:178962/40‑1 (MQ=255)
                                                                                                                       ttgCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGCCCCAGGCTGCGtt                   <  2:147898/133‑1 (MQ=255)
                                                                                                                        tgCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGgc                                                                        >  2:134766/1‑79 (MQ=255)
                                                                                                                        tgCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGgc                                                                        <  1:134766/79‑1 (MQ=255)
                                                                                                                                gCGTCGTGGGCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGCCCCAGGCTGCGTTGGAATGTCTgg        <  1:127344/135‑1 (MQ=255)
                                                                                                                                    cgtGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGGCGAGTCGTTGCCGCCAGAAAGCCCCAGGGTGCGGTGGTATGGCTGGTtgct   >  3:54789/1‑136 (MQ=255)
                                                                                                                                     gtGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGCCCCAGGCTGCGTTGGAATGTCTGGTtgctg  <  3:122141/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
CCTGAGGGCTTGCACCGCGCGGCTCGTCAAATTCGCGGCGACCACGGCTGTCGAGCAGGTTCAAGCCGGCCTCGAAGCGCTCACTGAAGGTGTGGCTCAGGCTCAGGTTCAACGCCTTGTTGCGGTAGGCGTCGTGGTCACCATCGCTGGCGAAGGATGGCCCGGTCGAGTCGATACCGGCAGTTTCATCCAGGCTGGCGCCAAGGTTGAAGCGGGTCGAGTCGTTGCCGCCAGAAAGCCCCAGGCTGCGTTGGAATGTCTGGTTGCTG  >  KT2440/603037‑603305

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: ‑