Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,628,125 C→G 23.3% 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.3% 85.4 / 25.3 43intergenic (+11/‑129)PP_3198/PP_31992Fe‑2S ferredoxin/monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (15/18);  new base G (4/6);  total (19/24)
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.86e-01

ACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTG  >  NC_002947/3627984‑3628265
                                                                                                                                             |                                                                                                                                            
aCGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTGCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATa                                                                                                                                      <  2:701378/150‑1 (MQ=255)
 cGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAg                                                                                                                                     <  2:333961/150‑1 (MQ=255)
       gcgcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCa                                                                                                                               <  1:152011/150‑1 (MQ=255)
           ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAg                                                                                                                                     >  1:529948/1‑140 (MQ=255)
           ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAg                                                                                                                                     <  2:529948/140‑1 (MQ=255)
               cGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcga                                                                                                                       <  2:322955/150‑1 (MQ=255)
                gCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgac                                                                                                                      <  2:656008/150‑1 (MQ=255)
                  ccGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCg                                                                                                                                 >  2:220064/1‑137 (MQ=255)
                  ccGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCg                                                                                                                                 <  1:220064/137‑1 (MQ=255)
                          ggCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTc                                                                                                            >  1:646268/1‑150 (MQ=255)
                            cACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCAt                                                                                                          <  1:158151/150‑1 (MQ=255)
                            cACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCAt                                                                                                          <  1:538793/150‑1 (MQ=255)
                                         ggCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCCGACCTATCGCGCGGCGCGATAg                                                                                                                                     <  2:544908/110‑1 (MQ=255)
                                         ggCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCCGACCTATCGCGCGGCGCGATAg                                                                                                                                     >  1:544908/1‑110 (MQ=255)
                                            cccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCTCGCGACGCGATTGGTCATGGACCGGCGCAGGCgg                                                                                          <  1:550100/150‑1 (MQ=255)
                                                  cGAGGCGAAGCCTTGGAGGCACTCGATCACAAGTAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCt                                                                                    <  2:699800/150‑1 (MQ=255)
                                                            ccTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGc                                                                          >  1:158754/1‑150 (MQ=255)
                                                                          gATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTgg                                                                                                              <  2:785256/100‑1 (MQ=255)
                                                                          gATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTgg                                                                                                              >  1:785256/1‑100 (MQ=255)
                                                                            tCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgccgcc                                                          >  2:164165/1‑150 (MQ=255)
                                                                               caAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCg                                                       >  2:240419/1‑150 (MQ=255)
                                                                               caAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCg                                                       <  2:295905/150‑1 (MQ=255)
                                                                                     aGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAg                                                                                                                                     <  2:624894/66‑1 (MQ=255)
                                                                                     aGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAg                                                                                                                                     >  1:624894/1‑66 (MQ=255)
                                                                                        gggcggCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACggg                                              >  1:777428/1‑150 (MQ=255)
                                                                                                 gctgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTGGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGCAAGCTGTGCATAGGCACAGGAAACGGCCGCACCGCATGACGGGGCGCacaac                                     >  2:38917/1‑150 (MQ=255)
                                                                                                   tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCa                                                                               >  2:678781/1‑106 (MQ=255)
                                                                                                   tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCa                                                                               <  1:678781/106‑1 (MQ=255)
                                                                                                   tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCAcaacaa                                   <  1:591623/150‑1 (MQ=255)
                                                                                                     ctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGa                                 >  2:561556/1‑150 (MQ=255)
                                                                                                       aTTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACgggg                                             <  2:224038/136‑1 (MQ=255)
                                                                                                       aTTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACgggg                                             >  1:224038/1‑136 (MQ=255)
                                                                                                            cGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAAc                          <  2:646268/150‑1 (MQ=255)
                                                                                                                tGAGCAGCTCGATAGCTGCCCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                      <  1:164165/150‑1 (MQ=255)
                                                                                                                   gcagcTAGATCGCTAACCTTTCGCGCCGCTCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCAGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGAc                   >  1:220095/1‑150 (MQ=255)
                                                                                                                       cTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGc               <  2:158754/150‑1 (MQ=255)
                                                                                                                          gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATg                                                                                                         >  1:783389/1‑57 (MQ=255)
                                                                                                                          gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATg                                                                                                         <  2:783389/57‑1 (MQ=255)
                                                                                                                                 gACCTATCGCGCGGCGCGATAGGTCACATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACt     >  2:161825/1‑150 (MQ=255)
                                                                                                                                   ccTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCa             <  1:596537/140‑1 (MQ=255)
                                                                                                                                   ccTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCa             >  2:596537/1‑140 (MQ=255)
                                                                                                                                   ccTATCGCGCGGCGCGATAGGTCACATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt   <  1:161825/150‑1 (MQ=255)
                                                                                                                                    cTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgtg  >  2:621784/1‑150 (MQ=255)
                                                                                                                                             cgcgcgATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACt     >  1:597286/1‑138 (MQ=255)
                                                                                                                                             cgcgcgATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACt     <  2:597286/138‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
ACGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTG  >  NC_002947/3627984‑3628265

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