Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,423,575 C→G 14.3% P564A (CCA→GCA)  PP_3038 → TraC domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,423,5750CG14.3% 73.3 / 6.5 28P564A (CCA→GCA) PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base G (2/2);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TTCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAGGCCAT  >  NC_002947/3423431‑3423723
                                                                                                                                                |                                                                                                                                                    
ttCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGGCACGt                                                                                                                                                 <  1:368990/150‑1 (MQ=255)
           ccGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGtttt                                                                                                                                      >  1:358921/1‑150 (MQ=255)
                                  tGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATg                                                                                                               <  1:82496/150‑1 (MQ=255)
                                   gggCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCCAGTCAGTGTGTATTGGCGGGTACCACCAACCAGGATGa                                                                                                              >  2:65102/1‑150 (MQ=255)
                                       gAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATAc                                                                                                          <  1:46100/150‑1 (MQ=255)
                                        aaTTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGGGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATAcc                                                                                                         >  1:316221/1‑150 (MQ=255)
                                                  cTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGc                                                                                               <  2:358921/150‑1 (MQ=255)
                                                       acaaGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCg                                                                                          >  2:344979/1‑150 (MQ=255)
                                                         aaGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGt                                                                                        >  1:14176/1‑150 (MQ=255)
                                                              agagagCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCg                                                                                   <  2:12255/150‑1 (MQ=255)
                                                                  agCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGt                                                                               <  2:456112/150‑1 (MQ=255)
                                                                      ccAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTAtt                                                                           <  2:2530/150‑1 (MQ=255)
                                                                           gCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCg                                                                      <  2:222626/150‑1 (MQ=255)
                                                                                  aaTTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGc                                                                                 >  1:452524/1‑132 (MQ=255)
                                                                                  aaTTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGc                                                                                 <  2:452524/132‑1 (MQ=255)
                                                                                     tcttctCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCg                                                            >  1:348107/1‑150 (MQ=255)
                                                                                       ttctCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTcgc                                                                                                                                   <  2:81933/77‑1 (MQ=255)
                                                                                       ttctCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTcgc                                                                                                                                   >  1:81933/1‑77 (MQ=255)
                                                                                                 cTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGAcc                                                <  2:348107/150‑1 (MQ=255)
                                                                                                       acacTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAg                                          >  2:203368/1‑150 (MQ=255)
                                                                                                                          tGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGa                       >  1:16197/1‑150 (MQ=255)
                                                                                                                               gaagaaCGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTaccac                                                                                                                          <  2:147618/46‑1 (MQ=255)
                                                                                                                               gaagaaCGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTaccac                                                                                                                          >  1:147618/1‑46 (MQ=255)
                                                                                                                               gaagaaCGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGc                  >  2:345677/1‑150 (MQ=255)
                                                                                                                                       actGACGTGGCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGCAGCATTGAGTCGCGCTCGTGACCAGCTGTggggt          >  2:208645/3‑148 (MQ=255)
                                                                                                                                       actGACGTGGCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCt          >  2:76063/3‑150 (MQ=255)
                                                                                                                                         tGACGTGGCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGa        <  1:208645/150‑1 (MQ=255)
                                                                                                                                              tGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAGGCCAt  <  1:404981/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTCCAAGCCATCCGAGGTAAGTGGATCGTCGAGCTGGGCGAATTGGACAGCTTCAACAAGGCAGAGAGCACCAAGGCCAAGCAATTCTTCTCGGCATCTATCGACACTTACCGCGAAAGCTATGGGCGAAGAACGAGTGACGTGCCACGTCAGTGTGTTTTCGCGGGTACCACCAACCAGGATGAATACCTTAAGGACGCCACCGGTAACCGGCGTTATTGGCCGGTCGCTTGCGTCAAGGTGGACCTTGAAGCATTGCGTCGCGCTCGTGACCAGCTGTGGGCTGAGGCCAT  >  NC_002947/3423431‑3423723

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: