Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,980,326 C→G 8.8% D117H (GAC→CAC) ‡ cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,980,3260CG 8.8% 133.4 / 3.3 45D117H (GAC→CAC) ‡cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base C (23/18);  new base G (2/2);  total (25/20)
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.68e-01

TCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCG  >  NC_002947/3980184‑3980475
                                                                                                                                              |                                                                                                                                                      
tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                                                                                                                                                 <  2:352028/150‑1 (MQ=255)
tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGACAGCTCCGGTGCCGGGCGATGATCGccc                                                                                                                                                 >  2:677858/1‑150 (MQ=255)
  gcgGCCCAGCCATTGGCTGGCCAGGAAGTTGATCCCGTGTCCAGTCTTCGGCTTTCCGCCCTGGCGCAGATAATGCTAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGATGACCGCCCgg                                                                                                                                               <  1:405321/150‑1 (MQ=255)
     gCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCAc                                                                                                                                            >  1:640317/1‑150 (MQ=255)
       ccAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGATGACCGCCCgga                                                                                                                                              >  1:445340/1‑145 (MQ=255)
       ccAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGATGACCGCCCgga                                                                                                                                              <  2:445340/146‑2 (MQ=255)
        cAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAg                                                                                                                                         >  1:159614/1‑150 (MQ=255)
        cAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAg                                                                                                                                         >  1:710800/1‑150 (MQ=255)
             aTTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCa                                                                                                                                    >  1:269108/1‑150 (MQ=255)
               tggctggcCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGc                                                                                                                                  <  2:235389/150‑1 (MQ=255)
                        aGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCa                                                                                                                        <  2:116291/150‑1 (MQ=255)
                              ttGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGggcg                                                                                                                  <  1:526581/150‑1 (MQ=255)
                                   cAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAg                                                                                                             >  1:185943/1‑150 (MQ=255)
                                     ggTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCt                                                                                                           >  2:149358/1‑150 (MQ=255)
                                      gTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTc                                                                                                          <  2:185943/150‑1 (MQ=255)
                                               gTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgttg                                                                                                 <  2:640317/150‑1 (MQ=255)
                                                   gCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCACAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCCCCCGAATCACCTGCACGCCACGGGCTGCCAGCTCGACCAGTTGTTCg                                                                                             >  2:94513/1‑150 (MQ=255)
                                                     cTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCg                                                                                                                                                <  1:197430/98‑1 (MQ=255)
                                                     cTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCg                                                                                                                                                >  2:197430/1‑98 (MQ=255)
                                                     cTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGAc                                                                                           >  1:588882/1‑150 (MQ=255)
                                                     cTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGAc                                                                                           >  2:215431/1‑150 (MQ=255)
                                                      ttGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCa                                                                                                                                    <  2:334951/109‑1 (MQ=255)
                                                      ttGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCa                                                                                                                                    >  1:334951/1‑109 (MQ=255)
                                                           gccCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTa                                                                                     <  1:215431/150‑1 (MQ=255)
                                                            cacTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGTTCCGGGCGATCATCGCCCCACAACTGAATAGCGTACTCTCAACGGCCGGCCAGCTCGACCAGTTGTTCGACGcctgag                                                                                    >  2:22561/3‑146 (MQ=255)
                                                             ccTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAgc                                                                                   >  2:157094/1‑150 (MQ=255)
                                                             ccTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAgc                                                                                   <  2:588882/150‑1 (MQ=255)
                                                                  cgcAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAg                                                                              <  1:157094/150‑1 (MQ=255)
                                                                         aaTGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgtt                                                                                                  <  2:65292/123‑1 (MQ=255)
                                                                         aaTGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAgtt                                                                                                  >  1:65292/1‑123 (MQ=255)
                                                                                  gCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGt                                                              >  1:470645/1‑150 (MQ=255)
                                                                                            cTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGacac                                                    >  1:322436/1‑150 (MQ=255)
                                                                                              gCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCa                                                  >  1:156025/1‑150 (MQ=255)
                                                                                                tCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAgg                                                >  1:735478/1‑150 (MQ=255)
                                                                                                  aCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCCGCAGCTCGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                              <  1:149358/150‑1 (MQ=255)
                                                                                                  aCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTc                                                                                              >  1:693599/1‑102 (MQ=255)
                                                                                                  aCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTc                                                                                              <  2:693599/102‑1 (MQ=255)
                                                                                                  aCCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGAt                                              >  1:265071/1‑150 (MQ=255)
                                                                                                               cccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                                 <  2:322436/150‑1 (MQ=255)
                                                                                                               cccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                                 >  2:447573/1‑150 (MQ=255)
                                                                                                               cccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtgg                                 <  2:841429/150‑1 (MQ=255)
                                                                                                                ccAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTtggcgtggc                                >  1:428117/1‑150 (MQ=255)
                                                                                                                   ggCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACCCCACGCTGACCATGGCGCGGCgca                             >  2:266081/1‑150 (MQ=255)
                                                                                                                                    tccgggCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATAc            <  1:628209/150‑1 (MQ=255)
                                                                                                                                            gTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGAc    <  1:276075/150‑1 (MQ=255)
                                                                                                                                              cATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCg  >  1:549081/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
TCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAGCTGCACGCCACGGGCGGCCAGCTCGACCAGTTGTTCGACGCCGTAGCGCCAGTAGCCAACGCCACCGTGCAGCGACACCAGGATGACCTTGGCGTGGCGCAGCACCTGGTCGACATACAGGTCGACCG  >  NC_002947/3980184‑3980475

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