Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,488,952 T→C S73G (AGT→GGT)  yfdE ← acetyl‑CoA:oxalate CoA‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,488,9520TC98.0% 169.7 / NA 50S73G (AGT→GGT) yfdEacetyl‑CoA:oxalate CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (30/19);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGCGGGTTGATTTCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACTCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGT  >  NC_000913/2488804‑2489097
                                                                                                                                                    |                                                                                                                                                 
cgcgGGTTGATTTCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACAccct                                                                                                                                                 >  1:530901/1‑151 (MQ=255)
   gggTTGATTTCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCtt                                                                                                                                              <  2:675542/151‑1 (MQ=255)
     gTTGATTTCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTg                                                                                                                                            <  1:169710/151‑1 (MQ=255)
       tGATTTCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGcc                                                                                                                                          >  2:524065/1‑151 (MQ=255)
       tGATTTCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGc                                                                                                                                           <  2:754254/150‑1 (MQ=255)
            tCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGa                                                                                                                                      <  1:114365/150‑1 (MQ=255)
                     gTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATaaaa                                                                                                                             >  1:108745/1‑150 (MQ=255)
                          ccATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTa                                                                                                                        >  1:66926/1‑150 (MQ=255)
                            aTGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAAt                                                                                                                      >  1:217968/1‑150 (MQ=255)
                                       aGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGtc                                                                                                          >  1:215299/1‑151 (MQ=255)
                                         tttttCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTATGAGCATATTTATAAATATACTTTTAGCGTGATTATTCTTTAAATCAAGAACCACTCCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACt                                                                                                             <  1:147645/146‑1 (MQ=255)
                                          ttttCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGtcca                                                                                                        <  2:759273/150‑1 (MQ=255)
                                                  tGTTCCTGGGCGAAAATTCTCAGCTAATACAGCAGCGTGTGTGAGCATATTTATAAATATACTTTTATCGTGATAATTCTTTAAATCAAGAACAACTCTCTCTTTGCCATTATTAATAAAACTGTAAGAGAGTGACTGTCCAACcacaaa                                                                                                <  1:669991/150‑3 (MQ=255)
                                                  tGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACAt                                                                                                 >  1:509396/1‑149 (MQ=255)
                                                   gTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAgg                                                                                              >  1:197899/1‑151 (MQ=255)
                                                   gTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAg                                                                                               <  2:218000/150‑1 (MQ=255)
                                                    tACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAggg                                                                                             >  1:83441/1‑151 (MQ=255)
                                                     aCCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAggg                                                                                             <  2:22807/150‑1 (MQ=255)
                                                       cTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGAcc                                                                                          >  2:322843/1‑151 (MQ=255)
                                                            cGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAAtg                                                                                     >  2:324028/1‑151 (MQ=255)
                                                            cGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAAtg                                                                                     <  1:235773/151‑1 (MQ=255)
                                                             gAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAAtgt                                                                                    <  2:197932/151‑1 (MQ=255)
                                                               aaaTTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGc                                                                                  >  1:238838/1‑151 (MQ=255)
                                                                aaTTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCg                                                                                 >  2:586553/1‑151 (MQ=255)
                                                                 aTTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCgg                                                                                >  2:411892/1‑151 (MQ=255)
                                                                 aTTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCgg                                                                                >  2:402388/1‑151 (MQ=255)
                                                                          tAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTatcatc                                                                        >  2:519238/1‑150 (MQ=255)
                                                                           aaTACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAACCTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTAtcatca                                                                       <  1:683868/150‑1 (MQ=255)
                                                                               cATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATg                                                                  <  1:563795/151‑1 (MQ=255)
                                                                                aTCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATg                                                                  >  2:97319/1‑150 (MQ=255)
                                                                                       tGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAACCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCcggc                                                           <  2:424149/150‑1 (MQ=255)
                                                                                          ttGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCt                                                       >  2:103130/1‑151 (MQ=255)
                                                                                                 tatTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACttt                                                >  1:462573/1‑151 (MQ=255)
                                                                                                    ttATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAAt                                             >  1:325538/1‑151 (MQ=255)
                                                                                                    ttATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAAt                                             >  2:716293/1‑151 (MQ=255)
                                                                                                     tataAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAAtt                                            >  2:184420/1‑151 (MQ=255)
                                                                                                      ataAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAAtt                                            >  2:79669/1‑150 (MQ=255)
                                                                                                            ataCTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCtt                                      <  2:181511/150‑1 (MQ=255)
                                                                                                            ataCTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCtt                                      <  1:478834/150‑1 (MQ=255)
                                                                                                            ataCTTTTATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTg                                     <  2:114679/151‑1 (MQ=255)
                                                                                                                   tATCGTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCat                              >  2:12980/1‑151 (MQ=255)
                                                                                                                       gTGATCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTa                          <  2:2499/151‑1 (MQ=255)
                                                                                                                          aTCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGAACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACa                        >  1:552488/1‑150 (MQ=255)
                                                                                                                           tCATTCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACaa                       >  2:31583/1‑150 (MQ=255)
                                                                                                                              ttCTTTAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAaaga                    >  2:399638/1‑150 (MQ=255)
                                                                                                                                  ttAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGt                 <  2:398966/149‑1 (MQ=255)
                                                                                                                                  ttAAATCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTg               >  2:321522/1‑151 (MQ=255)
                                                                                                                                    aaaTCAAGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTg               >  2:596724/1‑149 (MQ=255)
                                                                                                                                         aaGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCg        >  2:246549/1‑151 (MQ=255)
                                                                                                                                          aGAACCACACCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCg        >  1:659927/1‑150 (MQ=255)
                                                                                                                                               cacaCCCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGt  <  2:160993/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCGGGTTGATTTCTTGAAGCGTTTCCCATGAAAACCCCAGTTTTTCCATTGTACCTGGGCGAAAATTCTCAGCTAATACATCAGCTTGTTTGAGCATATTTATAAATATACTTTTATCGTGATCATTCTTTAAATCAAGAACCACACTCTCTTTGCCATGATTAATAAAACTGTAATAGAGTGACTGTCCATCCACATAGGGACCAAATGTGCGGGTATCATCACCATGACCCGGCGGCTCAACTTTAATTACCCTTGCGCCCATATTACAAAGAAGTTGAGTTCCGAAAGGT  >  NC_000913/2488804‑2489097

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