Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 461,983 C→T 9.7% G162D (GGC→GAC)  pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947461,9830CT 9.7% 107.1 / 3.3 41G162D (GGC→GAC) pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base C (16/21);  new base T (2/2);  total (18/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 = 4.68e-01

CCTGCTGCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCGG  >  NC_002947/461834‑462131
                                                                                                                                                     |                                                                                                                                                    
ccTGCTGCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGc                                                                                                                                                      <  1:929629/150‑1 (MQ=255)
 ctgctgCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATGTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGCCCGCCTCGTGCAAGGCCAGGGTGTATGGATTGcc                                                                                                                                                     >  2:549572/1‑150 (MQ=255)
                         ttCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATg                                                                                                                             >  2:716680/1‑150 (MQ=255)
                              ccAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCCTGAAAAGCGCCCAGTGGATGCCTgg                                                                                                                        >  2:104306/1‑150 (MQ=255)
                              ccAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTgg                                                                                                                        <  1:146797/150‑1 (MQ=255)
                                   ggCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGac                                                                                                                   >  1:313996/1‑150 (MQ=255)
                                      tGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGACAAGCGCCCAGTGGGTGCCTGGCTGACact                                                                                                                >  1:602901/1‑150 (MQ=255)
                                         ggACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCa                                                                                                                                                 <  2:471867/114‑1 (MQ=255)
                                         ggACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCa                                                                                                                                                 >  1:471867/1‑114 (MQ=255)
                                          gACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACactgact                                                                                                            <  1:529099/150‑1 (MQ=255)
                                             cGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCATGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgca                                                                                                         <  1:1069660/150‑1 (MQ=255)
                                             cGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATTAAAAGCGCCCAGGGGATGCCTGGCTGACACTGACGgca                                                                                                         >  1:269925/1‑150 (MQ=255)
                                               cAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgcagc                                                                                                       >  2:1031203/1‑150 (MQ=255)
                                               cAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgcagc                                                                                                       >  2:1031124/1‑150 (MQ=255)
                                                aaaGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTgcagca                                                                                                      >  1:318300/1‑150 (MQ=255)
                                                 aaGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGACAGCTAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCaa                                                                                                     <  1:996526/150‑1 (MQ=255)
                                                   ggCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCGCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCCCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAgg                                                                                                   <  1:1037496/150‑1 (MQ=255)
                                                      tCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGa                                                                                                >  1:610424/1‑150 (MQ=255)
                                                        aGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAAc                                                                                              <  1:549572/150‑1 (MQ=255)
                                                          gCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCGGCGAAAGCCTGGCTGACACTGACTGCAGCAAGGCGAActg                                                                                            >  1:66977/1‑150 (MQ=255)
                                                                       ccTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCACGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGtt                                                                               <  1:1031203/150‑1 (MQ=255)
                                                                       ccTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGATTGCGGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGTCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGtt                                                                               <  1:1031124/150‑1 (MQ=255)
                                                                            gTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGa                                                                                                                               <  2:1201161/97‑1 (MQ=255)
                                                                            gTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGa                                                                                                                               >  1:1201161/1‑97 (MQ=255)
                                                                              aGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTgc                                                                        <  2:792784/150‑1 (MQ=255)
                                                                                   cGTTGATGAAAGACCGCAAGTGTTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGAc                                                                   <  2:66977/150‑1 (MQ=255)
                                                                                        atgaAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGa                                                                                                                    >  1:35233/1‑96 (MQ=255)
                                                                                        atgaAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGa                                                                                                                    <  2:35233/96‑1 (MQ=255)
                                                                                         tgaAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGTAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCg                                                             >  1:300684/1‑150 (MQ=255)
                                                                                                   cAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGACTGCGCGACCAGGCGATGAATTCgg                                                   <  2:610424/150‑1 (MQ=255)
                                                                                                            gCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCa                                          <  1:192745/150‑1 (MQ=255)
                                                                                                              tGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACt                                        <  2:313996/150‑1 (MQ=255)
                                                                                                                     ccGCATCGTGCAAGGCAAGCGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCATGGCTGACCCTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGtt                                 <  2:269925/150‑1 (MQ=255)
                                                                                                                       gCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGTGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCa                               <  1:1068634/150‑1 (MQ=255)
                                                                                                                             tGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgt                         <  2:339437/150‑1 (MQ=255)
                                                                                                                                aaGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgttgg                      >  1:110281/1‑150 (MQ=255)
                                                                                                                                     cAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgt                         <  2:1113750/142‑1 (MQ=255)
                                                                                                                                     cAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACggcgt                         >  1:1113750/1‑142 (MQ=255)
                                                                                                                                         gtgtATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATg             >  1:1089860/1‑150 (MQ=255)
                                                                                                                                                 aTTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATc     >  1:794690/1‑150 (MQ=255)
                                                                                                                                                  ttGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATcc    <  1:803950/150‑1 (MQ=255)
                                                                                                                                                    gCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCgg  <  2:1089860/150‑1 (MQ=255)
                                                                                                                                                    gCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCgg  <  1:93560/150‑1 (MQ=255)
                                                                                                                                                     ccAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGAt              <  2:1234499/137‑1 (MQ=255)
                                                                                                                                                     ccAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGAt              >  1:1234499/1‑137 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGCTGCCGGGCAAGGCGTTCAAGTTCATCCAGCGGCTGCGGACCGCAAAGGCTCAGGCAGATCTGGCCGCCTTGGTAGAAGCGTTGATGAAAGCCCGCAAGCGCTTGCTGGAAGGCCGCATCGTGCAAGGCCAGGGTGTATCGATTGCCAGCATGAAAAGCGCCCAGTGGATGCCTGGCTGACACTGACTGCAGCAAGGCGAACTGCTGTTGGGCCGTTGGATTGCGCGACCAGGCGATGAATTCGGCATGGATCACTTCGCGTTCACGGCGTTGGCGTTCGATGCCCAGATCCGG  >  NC_002947/461834‑462131

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