Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,356,602 G→A 10.8% intergenic (+28/‑30) modC → / → PP_3831 molybdate ABC transporter ATP‑binding protein ModC/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,356,6020GA10.8% 95.7 / 2.6 37intergenic (+28/‑30)modC/PP_3831molybdate ABC transporter ATP‑binding protein ModC/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/13);  new base A (2/2);  total (22/15)
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.00e-01

CTGGATGCTGGCGGCAATGCCCTGCTTGCACGCATCACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGG  >  NC_002947/4356455‑4356746
                                                                                                                                                   |                                                                                                                                                
cTGGATGCTGGCGGCAATGCCCTGCTTGCACGCATCACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGa                                                                                                                                                 >  8:61829/1‑149 (MQ=255)
   gATGCTGGCGGCAATGCCCTGCTTGCACGCATCACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGAt                                                                                                                                                >  1:54103/1‑147 (MQ=255)
   gATGCTGGCGGCAATGCCCTGCTTGCACGCATCACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCa                                                                                                                                              >  7:26085/1‑149 (MQ=255)
           cggcAATGCCCTGCTTGCACGCATCACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTg                                                                                                                                      <  6:80465/149‑1 (MQ=255)
               aaTGCCCTGCTTGCACGCATCACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATc                                                                                                                                   >  2:114029/1‑148 (MQ=255)
               aaTGCCATGCTTGCACGCATAACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATc                                                                                                                                   >  6:236263/1‑148 (MQ=255)
                                            ctcGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGcccc                                                                                                      <  8:255357/148‑1 (MQ=255)
                                                gCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATg                                                                                                                  >  6:320565/1‑132 (MQ=255)
                                                gCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATg                                                                                                                  <  5:320565/132‑1 (MQ=255)
                                                      cAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCa                                                                                                                    <  8:73115/124‑1 (MQ=255)
                                                      cAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCa                                                                                                                    >  7:73115/1‑124 (MQ=255)
                                                      cAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTAATCAGTGACCTGATCGAGGCccgc                                                                                                                          >  7:321187/1‑118 (MQ=255)
                                                      cAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTAATCAGTGACCTGATCGAGGCccgc                                                                                                                          <  2:64157/118‑1 (MQ=255)
                                                      cAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTAATCAGTGACCTGATCGAGGCccgc                                                                                                                          <  8:321187/118‑1 (MQ=255)
                                                      cAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTAATCAGTGACCTGATCGAGGCccgc                                                                                                                          >  1:64157/1‑118 (MQ=255)
                                                            ggCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTg                                                                                       <  6:59647/147‑1 (MQ=255)
                                                            ggCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCa                                                                                     >  2:1362/1‑149 (MQ=255)
                                                                         ggCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTgacga                                                                         >  4:174014/1‑148 (MQ=255)
                                                                          gCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTgacga                                                                         >  8:252817/1‑147 (MQ=255)
                                                                                aCTGTTCGCCCAGATCAAGTCGGTGGCCCTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCa                                                                                                                                              >  5:45411/1‑71 (MQ=255)
                                                                                aCTGTTCGCCCAGATCAAGTCGGTGGCCCTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCa                                                                                                                                              <  6:45411/71‑1 (MQ=255)
                                                                                      cGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGc                                                            <  1:324398/148‑1 (MQ=255)
                                                                                      cGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGc                                                            >  2:324398/1‑148 (MQ=255)
                                                                                       gCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCt                                                          <  2:59641/149‑1 (MQ=255)
                                                                                         ccAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGa                                                        <  1:1362/149‑1 (MQ=255)
                                                                                          cAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGAc                                                       <  3:174014/149‑1 (MQ=255)
                                                                                              tCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGt                                                   >  4:138476/1‑149 (MQ=255)
                                                                                               cAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTc                                                  >  8:114744/1‑149 (MQ=255)
                                                                                                          cccTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGgcgc                                       <  7:252817/149‑1 (MQ=255)
                                                                                                           ccTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGgcgc                                       >  8:56088/1‑148 (MQ=255)
                                                                                                                                    cgcgTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGa             >  4:17369/1‑149 (MQ=255)
                                                                                                                                      cgTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATg           <  3:17369/149‑1 (MQ=255)
                                                                                                                                             tCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACg    >  3:39272/1‑149 (MQ=255)
                                                                                                                                             tCCATTGATAACTGACCTGATCGAGCCCCGCAGCCCATGCTTGACTGCCCCCTGACCCGCACCCAGCACTACGTTCACCACAGTCCGCCGGGCCAGACCCGTCGGCGCTGGCGCGCCCGTTTCATCTACCTGGCTGCCGCAGGCCAACg    >  6:267649/1‑149 (MQ=255)
                                                                                                                                                aTTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACggg  >  8:245005/1‑148 (MQ=255)
                                                                                                                                                aTTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCAGGc                                                            <  7:50783/90‑1 (MQ=255)
                                                                                                                                                aTTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCAGGc                                                            >  8:50783/1‑90 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGGATGCTGGCGGCAATGCCCTGCTTGCACGCATCACCCGCTTCTCGGCAGACCAGCTCGGCTTGCACACAGGGCAGATACTGTTCGCCCAGATCAAGTCGGTGGCCCTTTTGGGTTGAGCATCCTCAGCGCGCGTATTGTCCATTGATCAGTGACCTGATCGAGGCCCGCCGCCCATGCTTGACTGCCCCCTGCCGCGCACCCTGCACTACGTTGACGACAGTCAGCCGGGCCTGACCCGTCGGCGCTGGCGCGACCGTTTCATCTACCTGGATGCCGATGGCCAACGGG  >  NC_002947/4356455‑4356746

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