Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 578,044 C→A 8.7% Q4K (CAA→AAA)  fdoH → formate dehydrogenase‑O subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947578,0440CA 8.7% 138.3 / 2.1 46Q4K (CAA→AAA) fdoHformate dehydrogenase‑O subunit beta
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base A (2/2);  total (23/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 = 9.02e-01

ATCGGCATTCCGCTGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGGAATGAACTGCGTGACG  >  NC_002947/577895‑578191
                                                                                                                                                     |                                                                                                                                                   
aTCGGCATTCCGCTGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCCCGTGAAAGTGGAGAAACTCTAATGgccagcc                                                                                                                                                     >  1:292605/1‑150 (MQ=255)
    gCATTCCGCTGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGa                                                                                                                                                 >  2:8755/1‑150 (MQ=255)
     cATTCCGCTGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGAc                                                                                                                                                <  1:858554/150‑1 (MQ=255)
          cGCTGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGAcatcat                                                                                                                                           >  1:491945/1‑150 (MQ=255)
           gCTGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACatcatc                                                                                                                                          <  1:863241/150‑1 (MQ=255)
             tGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGc                                                                                                                                        <  1:1084164/150‑1 (MQ=255)
                  cggggCCTCACCGGCGCTACCCGGCATGGCTAACTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCACACACCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCt                                                                                                                                   >  2:443722/2‑150 (MQ=255)
                                          cATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCt                                                                                                           <  1:635008/150‑1 (MQ=255)
                                            tGGCTACCGGACCAATACCTTGGTGCCGTTCCTCGGTGACAGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCg                                                                                                         <  2:841296/150‑1 (MQ=255)
                                                      aCCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGt                                                                                                                   <  1:39212/130‑1 (MQ=255)
                                                      aCCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGt                                                                                                                   >  2:39212/1‑130 (MQ=255)
                                                          aTACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCagga                                                                                           <  1:587311/150‑1 (MQ=255)
                                                               ttGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCCGCaggaggagg                                                                                      >  1:27563/1‑150 (MQ=255)
                                                                 ggTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTc                                                                                    <  2:424467/150‑1 (MQ=255)
                                                                 ggTGCCGTTCCGCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGACGCCCTCGGTGCGTCAGCAGGAGGAGGTc                                                                                    <  2:468666/150‑1 (MQ=255)
                                                                        ttCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGc                                                                             <  1:146577/150‑1 (MQ=255)
                                                                         tCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCt                                                                            <  2:179584/150‑1 (MQ=255)
                                                                                    ggCACCACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTa                                                                 <  2:139845/150‑1 (MQ=255)
                                                                                    ggCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGCAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTa                                                                 >  1:202431/1‑150 (MQ=255)
                                                                                        acacGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCg                                                                                                                    <  2:957344/95‑1 (MQ=255)
                                                                                        acacGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCg                                                                                                                    >  1:957344/1‑95 (MQ=255)
                                                                                        acacGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCaa                                                             <  1:8755/150‑1 (MQ=255)
                                                                                             cAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCa                                                        >  1:572061/1‑150 (MQ=255)
                                                                                              aGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCAt                                                       >  2:1085206/1‑150 (MQ=255)
                                                                                                  gCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGGGGACGGCAAGCTGATCGACACTACCAAGTGCATCGGc                                                   >  2:640565/1‑150 (MQ=255)
                                                                                                  gCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTc                                                                                    >  1:1064816/1‑117 (MQ=255)
                                                                                                  gCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTc                                                                                    <  2:1064816/117‑1 (MQ=255)
                                                                                                  gCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGc                                                   <  2:491945/150‑1 (MQ=255)
                                                                                                    cAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTg                                                 >  1:414381/1‑150 (MQ=255)
                                                                                                    cAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTg                                                 >  1:128778/1‑150 (MQ=255)
                                                                                                            aaGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCt                                         >  2:552070/1‑150 (MQ=255)
                                                                                                               tCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGcc                                      >  2:199717/1‑150 (MQ=255)
                                                                                                                  ttCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGc                                                                             >  1:916510/1‑108 (MQ=255)
                                                                                                                  ttCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGc                                                                             <  2:916510/108‑1 (MQ=255)
                                                                                                                      tCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCTCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGc                               <  2:1122229/150‑1 (MQ=255)
                                                                                                                      tCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGc                               >  1:1122229/1‑150 (MQ=255)
                                                                                                                        gTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCgg                             >  1:549435/1‑149 (MQ=255)
                                                                                                                           aaaGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCt                          >  1:100843/1‑150 (MQ=255)
                                                                                                                               tGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGa                      <  2:100843/150‑1 (MQ=255)
                                                                                                                                    aaaCTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGGa                 >  1:173677/1‑150 (MQ=255)
                                                                                                                                    aaaCTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGGa                 >  2:9022/1‑150 (MQ=255)
                                                                                                                                      aCTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGGAAt               <  2:173677/150‑1 (MQ=255)
                                                                                                                                         ctGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGTAATGaa            <  1:1095158/150‑1 (MQ=255)
                                                                                                                                          tGATGGCCAGCAAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCg                        <  2:958823/137‑1 (MQ=255)
                                                                                                                                          tGATGGCCAGCAAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCg                        <  2:961975/137‑1 (MQ=255)
                                                                                                                                          tGATGGCCAGCAAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCg                        >  1:961975/1‑137 (MQ=255)
                                                                                                                                          tGATGGCCAGCAAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCg                        >  1:958823/1‑137 (MQ=255)
                                                                                                                                                 cagccaAGACATCATCGCCCGCTCGGCAACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGGAATGAACTGCGTGa    >  2:404229/1‑150 (MQ=255)
                                                                                                                                                   gccaAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGGAATGAACTGCGTGACg  >  2:264960/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATCGGCATTCCGCTGCACTGGGGCTTCACCGGCGCTACCCGGCATGGCTACCTGACCAATACCTTGGTGCCGTTCCTCGGTGACGGCAACACGCAGACGCCAGAGTCCAAGTCGTTCCTCGTGAAAGTGGAGAAACTCTGATGGCCAGCCAAGACATCATCGCCCGCTCGGCCACCACTACCGTGCCGCCCTCGGTGCGTCAGCAGGAGGAGGTCGCCAAGCTGATCGACACTACCAAGTGCATCGGCTGCAAGGCCTGCCAGGTGGCGTGCTCGGAGTGGAATGAACTGCGTGACG  >  NC_002947/577895‑578191

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