Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 337,061 G→C 8.0% G95A (GGC→GCC)  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,0610GC 8.0% 163.2 / 6.6 50G95A (GGC→GCC) yahFputative NAD(P)‑binding succinyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base C (2/2);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCC  >  NC_000913/336913‑337209
                                                                                                                                                    |                                                                                                                                                    
gTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCgg                                                                                                                                                      <  4:535377/149‑1 (MQ=255)
  gCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCt                                                                                                                                                    <  4:8351/149‑1 (MQ=255)
     gaagaaTTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGc                                                                                                                                                 >  2:513214/1‑149 (MQ=255)
      aagaaTTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACC‑AAGCGCAAAGCGGCTCGCAc                                                                                                                                               >  4:518923/1‑149 (MQ=255)
          aTTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAg                                                                                                                                            <  1:285689/149‑1 (MQ=255)
           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                                                                                                                              <  4:489430‑M2/77‑14 (MQ=255)
                                                                                                gAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGCCTCGTGCGAGCcgctttgcgcttt                                                                                                                              >  3:489430‑M2/1‑64 (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                                                                        <  2:128355/107‑1 (MQ=255)
                                                                                                                        gaaCTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGt                                                                        >  1:128355/1‑107 (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:518923/149‑1 (MQ=255)
                                                                                                                                         agcgcaaagcgGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGa             <  3:410288/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)
                                                                                                                                                    |                                                                                                                                                    
GTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATTTCGGTCAACGGTCTGTTTGCCGCTCGCGAAGCGCGTCAGGCGCTGCAAAACGATCTCAACGTGATGCTGTTTTCC  >  NC_000913/336913‑337209

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