Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 435,616 C→A 13.8% intergenic (‑414/‑36) glcB ← / → PP_0357 malate synthase G/ACT domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947435,6160CA13.8% 92.0 / 2.6 29intergenic (‑414/‑36)glcB/PP_0357malate synthase G/ACT domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (15/10);  new base A (2/2);  total (17/12)
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.05e-01

GGGAGTAGGTATGTTCGCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTGGTAACTGGCT  >  NC_002947/435484‑435746
                                                                                                                                    |                                                                                                                                  
gggAGTAGGTATGTTCGCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTgc                                                                                                                                 >  1:440396/1‑136 (MQ=255)
          aTGTTCGCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGgcgc                                                                                                                       >  2:253054/1‑136 (MQ=255)
              tCGCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGACCGGCGCGGTc                                                                                                                   >  1:57581/1‑136 (MQ=255)
               cGCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTcc                                                                                                                  <  1:190581/136‑1 (MQ=255)
                gCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCg                                                                                                                 >  1:253780/1‑136 (MQ=255)
                gCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCg                                                                                                                 >  2:242165/1‑136 (MQ=255)
                gCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTATGCCGCCCGGCGCGGTcc                                                                                                                  >  2:469491/1‑135 (MQ=255)
                gCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTATGCCGCCCGGCGCGGTcc                                                                                                                  >  2:306857/1‑135 (MQ=255)
                gCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTATGCCGCCCGGCGCGGTcc                                                                                                                  <  1:306857/135‑1 (MQ=255)
                gCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTATGCCGCCCGGCGCGGTcc                                                                                                                  <  1:469491/135‑1 (MQ=255)
                      aaCAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGa                                                                                                           <  2:253780/136‑1 (MQ=255)
                        cAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAgg                                                                                                         <  2:81678/136‑1 (MQ=255)
                                            aaTCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCtt                                                                                     <  2:418892/136‑1 (MQ=255)
                                                       aGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGt                                                                          <  1:253054/136‑1 (MQ=255)
                                                            acGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAAttt                                                                     >  1:397702/1‑135 (MQ=255)
                                                                      gtgCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACaa                                                           >  1:416687/1‑136 (MQ=255)
                                                                      gtgCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACaa                                                           >  1:332417/1‑136 (MQ=255)
                                                                           gATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAggccg                                                      <  1:12967/136‑1 (MQ=255)
                                                                                       gCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAAc                                          <  1:242165/136‑1 (MQ=255)
                                                                                        cAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTAGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACg                                         >  1:390676/1‑136 (MQ=255)
                                                                                                                     cGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTAtt                                     >  2:38496/1‑111 (MQ=255)
                                                                                                                     cGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTAtt                                     <  1:38496/111‑1 (MQ=255)
                                                                                                                     cGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTg            >  1:193963/1‑136 (MQ=255)
                                                                                                                     cGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTg            >  2:480293/1‑136 (MQ=255)
                                                                                                                      gATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTgg           <  2:332417/136‑1 (MQ=255)
                                                                                                                      gATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTgg           <  2:416687/136‑1 (MQ=255)
                                                                                                                      gATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTgg           >  2:430527/1‑136 (MQ=255)
                                                                                                                          ccccGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTGGTAAc       >  1:171441/1‑136 (MQ=255)
                                                                                                                               aTTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTGGTAActggct  >  2:348753/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GGGAGTAGGTATGTTCGCAAAAAACAGGTGGATTGTTCCAGATAAATCTTGAAACAGTACACGATTGGCTGTGCAGATTGCCCTGCGGCAGGTTGCCGCGCCTTGCCTATACTGACTCGATCCCCCGATTTTCTGCCGCCCGGCGCGGTCCGACACGAGGAAGGGCTTGTGGATCATCTTGTACTGACTGTAATTGCCCCTGACAAGGCCGGCCAGGTCGAACGTATTGCCCAGTGCATCGCCGACCACAGTGGTAACTGGCT  >  NC_002947/435484‑435746

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