Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,653,045 C→A 13.8% intergenic (‑95/+66) PP_4961 ← / ← PP_4962 lipoprotein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,653,0450CA13.8% 80.5 / 3.7 29intergenic (‑95/+66)PP_4961/PP_4962lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (13/12);  new base A (2/2);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

GCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCG  >  NC_002947/5652896‑5653193
                                                                                                                                                     |                                                                                                                                                    
gcagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCcgc                                                                                                                                                      <  2:462092/150‑1 (MQ=255)
                  gtgtCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGc                                                                                                                                    <  2:220711/150‑1 (MQ=255)
                                                      ttcaGTTTCCTCGATTTGGGGCTGCTGTGCCGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGc                                                                                                <  1:164962/150‑1 (MQ=255)
                                                           tttCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTAtt                                                                                           >  2:476617/1‑150 (MQ=255)
                                                                tCGATTTGGGGCTGCTGTGCAGCCCTTCTCTGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAccc                                                                                       >  2:65703/1‑149 (MQ=255)
                                                                tCGATTTGGGGCTGCTGTGCAGCCCTTCTCTGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAccc                                                                                       <  1:65703/149‑1 (MQ=255)
                                                                    tttGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAgg                                                                                                    <  2:440414/132‑1 (MQ=255)
                                                                    tttGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAgg                                                                                                    >  1:440414/1‑132 (MQ=255)
                                                                                      cccTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt                                                                >  1:7912/1‑150 (MQ=255)
                                                                                                    gCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTggtggt                                                     <  1:294294/147‑1 (MQ=25)
                                                                                                     cccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc                                                 >  2:277613/1‑150 (MQ=255)
                                                                                                     cccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAAGACGGCACAGCGGCCCCTTCTGCAATACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc                                                 >  1:184204/1‑150 (MQ=255)
                                                                                                      ccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGCAGcc                                                                                                                                   >  1:462869/1‑67 (MQ=17)
                                                                                                      ccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGCAGcc                                                                                                                                   <  2:462869/67‑1 (MQ=18)
                                                                                                        gCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                                                                               <  1:363355/117‑1 (MQ=255)
                                                                                                        gCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                                                                               >  2:363355/1‑117 (MQ=255)
                                                                                                        gCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcgtc                                              <  1:277613/150‑1 (MQ=255)
                                                                                                          tCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGc                                                                                                <  1:405050/98‑1 (MQ=255)
                                                                                                          tCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGc                                                                                                >  2:405050/1‑98 (MQ=255)
                                                                                                           cccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTa                                           <  1:148572/150‑1 (MQ=255)
                                                                                                            ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTAc                                          <  2:7912/150‑1 (MQ=255)
                                                                                                              acaGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGCGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACtt                                        >  1:61378/1‑150 (MQ=255)
                                                                                                                  gTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGAGAAgg                                                                                                                                                 <  2:303282/41‑1 (MQ=255)
                                                                                                                  gTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGAGAAgg                                                                                                                                                 >  1:303282/1‑41 (MQ=255)
                                                                                                                     cAGGGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTAGTACTTTTCCaag                                 <  2:177567/150‑3 (MQ=255)
                                                                                                                            gTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCaa                                                                                                            >  2:19322/1‑68 (MQ=255)
                                                                                                                            gTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCaa                                                                                                            <  1:19322/68‑1 (MQ=255)
                                                                                                                                          gggCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTg            >  2:443727/1‑150 (MQ=255)
                                                                                                                                                  ccGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGAc    >  2:347226/1‑150 (MQ=255)
                                                                                                                                                  ccGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGAc    >  1:488422/1‑150 (MQ=255)
                                                                                                                                                    gcgAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCg  >  1:485278/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCCCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCCTTCCGGGCCCACGCGCTTGCCGATGACCG  >  NC_002947/5652896‑5653193

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: