Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,489,152 C→T 59.7% S6N (AGC→AAC)  yfdE ← acetyl‑CoA:oxalate CoA‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,489,1520CT59.7% 14.5 / 24.2 37S6N (AGC→AAC) yfdEacetyl‑CoA:oxalate CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base C (9/5);  new base T (8/13);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCACAAGGGA  >  NC_000913/2489018‑2489278
                                                                                                                                      |                                                                                                                              
gggTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCCAAAGGTCCATTAAAGACATTTGTCATATCGAAAATAATAAACATTCAAACGGCCCTTTGCTTTCattat                                                                                                                       >  1:446193/1‑143 (MQ=255)
        tcaCCATGACCCGGCGACTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGATCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAAATAATAAGCCTTCAAACGACCCTTTGATTTCATTATTTGTCAt                                                                                                                >  2:458438/1‑143 (MQ=255)
               gACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCACACCGCCCTTTGCTTTCATTATTTTTCATTTGATTa                                                                                                         >  1:343993/1‑143 (MQ=255)
                  ccGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTAtttcttattttcttatac                                                                                                      >  1:386115/1‑141 (MQ=255)
                  ccGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTTTCATTTGTTTATGc                                                                                                      >  1:128528/1‑143 (MQ=255)
                           cAACTTTAATTACCCTTGCTCCGATATTACAAAGAAGTTGAGTTCCGAAATGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGaaaaa                                                                                             <  2:283708/143‑1 (MQ=255)
                             tcTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATATCTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAg                                                                                           <  2:252515/142‑1 (MQ=255)
                             gcTTTATTTACGCTTGCGCCTATATTACAACGTTTTTGGTTTCTGTAATGTCCATTAAGGACATGTGTCATTTCGATAACTTATAAGCTTTCAAACTGCCCTTTTTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAg                                                                                           <  2:560778/142‑1 (MQ=255)
                               tttAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTACGGACCTGTGTCATATCGATCACTACTAAGCCTTCAAACGGCACTTTATTTTAATAATTTGTCATTAGATTATGCCTTGAAAAAAGcc                                                                                         >  1:458609/1‑143 (MQ=255)
                                        ccTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTAATTATTTGTCATTTGCTTATGACTTGAAAAAAGCCGCAGGTGaa                                                                                >  2:281221/1‑143 (MQ=255)
                                          ttGCGCTCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAAGAAGCCTTCAAACGGCGCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGaaaa                                                                              <  2:166721/143‑1 (MQ=255)
                                            catgctatatTCCAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCGTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGa                                                                            <  1:185259/137‑1 (MQ=255)
                                                  atatTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTc                                                                      <  1:246527/143‑1 (MQ=255)
                                                    atTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACTGCCCTTTGCTTTCATTATTTGTCATTTGCTTATACCTTGAAAAAAGCCGAAGATGAAAAAACACTTCga                                                                    >  1:173015/1‑142 (MQ=255)
                                                         aaaGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAACAAGCCGAAGGTGAAAAGACACTTCTGCAtat                                                               >  1:238779/1‑141 (MQ=255)
                                                         aaaGAAGTTAAGTTCCGAAAAGTGCATTAAGGGTAAGAGGTATATTTATAACAAATAAGCCTTCAAACGGCCCTTTATATTCAATATTTCTCATTTGCTCATAGCGTGATAAAAAGCGAAGGTGAAAAACAACTCCCGCAttt                                                               >  2:341890/1‑143 (MQ=255)
                                                                        cGAAAGGTCCATTAAGGTCATTTGTCATTTCGATTACTAATAAGCCTTCAAACGGCCCTTGGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                                <  1:176514/143‑1 (MQ=255)
                                                                        cGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCAGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATaa                                                >  1:213890/1‑143 (MQ=255)
                                                                           aaGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACAt                                             <  1:498222/143‑1 (MQ=255)
                                                                            atgTCCATTATGGAGATGTGTCATATCGATAACTAATAAGCTTTCAAACGGCCCTTTGGTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATg                                            <  2:146624/141‑1 (MQ=255)
                                                                            acgTCCATTAAGGACATGTGTCATATCGATAACTCATAAGCCTTCAAACGGCACTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTTAAAAGACACTTCGGCATTTATGTAAAATCAATAACATg                                            <  2:250127/141‑1 (MQ=255)
                                                                                ccATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCCTGAAAAAAGCGGAAGGTGACACGACACTTCGGCATTTCCGTAAACTCAAGAACAGGaata                                        >  2:17546/1‑143 (MQ=255)
                                                                                 cATTAAGTACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCTCTTTTTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaa                                       <  1:353432/143‑1 (MQ=255)
                                                                                   ttAAGGACATGTGTCATTTCAATAACTAATAAGTCTTTAAACGGGACTTTGTTTTCATTATCTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGattaatt                                     <  1:251123/143‑1 (MQ=255)
                                                                                   ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACTGCCCTTAGTAAAAATTATTTGTCATTTGTAAATTCCTAGAAAAAAGCCGAAGGTCAAAAGACACTTCTGCATTTATGTAAAATCAATAACTTGattaatt                                     >  2:64495/1‑143 (MQ=255)
                                                                                   ttAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAGATGTTATTACGACAACTATGTAAAATCAATAACATGattaatt                                     >  1:545790/1‑143 (MQ=255)
                                                                                     aaGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTAt                                   <  1:312058/143‑1 (MQ=255)
                                                                                         aCATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGaaa                               <  2:449123/143‑1 (MQ=255)
                                                                                                ttatatCGATAGCTAATACGCCTCTATACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGACTTGAAACAAGTCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCg                        <  1:524007/141‑1 (MQ=255)
                                                                                                acATATCGATAACTAATAGGCCTTCAAACGGCCCTTTGTTTTCATGATTTGTCATTTGCTTATGCCTTGAAAAAAGTCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCg                        <  2:27946/142‑1 (MQ=255)
                                                                                                     tCGATAACTAATAAGCCTTCTCAGGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGa                   <  2:491322/143‑1 (MQ=255)
                                                                                                     tCGATAACTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGa                   <  2:12875/143‑1 (MQ=255)
                                                                                                            cTAATAAGCCTTCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGTCTTGAAAAAAGCCGAAGGTGAAAAGAGACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATaagt            <  1:198647/143‑3 (MQ=255)
                                                                                                                 aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTTTCATTTGCATATGCCTTGAACAAAGCCGAAGGTGAAAAGGCACTTCGGCATGTATGTAAAATCAATAGCATGATTAATTAAGAAAACAGACAACTGACATAAATCcaca       >  2:131473/1‑143 (MQ=255)
                                                                                                                 aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTAGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAGTAACATGATTAATTATAAAACCAGCCGACTGACATAAATCcaca       >  2:131748/1‑143 (MQ=255)
                                                                                                                 aaGCCTTCAAACGGCCCTTTGCTTTCATTATTTGACATTTGCTTATGCCTTGAAAAAAGCCGAGGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGACAACTGACATAAATCcaca       >  2:131433/1‑143 (MQ=255)
                                                                                                                      ttCAAACGGCCCTTTGTTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAACTCAATAACATGATTAATTATGAAACAAGTCGACTGACATAAATCCACAAggaa  >  1:24363/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                              
GGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGTCCATTAAGGACATGTGTCATATCGATAACTAATAAGCCTTCAAACGGCCCTTTGCTTTCATTATTTGTCATTTGCTTATGCCTTGAAAAAAGCCGAAGGTGAAAAGACACTTCGGCATTTATGTAAAATCAATAACATGATTAATTATGAAACCAGTCGACTGACATAAATCCACAAGGGA  >  NC_000913/2489018‑2489278

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