Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,164,740 G→A V368I (GTC→ATC)  btuB → cobalamin/cobinamide outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,164,7400GA97.1% 115.6 / ‑3.9 34V368I (GTC→ATC) btuBcobalamin/cobinamide outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (17/16);  total (18/16)
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.84e-01

ACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAAC  >  NC_000913/4164662‑4164816
                                                                              |                                                                            
acgacgCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATc                                                                            >  2:827643/1‑81 (MQ=255)
   acgCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGc                                                                         <  1:827610/81‑1 (MQ=255)
              aGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCtt                                                              >  1:827655/1‑81 (MQ=255)
                  tATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTg                                                            >  2:827648/1‑79 (MQ=255)
                   aTGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGaa                                                          <  1:827611/80‑1 (MQ=255)
                    tGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAgg                                                        <  2:827612/81‑1 (MQ=255)
                        gAGGATGGATATGATCAACGTAATACAGGCGTCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGgcgc                                                     <  2:827613/80‑1 (MQ=255)
                         aGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCgcag                                                   <  2:827615/81‑1 (MQ=255)
                         aGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCgca                                                    <  1:827614/80‑1 (MQ=255)
                                   tGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTg                                          <  1:827617/80‑1 (MQ=255)
                                   tGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTg                                          <  2:827616/80‑1 (MQ=255)
                                    gATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGAc                                        >  1:827631/1‑81 (MQ=255)
                                      tCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGa                                      >  1:827653/1‑81 (MQ=255)
                                      tCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGa                                      >  2:827665/1‑81 (MQ=255)
                                       cAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGa                                      >  2:827662/1‑80 (MQ=255)
                                        aaCGTAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGAt                                     >  2:827673/1‑80 (MQ=255)
                                            tAATACCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTc                                <  1:827618/81‑1 (MQ=255)
                                                aCCGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAg                            >  1:827649/1‑81 (MQ=255)
                                                 ccGGCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATGACTCACAGt                           <  2:827619/81‑1 (MQ=255)
                                                   ggCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGttt                         >  2:827694/1‑81 (MQ=255)
                                                   ggCATCTATCTGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGttt                         <  2:827620/81‑1 (MQ=255)
                                                      atctatctGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGgt                      >  2:827658/1‑81 (MQ=255)
                                                          atctGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGgtcgt                   >  1:827672/1‑80 (MQ=255)
                                                           tctGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGgtc                     <  2:827621/77‑1 (MQ=255)
                                                            ctGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCAt                >  2:827654/1‑81 (MQ=255)
                                                            ctGACCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCAt                >  2:827656/1‑81 (MQ=255)
                                                               aCCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCAt                <  2:827623/78‑1 (MQ=255)
                                                               aCCGGGCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGa             <  2:827622/81‑1 (MQ=255)
                                                                  gggCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAAcc          >  1:827660/1‑81 (MQ=255)
                                                                   ggCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAACCt         <  2:827624/81‑1 (MQ=255)
                                                                   ggCTGCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAACCt         <  2:827625/81‑1 (MQ=255)
                                                                       gCAACAAATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAACCTGGc      >  2:827681/1‑80 (MQ=255)
                                                                         aacaaGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAACCTGGCaaa   >  2:827663/1‑81 (MQ=255)
                                                                          acaaATCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAAc  >  2:827666/1‑81 (MQ=255)
                                                                              |                                                                            
ACGACGCCGGGTACAGGTTATGTTGAGGATGGATATGATCAACGTAATACCGGCATCTATCTGACCGGGCTGCAACAAGTCGGCGATTTTACCTTTGAAGGCGCAGCACGCAGTGACGATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAAC  >  NC_000913/4164662‑4164816

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