Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,239,707 G→C 13.4% G101A (GGG→GCG)  PP_1080 → molybdopterin‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,239,7070GC13.4% 72.2 / 6.1 30G101A (GGG→GCG) PP_1080molybdopterin‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (11/15);  new base C (2/2);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

CAAGGCCGTGGCGTTGCAGGACATCCAGAACGACCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAACCGTTTCTCAGCCACTTCAGTGGACCAGTTGCCCCAACACCTGGTCA  >  NC_002947/1239558‑1239853
                                                                                                                                                     |                                                                                                                                                  
cAAGGCCGTGGCGTTGCAGGACATCCAGAACGACCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATgg                                                                                                                                                    >  2:463245/1‑150 (MQ=255)
             ttGCAGGACATCCAGAACGACCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTg                                                                                                                                       <  2:81374/150‑1 (MQ=255)
             ttGCAGGACATCCAGAACGACCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTg                                                                                                                                       <  2:81200/150‑1 (MQ=255)
              gatctgACATCCAGAACGACCCGGGTCGCCTGCGCCAGCCCCACAAGCGGGTCGGTGAGCATTGGCAGGCTATTGGCTCGGATGAGGCTTTCGCACTGTCAGCCGATGAATGGGGTGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTgg                                                                                                                                      <  1:320669/145‑1 (MQ=255)
                    aCATCCAGAACGACCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGcagatc                                                                                                                                >  2:320669/1‑145 (MQ=255)
                         cAGAACGACCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCaa                                                                                                                           >  1:305223/1‑150 (MQ=255)
                                aCCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAAcc                                                                                                                    <  1:212260/150‑1 (MQ=255)
                                  ccGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCg                                                                                                                  >  1:306006/1‑150 (MQ=255)
                                            tGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGcaca                                                                                                        >  1:597307/1‑150 (MQ=255)
                                                 cAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTAc                                                                                                   >  2:416835/1‑150 (MQ=255)
                                                    ccccACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGc                                                                                                <  2:305223/150‑1 (MQ=255)
                                                            gCGGGTCGTCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGAc                                                                                        <  1:490889/150‑1 (MQ=255)
                                                                               cAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAg                                                                                                                          >  1:587012/1‑97 (MQ=255)
                                                                               cAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAg                                                                                                                          <  2:587012/97‑1 (MQ=255)
                                                                                    tATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCt                                                                >  1:192005/1‑150 (MQ=255)
                                                                                          cTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCaa                                                          <  2:549670/150‑1 (MQ=255)
                                                                                           tgagATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGCGCGCAATGCCGTg                                                                                                                                       <  2:180006/69‑1 (MQ=255)
                                                                                           tgagATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGCGCGCAATGCCGTg                                                                                                                                       >  1:180006/4‑72 (MQ=255)
                                                                                           tGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGCGCGCAATGCCGTg                                                                                                                                       <  2:180049/72‑1 (MQ=255)
                                                                                           tGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGCGCGCAATGCCGTg                                                                                                                                       >  1:180049/1‑72 (MQ=255)
                                                                                            gggATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGa                                                        <  2:532919/150‑1 (MQ=255)
                                                                                                 gAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGt                                                   >  1:239958/1‑150 (MQ=255)
                                                                                                    gCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCATCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAAc                                                >  1:171052/1‑150 (MQ=255)
                                                                                                             atGGCAGCGGAAAAATTGTGGGCGTTCCAGCAAGCCCAGGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAACCGTTTCTCa                                       <  2:239958/149‑1 (MQ=255)
                                                                                                                  aGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAACCGTTTCTCAGCCAc                                  >  2:24591/1‑150 (MQ=255)
                                                                                                                            ttGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGCCCCGTAACCGTTTCTCAGCCACTTCAGTGGAc                        <  1:182648/150‑1 (MQ=255)
                                                                                                                                 ggCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAACCGTTTCTCAGCCACTTCAGTGGACCAGtt                   <  2:597307/150‑1 (MQ=255)
                                                                                                                                                  aTGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAACCGTTTCTCAGCCACTTCAGTGGACCAGTTGCCCCAACACCTGGTCa  <  2:171052/150‑1 (MQ=255)
                                                                                                                                                  aTGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAACCGTTTCTCAGCCACTTCAGTGGACCAGTTGCCCCAACACCTGGTCa  <  2:552510/150‑1 (MQ=255)
                                                                                                                                                   tGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGATCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTaa                                                 >  1:502463/1‑102 (MQ=255)
                                                                                                                                                   tGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTaa                                                 <  2:502463/102‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAAGGCCGTGGCGTTGCAGGACATCCAGAACGACCCGGATCGCCTGCGCCAGCCCCACAAGCGGGTCGGCGAGCATTGGCAGGCTATTGGCTGGGATGAGGCTTTCGCACTGGCAGCCGAAAAATTGTGGGCGGTCCAGCAAGCCCATGGGCGCAATGCCGTGGCGGTCTACCAAGGCAACCCGAGCGTGCACAATTACGGCTTGATGACCCACAGCAACTACTTTCTAGGCCTGCTCAAGACCCGTAACCGTTTCTCAGCCACTTCAGTGGACCAGTTGCCCCAACACCTGGTCA  >  NC_002947/1239558‑1239853

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: ‑