Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,364,659 A→T 9.7% intergenic (+29/+120) dctA‑I → / ← PP_1189 C4‑dicarboxylate transport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,364,6590AT 9.7% 114.1 / 3.5 41intergenic (+29/+120)dctA‑I/PP_1189C4‑dicarboxylate transport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (17/20);  new base T (2/2);  total (19/22)
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.79e-01

TGGGTCAAGGAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCTCAGCCCAGTACCACTTCCAACACCCGG  >  NC_002947/1364514‑1364806
                                                                                                                                                 |                                                                                                                                                   
tGGGTCAAGGAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAgg                                                                                                                                                  <  5:184319/149‑1 (MQ=255)
  ggTCAAGGAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTc                                                                                                                                                <  3:56979/149‑1 (MQ=255)
    tCAAGGAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGc                                                                                                                                              <  8:251740/149‑1 (MQ=255)
    tCAAGGAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGc                                                                                                                                              >  7:251740/1‑149 (MQ=255)
      aaGGAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCAGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCt                                                                                                                                            >  5:261383/1‑149 (MQ=255)
         gAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCtttt                                                                                                                                         <  5:69238/149‑1 (MQ=255)
           aaTGGACAACGACAAGCTGGCTTCGGAACTGGGCTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTc                                                                                                                                       >  2:21145/1‑149 (MQ=255)
            aTGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCTTTGGTCGCCTTtctca                                                                                                                                      <  5:82851/149‑2 (MQ=255)
             tGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCgt                                                                                                                                     <  2:268472/149‑1 (MQ=255)
                       cAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGcgcg                                                                                                                            >  5:316139/1‑148 (MQ=255)
                           cTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGAccc                                                                                                                       <  2:78974/149‑1 (MQ=255)
                                         ggACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCTTTGGTCGCCttt                                                                                                                                          <  6:267280/116‑1 (MQ=255)
                                         ggACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCTTTGGTCGCCttt                                                                                                                                          >  5:267280/1‑116 (MQ=255)
                                          gACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGc                                                                                                        >  3:23338/1‑149 (MQ=255)
                                                       cgccgcTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTaaaa                                                                                           <  4:23338/149‑1 (MQ=255)
                                                           gcTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTg                                                                                       <  1:21145/149‑1 (MQ=255)
                                                                  gATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCtt                                                                                                          <  3:199721/123‑1 (MQ=255)
                                                                  gATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCtt                                                                                                          >  4:199721/1‑123 (MQ=255)
                                                                         gcccgACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCaaa                                                                                                                                                    <  8:314258/74‑1 (MQ=255)
                                                                                      ccTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCtt                                                            >  4:184724/1‑149 (MQ=255)
                                                                                       cTGGGCGTGGCTGAAGTCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCtt                                                            <  2:303734/148‑1 (MQ=255)
                                                                                       cTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTg                                                           >  5:285472/1‑149 (MQ=255)
                                                                                       cTGGGCGTGGCTGAAGCCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCtt                                                            >  1:303734/1‑148 (MQ=255)
                                                                                          ggCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGcc                                                         <  4:136641/148‑1 (MQ=255)
                                                                                             gTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAAt                                                                                                   >  1:135444/1‑103 (MQ=255)
                                                                                             gTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAAt                                                                                                   <  2:135444/103‑1 (MQ=255)
                                                                                             gTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCcgc                                                      >  2:220277/1‑148 (MQ=255)
                                                                                                    aaGGCCCGGCTCGTTGAGAGCGGGTTGAAATGAGAAAGGCGACCAACGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACgg                                              <  7:231978/149‑1 (MQ=255)
                                                                                                     aGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACggg                                             <  3:184724/149‑1 (MQ=255)
                                                                                                      ggCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGgcg                                                                                    <  5:344119/109‑1 (MQ=255)
                                                                                                      ggCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGgcg                                                                                    >  6:344119/1‑109 (MQ=255)
                                                                                                                gTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACggg                                             >  1:116061/1‑138 (MQ=255)
                                                                                                                gTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACggg                                             <  2:116061/138‑1 (MQ=255)
                                                                                                                gTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGt                                   <  6:316139/148‑1 (MQ=255)
                                                                                                                gTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTg                                  <  1:220277/149‑1 (MQ=255)
                                                                                                                gTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTg                                  <  6:285472/149‑1 (MQ=255)
                                                                                                                      agCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCt                             >  2:192827/1‑148 (MQ=255)
                                                                                                                           gTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCTCAGcc                        >  5:220104/1‑148 (MQ=255)
                                                                                                                           gTTGAAATGACACAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCTCAGcc                        <  4:201931/148‑1 (MQ=255)
                                                                                                                                 aTGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCTCAGCCCAGTAc                  >  8:89020/1‑148 (MQ=255)
                                                                                                                                  tGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTaaaa                                                                                           >  2:82/1‑74 (MQ=255)
                                                                                                                                     aaaaGGCGACCTTTGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCTCAGCCCAGTACCACtt             >  4:155043/2‑149 (MQ=255)
                                                                                                                                                aaaGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCTCAGCCCAGTACCACTTCCAACACCCgg  <  7:89020/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TGGGTCAAGGAAATGGACAACGACAAGCTGGCTTCGGAACTGGACTCCGGTGGTGCGCCGCTGATCGATACCCGCCCGACGGACGACCTGGGCGTGGCTGAAGGCCCGGCTCGTTGAGAGCTGGTTGAAATGAGAAAGGCGACCAAAGGTCGCCTTTTTCGTGTTAGCGCGGACCCTTGTGTGGGCCTTGCTCAATTTCTAAAAGCTGGCGCGATTGGACTGTGGGAGCGGGCTTGCCCGCGAACACGGGCGAAGCCCGTGGCCCTCAGCCCAGTACCACTTCCAACACCCGG  >  NC_002947/1364514‑1364806

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: