Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,795,229 A→T 18.6% intergenic (‑77/+24) btuC ← / ← ihfA vitamin B12 ABC transporter membrane subunit/integration host factor subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,795,2290AT18.6% 51.6 / 6.6 28intergenic (‑77/+24)btuC/ihfAvitamin B12 ABC transporter membrane subunit/integration host factor subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base T (3/2);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTA  >  NC_000913/1795090‑1795364
                                                                                                                                           |                                                                                                                                       
accTGCCAGGCATAATAACCAGCGAATCTTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCg                                                                                                                                      <  1:113419/141‑1 (MQ=255)
   tGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGcc                                                                                                                                   >  1:154650/1‑143 (MQ=255)
     aCAGGCATAATAACCAGCGAATATTGTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCtt                                                                                                                                 <  1:457377/143‑1 (MQ=255)
      cAGGCATAACAACCGGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCAGGGTATCCGTTCGGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCttt                                                                                                                                <  1:165876/143‑1 (MQ=255)
       gggCATAGGAACCAGTGAAGTGTTGGTCGCTGTTGTTGGCGGGCACGTGTCAGCATGTTATCCGTTCTGCTGAAGTGTCAGGGCGTTGATTTTCCGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCtttt                                                                                                                               <  2:466194/142‑1 (MQ=255)
             aagctgcAGCGGATAGTTTGTCGCGGTTGTTGGCGGGCAAGTGTCAGCAGGGTATCCGTTCTGCTGAAGTGTCAGGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtta                                                                                                                         <  1:96477/137‑1 (MQ=255)
                       cgatatTTTCTTGTTGTTGTTGGCAGTGAAGTGTCAGTATGTTATTCGTTCTGCTGTAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGATAGTTAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTa                                                                                                               <  2:279383/141‑1 (MQ=255)
                          tatTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTc                                                                                                            >  2:73478/1‑143 (MQ=255)
                          gatTTGGGCGCTGGTGGTGGGGGGCAAGTGTCAGCATGGGAGCGGTGCTGCTGAAGTGTGATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGTGCGGCCTTTTTAGTTAGATCAGATTACTc                                                                                                            <  2:132087/142‑1 (MQ=255)
                                  cgggggtgttgGCGGGCAAGTGTCAGCATGGCATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTgg                                                                                                    <  2:461637/137‑1 (MQ=255)
                                              cGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGtttt                                                                                        <  1:279798/143‑1 (MQ=255)
                                               gggCACGTGTCAGCCTGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAACAGAACAAAGGCCGCAGAGCGGCCTTTTTCGTTCGCTCCGATTACTCCTCTTTGGGCGAAGCGttata                                                                                       >  2:253782/1‑140 (MQ=255)
                                               gggCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGGAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGTGCGAAGCGTTTTc                                                                                       >  1:439067/1‑143 (MQ=255)
                                                   aaGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGAcc                                                                                   >  1:483717/1‑143 (MQ=255)
                                                        tCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGACACGTTTTCTACCCGGCt                                                                              >  2:500442/1‑143 (MQ=255)
                                                                   tCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCt                                                                   >  2:310581/1‑143 (MQ=255)
                                                                     cGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAAcctattt                                                                 >  1:367293/1‑137 (MQ=255)
                                                                                 ggggcATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg                                                     <  2:374833/139‑1 (MQ=255)
                                                                                    tCATGGCGTTGATTTTACGGTGACTCTTCGACAGTAAAAAGAAAAAAGGACGCAGCGCGGCCTTTTACGTTAGATCAGATTACTAGTCATTGGGCGAAGCATTTTCGAACCGGATCTTTAAGTTCTGACCGGGCCTGAATGTc                                                  >  2:473214/1‑143 (MQ=255)
                                                                                       gggCGTTGATTTGACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTcacc                                               <  1:419564/142‑1 (MQ=255)
                                                                                                 tttACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAAAGCGGACTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGACGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGACGGTCACCACCCGACGTg                                     >  2:165832/1‑143 (MQ=255)
                                                                                                            tctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCAGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGa                          >  2:82867/1‑143 (MQ=255)
                                                                                                                             taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCTTGCGGTAAGAGGAATAGCCTGGCCCGattg         >  2:66315/2‑142 (MQ=255)
                                                                                                                             taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGGCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCGCGCCCgcact         >  2:456494/2‑139 (MQ=255)
                                                                                                                              aaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGTGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTc        >  1:420729/1‑143 (MQ=255)
                                                                                                                                aaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCGTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGGTTTCgg      >  1:491182/1‑143 (MQ=255)
                                                                                                                                   ggCCGCAGAGCGTCCTTTTTAGTTAGATCAGATTACTCTTCTGTGGGCGAAGCGTGTTCGACCCGGCTTTTTGCCTTCTGCCCGTGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGtt   <  2:413338/143‑1 (MQ=255)
                                                                                                                                    cccGCAGGGCGGCCTTTCGAGTTAGATCCGATTACTCGTCTCTGGGCGAAGCGTTTTCGACCCGGCTTGTTAACTACTGCCCGGGGCTGAAGGCCACCACGCGCCGGGCTGTAATGGGAATATTCTCGCCCGTTTTCGGGTTa  <  2:232690/142‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                       
AACTGACAGGCATAATAACCAGCGAATATTTTGTCGCTGTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCGTGCTGTAATGGGAATATCCTCGCCCGTTTTCGGGTTA  >  NC_000913/1795090‑1795364

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: