Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,628,125 C→G 23.9% intergenic (+11/‑129) PP_3198 → / → PP_3199 2Fe‑2S ferredoxin/monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,628,1250CG23.9% 32.2 / 9.7 21intergenic (+11/‑129)PP_3198/PP_31992Fe‑2S ferredoxin/monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base G (3/2);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGA  >  NC_002947/3627985‑3628266
                                                                                                                                            |                                                                                                                                             
cGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAg                                                                                                                                      <  2:188351/150‑1 (MQ=255)
cGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGCCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAg                                                                                                                                      <  1:140694/150‑1 (MQ=255)
                 ccGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgacgc                                                                                                                     >  2:185458/1‑150 (MQ=255)
                          gCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCa                                                                                                            <  1:185458/150‑1 (MQ=255)
                                          ccccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                            <  1:120456/150‑1 (MQ=255)
                                              tGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTg                                                                                                                <  1:274151/126‑1 (MQ=255)
                                              tGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTg                                                                                                                >  2:274151/1‑126 (MQ=255)
                                                      cGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCa                                                                                >  1:186553/1‑150 (MQ=255)
                                                         aGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAg                                                                             <  2:186553/150‑1 (MQ=255)
                                                                        cGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgcc                                                              <  1:246996/150‑1 (MQ=255)
                                                                           tCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATa                                                                                                                                       >  2:412686/1‑74 (MQ=255)
                                                                           tCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATa                                                                                                                                       <  1:412686/74‑1 (MQ=255)
                                                                                     gcgggcggCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACg                                                 >  2:256884/1‑150 (MQ=255)
                                                                                        ggcggcTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACgggg                                              <  1:43156/150‑1 (MQ=255)
                                                                                                   gctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAg                                   >  2:156873/1‑150 (MQ=255)
                                                                                                             cTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAg                                                                             >  2:190357/1‑98 (MQ=255)
                                                                                                             cTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAg                                                                             <  1:190357/98‑1 (MQ=255)
                                                                                                              ttGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCNCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAGCAAAACGTg                        >  1:274663/1‑150 (MQ=255)
                                                                                                                                   ccaTCGCGCCCCGCAATAGGTCGCATCCCCCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgtg   <  2:274663/148‑1 (MQ=255)
                                                                                                                                   cTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgtg   >  1:370690/1‑150 (MQ=255)
                                                                                                                                    tATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGAAAGGTGGGCAGTTGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGa  >  2:366137/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                             
CGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGA  >  NC_002947/3627985‑3628266

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