Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,295 C→A 21.1% R233S (CGT→AGT)  PP_0159 → family III CoA‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947170,2950CA21.1% 31.5 / 4.6 19R233S (CGT→AGT) PP_0159family III CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base C (9/6);  new base A (2/2);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GGTTGTATTCCGCGGTCGCCATGCTGGCTGCCCTCGCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTCGGTAATGACAGCCAGTTCCGCAAGTTCGCTGAGGTGGCCGGGCAGCCGCAATGGGCGGAG  >  NC_002947/170146‑170435
                                                                                                                                                       |                                                                                                                                            
ttctGTATTCCGCGGTCGCCATGCTGGCTGCCCTCGCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTcccc                                                                                                                                              <  2:596881/147‑1 (MQ=255)
     tATTCCGCGGTCGCCATGCTGGCTGCCCTCGAACATCGTGATCATACAGGGGGCGGCCCCGCTGTCGACATGGCGTCGATCCATTCGCGCTCGGGCACGTTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCgagat                                                                                                                                         >  1:596881/1‑146 (MQ=255)
                  ccATGCTGGCTGCCCTCGCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATgcgc                                                                                                                            >  1:244074/1‑150 (MQ=255)
                  ccATGCTGGCTGCCCTCGCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATgcgc                                                                                                                            <  2:491715/150‑1 (MQ=255)
                                 tCGCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTcc                                                                                                             >  1:520777/1‑150 (MQ=255)
                                   gCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCAGTCgg                                                                                                                                         >  1:80394/1‑120 (MQ=255)
                                   gCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCAGTCgg                                                                                                                                         <  2:80394/120‑1 (MQ=255)
                                   gCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCAGTCgg                                                                                                                                         <  2:80377/120‑1 (MQ=255)
                                   gCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCAGTCgg                                                                                                                                         >  1:80377/1‑120 (MQ=255)
                                                             ccAGCATGTCGACATGGCATTGCTTGATGTGGAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATgg                                                                                 <  2:40419/150‑1 (MQ=255)
                                                                       gACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCAtt                                                                       >  2:141812/1‑150 (MQ=255)
                                                                                  gCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCCCCCCACTACGTCTGTCAGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTCgg                                                            <  1:380732/150‑1 (MQ=255)
                                                                                            cAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTCGGTAATGAAAGc                                                  >  1:196386/1‑150 (MQ=255)
                                                                                                             gCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACGAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTAGGTAATGACAGCCAGTTCCGCAAGTTc                                   >  2:223939/1‑150 (MQ=255)
                                                                                                              cGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTCGGTAATGACAGCCAGTTCCGCAAGTTCg                                  <  1:293015/150‑1 (MQ=255)
                                                                                                                    aGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATgc                                                                                                                              <  1:59009/52‑1 (MQ=255)
                                                                                                                    aGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATgc                                                                                                                              >  2:59009/1‑52 (MQ=255)
                                                                                                                           gAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTCGGTAATGACAGCCAGTTCCGCAAGTTCGCTGAGGTGGCCgg                     >  1:4455/1‑150 (MQ=255)
                                                                                                                                            gCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTCGGTAATGACAGCCAGTTCCGCAAGTTCGCTGAGGTGGCCGGGCAGCCGCAATGGGCgg    >  2:484989/1‑150 (MQ=255)
                                                                                                                                              aCCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGTCAGATGGCGATTTCATTCTTACCGTCGGTAATGACAGCCAGTTCCGCAAGTTCGCTGAGGTGGCCGGGCAGCCGCAATGGGCggag  <  1:18441/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                            
GGTTGTATTCCGCGGTCGCCATGCTGGCTGCCCTCGCCCATCGGGATCAGACAGGGGTCGGCCAGCATGTCGACATGGCATTGCTTGATGTGCAGGT‑‑GGCCTGTCTGGCGAACCAGGCGATGAACTACCTGACCACGGGCACCCCTCCCCGTCGGCTGGGTAATGCGCATCCTAATATAGTCCCTTACCAGGACTTTCCGACGGCAGATGGCGATTTCATTCTTACCGTCGGTAATGACAGCCAGTTCCGCAAGTTCGCTGAGGTGGCCGGGCAGCCGCAATGGGCGGAG  >  NC_002947/170146‑170435

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