Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 346,909 C→A 8.5% L157L (CTG→CTT)  btuD_1 ← Vitamin B12 import ATP‑binding protein BtuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143346,9090CA 8.5% 145.4 / 3.0 47L157L (CTG→CTT) btuD_1Vitamin B12 import ATP‑binding protein BtuD
Reads supporting (aligned to +/- strand):  ref base C (16/27);  new base A (2/2);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTGGAACACCA  >  ACB143/346779‑347042
                                                                                                                                  |                                                                                                                                     
gTGGATCCAGGTCCGAGAACGGTTCGTCCATCAGCCGGAGGTCTGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTTCTGCATGCCGCCAGAAAGCTCATCCCGCTTCTTGTTGCGCCACTGGGTCAGGcc                                                                                                                                  <  8:254791/136‑1 (MQ=255)
gTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTAAGGcc                                                                                                                                  >  5:194898/1‑136 (MQ=255)
gTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTAAGGcc                                                                                                                                  >  5:188657/1‑136 (MQ=255)
 tGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGcc                                                                                                                                  >  1:260446/1‑135 (MQ=255)
             cGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGc                                                                                                                     >  7:32/1‑136 (MQ=255)
                  aCGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTtctc                                                                                                                >  7:297588/1‑136 (MQ=255)
                  aCGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCGTGTTGCGCCACTGGGTCAGGCCCACCCGCTCAAGCTtctc                                                                                                                >  5:215173/1‑136 (MQ=255)
                  aCGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCAACCGGCTTCATGTTGCGCCACAGGGTCAGGCCAACCAGCTCAAGCTtctc                                                                                                                >  3:148232/1‑136 (MQ=255)
                    ggTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCtcgt                                                                                                              <  2:368243/136‑1 (MQ=255)
                    ggTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCtcgt                                                                                                              <  7:173485/136‑1 (MQ=255)
                    ggTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCtcgt                                                                                                              <  5:234332/136‑1 (MQ=255)
                    ggTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCtcgt                                                                                                              <  4:34649/136‑1 (MQ=255)
                    ggTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCtcgt                                                                                                              <  4:365841/136‑1 (MQ=255)
                     gTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCtcgtc                                                                                                             <  4:237841/136‑1 (MQ=255)
                                agcagGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTg                                                                                                   <  6:119151/135‑1 (MQ=255)
                                agcagGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTg                                                                                                   <  5:29595/135‑1 (MQ=255)
                                          cGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCtt                                                                                        >  3:105492/1‑136 (MQ=255)
                                              atccatcGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTcc                                                                                    <  8:80496/136‑1 (MQ=255)
                                                  atcGCCAGCGCCCGTGCCAGGCCCACTCGCTGCTGGATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGTGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGc                                                                                <  8:155998/136‑1 (MQ=255)
                                                   tcGCCAGCGCCCGGGCCAGGCCATCGCGCTGCTTCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGGCGACCAGCTTGCGCAGTTCCTTTTCCGGGCg                                                                               <  6:36951/136‑1 (MQ=255)
                                                             ccGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCAtg                                                                     >  5:298320/1‑135 (MQ=255)
                                                                 gCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCa                                                                 <  8:285825/136‑1 (MQ=255)
                                                                 gCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCa                                                                 <  8:51766/136‑1 (MQ=255)
                                                                 gCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCa                                                                 <  8:62185/136‑1 (MQ=255)
                                                                         cACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGCCCCTGCATTTCCAGGCCGAAg                                                         >  1:300076/1‑136 (MQ=255)
                                                                         cACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTGGCGCCACTGGGTCAGGCCCACCAGCGCAAGCTTCTCGTCGACCAGCTTGCGCCGTGCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGCAg                                                         >  7:249125/1‑136 (MQ=255)
                                                                          aCGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTAAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGc                                                        <  6:188657/136‑1 (MQ=255)
                                                                          aCGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTAAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGc                                                        <  6:194898/136‑1 (MQ=255)
                                                                            gcgcTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGGCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTg                                                      <  8:72390/136‑1 (MQ=255)
                                                                              gctgctGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGCAGCTGAt                                                    >  3:361855/1‑136 (MQ=255)
                                                                              gctgctGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGag                                                    >  6:321226/1‑135 (MQ=255)
                                                                              gctgctGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGAt                                                    <  3:100014/136‑1 (MQ=255)
                                                                              gctgctGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGGTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGAt                                                    >  2:46072/1‑136 (MQ=255)
                                                                                     cATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTtctc                                              <  5:82498/135‑1 (MQ=255)
                                                                                             cAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACGCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGt                                     >  1:132679/1‑136 (MQ=255)
                                                                                             cAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGt                                     >  8:241433/1‑136 (MQ=255)
                                                                                                       aTCCGGCTTCTTGTTGCGCCAGTGGGTCAGGCACACCAGCTCAAGCTTATCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAgg                            <  8:249125/135‑1 (MQ=255)
                                                                                                          cGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCAGTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCAt                         >  3:308407/1‑135 (MQ=255)
                                                                                                                   gttgCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCa                  <  6:279595/133‑1 (MQ=255)
                                                                                                                       cgcCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTc            <  5:242785/135‑1 (MQ=255)
                                                                                                                       cgcAACTGGGTCAGGACCACCAGCTCAAGCTTCTCGGCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCGAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCTAGGCATCAGGGCAAACTTc            <  6:41505/135‑1 (MQ=255)
                                                                                                                        gcCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTg          <  8:230512/136‑1 (MQ=255)
                                                                                                                        gcCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTg          <  5:162301/136‑1 (MQ=255)
                                                                                                                        gcCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTg          <  4:153489/136‑1 (MQ=255)
                                                                                                                          cACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTGGa        <  5:75278/136‑1 (MQ=255)
                                                                                                                               ggTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTGGAaca     <  6:11585/134‑1 (MQ=255)
                                                                                                                                gTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTGGAACACCa  >  1:59376/1‑136 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGATGTTCTCGCGCACCGTCAGCCAAGGCATCAGGGCAAACTTCTGGAACACCA  >  ACB143/346779‑347042

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