Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 337,072 G→C 8.2% A99P (GCG→CCG)  yahF → putative NAD(P)‑binding succinyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913337,0720GC 8.2% 160.1 / 6.5 49A99P (GCG→CCG) yahFputative NAD(P)‑binding succinyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base C (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

TTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTC  >  NC_000913/336924‑337218
                                                                                                                                                    |                                                                                                                                                  
tttAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAgg                                                                                                                                                    >  4:63674/1‑149 (MQ=255)
     gACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGt                                                                                                                                               >  1:428918/1‑149 (MQ=255)
          ctgACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGc                                                                                                                                          <  3:336723/149‑1 (MQ=255)
           tgACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGcc                                                                                                                                         >  1:100710/1‑149 (MQ=255)
            gACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCa                                                                                                                                        <  1:190571/149‑1 (MQ=255)
                 cGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTAtt                                                                                                                                   >  1:230749/1‑149 (MQ=255)
                         gAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAgcgc                                                                                                                           >  4:107647/1‑149 (MQ=255)
                               gCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCaaaaa                                                                                                                     <  4:318787/149‑1 (MQ=255)
                                              cTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGCCTCGTGCGAGCcgcttt                                                                                                                                              >  1:18318‑M2/1‑103 (MQ=255)
                                              cTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGCCTCGTGCGAGCcgcttt                                                                                                                                              <  2:18318‑M2/109‑7 (MQ=255)
                                                   gatTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAAcc                                                                                                 <  2:205449/149‑1 (MQ=255)
                                                      tGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTgg                                                                                              >  1:359671/1‑149 (MQ=255)
                                                           tcaATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTg                                                                                         >  4:547039/1‑149 (MQ=255)
                                                              aTGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGAtt                                                                                      <  4:547446/149‑1 (MQ=255)
                                                                 gTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCg                                                                                   <  3:107647/149‑1 (MQ=255)
                                                                  tAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCgg                                                                                  >  3:542302/1‑149 (MQ=255)
                                                                     aTCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCa                                                                               >  4:102982/1‑149 (MQ=255)
                                                                        gggTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACg                                                                            <  4:336213/149‑1 (MQ=255)
                                                                        gggTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACg                                                                            >  4:177402/1‑149 (MQ=255)
                                                                              ggACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGt                                                                      >  2:252612/1‑149 (MQ=255)
                                                                                 cAACGAGCAGTCACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTg                                                                   >  3:342571/1‑149 (MQ=255)
                                                                                  aaCGAGCAGTTACTGGTGGAGATCGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGc                                                                  >  4:100071/1‑149 (MQ=255)
                                                                                     gAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGc                                                               >  2:248687/1‑149 (MQ=255)
                                                                                     gAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGCCTCGTGCGAGCcgctttgcgcttt                                                                                                                                       >  3:489430‑M2/1‑64 (MQ=255)
                                                                                     gAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGCCTCGTGCGAGCcgctttgcgcttt                                                                                                                                       <  4:489430‑M2/77‑14 (MQ=255)
                                                                                          gTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTcgcg                                                          >  3:96561/1‑149 (MQ=255)
                                                                                          gTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGCTGATTTAGGTCAACGGTCTGTTTGCCGCTcgcg                                                          >  4:80978/1‑149 (MQ=255)
                                                                                           ttACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGttt                                                                    <  4:334530/138‑1 (MQ=255)
                                                                                           ttACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCAGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGttt                                                                    >  3:334530/1‑138 (MQ=255)
                                                                                                  tggAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGt                                                  <  1:248687/149‑1 (MQ=255)
                                                                                                   ggAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTc                                                 <  3:547039/149‑1 (MQ=255)
                                                                                                          gaagaaCTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCt                                          <  3:279754/149‑1 (MQ=255)
                                                                                                           aagaaCTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAgc                                                                                                                             <  1:563938/65‑1 (MQ=255)
                                                                                                           aagaaCTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAgc                                                                                                                             >  2:563938/1‑65 (MQ=255)
                                                                                                             gaaCTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGt                                                                                 >  1:128355/1‑107 (MQ=255)
                                                                                                             gaaCTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGt                                                                                 <  2:128355/107‑1 (MQ=255)
                                                                                                                   ttCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGa                                 <  4:155216/149‑1 (MQ=255)
                                                                                                                      aaCACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGAtct                              >  3:499532/1‑149 (MQ=255)
                                                                                                                       acacCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGAtctc                             >  2:316718/1‑149 (MQ=255)
                                                                                                                              agcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGa                      <  3:410288/149‑1 (MQ=255)
                                                                                                                              agcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGa                      <  3:518923/149‑1 (MQ=255)
                                                                                                                               gcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGAt                     <  3:411716/149‑1 (MQ=255)
                                                                                                                                 gcaaagctGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGc                   >  2:324214/1‑149 (MQ=255)
                                                                                                                                  caaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCt                  <  2:428918/149‑1 (MQ=255)
                                                                                                                                      gcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGa                                                                                        <  2:83261/75‑1 (MQ=255)
                                                                                                                                      gcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGa                                                                                        >  1:83261/1‑75 (MQ=255)
                                                                                                                                         gCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTcc           <  3:252784/149‑1 (MQ=255)
                                                                                                                                               aCGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAAc     <  3:61216/149‑1 (MQ=255)
                                                                                                                                                 gAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGt   <  3:321656/149‑1 (MQ=255)
                                                                                                                                                  aGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGtc  <  3:100071/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCCGATAACGTC  >  NC_000913/336924‑337218

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