Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,622,333 C→T 12.1% intergenic (+124/‑224) PP_5515 → / → PP_5516 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,622,3330CT12.1% 82.4 / 2.9 33intergenic (+124/‑224)PP_5515/PP_5516hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base T (2/2);  total (16/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 = 1.00e+00

AAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTCACAGGCGGCCATCTTCGGCCCTTTGTGGGGCTTTTTTGG  >  NC_002947/2622189‑2622476
                                                                                                                                                |                                                                                                                                               
aaGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACTCTCCGGGCAGGCCCCAAGTGGCCCAGCGGc                                                                                                                                              <  8:58494/148‑1 (MQ=255)
    gCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGcc                                                                                                                                          >  7:325716/1‑148 (MQ=255)
         tGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCCGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAAtgt                                                                                                                                    <  8:102671/149‑1 (MQ=255)
                              aaaaaTTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGTGGCAGCCCAATGTGTCCCGGAATGTGTGCAACAt                                                                                                               >  6:105236/1‑149 (MQ=255)
                              aaaaaTTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGTGGCAGCCCAATGTGTCCCGGAATGTGTCCCACAt                                                                                                               >  6:110858/1‑149 (MQ=255)
                              aaaaaTTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACAt                                                                                                               >  8:30107/1‑149 (MQ=255)
                                   tttGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGGGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTg                                                                                                          <  1:323323/149‑1 (MQ=255)
                                       agagCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGaaa                                                                                                      <  5:372543/149‑1 (MQ=255)
                                                               ccGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCa                                                                              >  4:299915/1‑149 (MQ=255)
                                                                 gcgACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCAc                                                                             >  1:10788/1‑148 (MQ=255)
                                                                 gcgACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCAc                                                                             >  5:213714/1‑148 (MQ=255)
                                                                   gACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAAccc                                                                                                   <  8:129991/124‑1 (MQ=255)
                                                                   gACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAAccc                                                                                                   >  7:129991/1‑124 (MQ=255)
                                                                      tACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGAtt                                                                       <  3:68306/149‑1 (MQ=255)
                                                                         ccccGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGAtt                                                                       >  1:371003/1‑146 (MQ=255)
                                                                           ccGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCTGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGaa                                                                  <  8:325716/149‑1 (MQ=255)
                                                                                  cTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCcgcg                                                           <  2:371003/149‑1 (MQ=255)
                                                                                     cAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCgcgc                                                          <  6:59494/147‑1 (MQ=255)
                                                                                             agcCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGAACCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCCCCAGATTAAGAATCCCGCGCTCCACCTTTc                                                >  7:87203/1‑149 (MQ=255)
                                                                                               ccAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGTGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCa                                              <  5:110858/149‑1 (MQ=255)
                                                                                               ccAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGAGTAAGAATCCCGCGCTCCACCTTTc                                                >  4:113169/1‑147 (MQ=255)
                                                                                               ccAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGACCCAAGTGGCCCAGTGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCa                                              <  5:105236/149‑1 (MQ=255)
                                                                                                cAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAg                                             >  1:174765/1‑149 (MQ=255)
                                                                                                cAAAATCCGCACTATGAGAACCATCCGGGCAGGCCCCAAGAGGACCCGCGGAAGCCCAAGCTGACCCGGAATGTGTCCCACATGCCTGGAAACCCCAGTAAACACGGGACCTGTAACCAGATTAAGAATCCCGCGCACAACCTTTCAAg                                             >  2:236926/1‑149 (MQ=255)
                                                                                                     tCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCTGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTca                                        <  4:93712/149‑1 (MQ=255)
                                                                                                       cGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTca                                        >  2:245502/1‑147 (MQ=255)
                                                                                                                               ggCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCATCCTCACATGCGGCCATCTTCGGCCCTTtttg              <  2:351132/149‑1 (MQ=255)
                                                                                                                                gCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTCACAGGCGGCCATCTTCGGCCCTTTGTgg             >  2:7997/1‑149 (MQ=255)
                                                                                                                                 ccccAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGACTGTGTCCCACATGGCTGGAACCCCCCGTAAACACGGGACCTGTCACCAGCTGAAGACTCCCGCGCGCCACCTGTCCAGCCTCACAGGCGGCCAGCTGCGGCCCTTGTgggg           >  5:78012/1‑149 (MQ=255)
                                                                                                                                 ccccAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTCACAGGCGGCCATCTTCGGCCCTTTGTggg            <  6:213714/149‑1 (MQ=255)
                                                                                                                                       gTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTCACAGGCGGCCATCTTCGGCCCTTTGTGGGGCtttt      <  3:37465/149‑1 (MQ=255)
                                                                                                                                           cccAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTCACAGGCGGCCATCTTCGGCCCTTTGTGGGGCTTTTTTgg  <  3:299915/149‑1 (MQ=255)
                                                                                                                                           cccAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTCACAGGCGGCCATCTTCGGCCCTTTGTGGGGCTTTTTTgg  <  2:10788/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AAGGGCAAATGGCCCAAATGATTTTGTTGCAAAAATTTGAGAGCCCACCTCGAAATTGATAGGCCGCGACTACCCCCGGTGCCTCCAGCGAGCAGCCAAAATCCGCCCTCTGAGAACCCTCCGGGCAGGCCCCAAGTGGCCCAGCGGCAGCCCAATGTGTCCCGGAATGTGTCCCACATGGCTGGAAACCCCAGTAAACACGGGACCTGTCACCAGATTAAGAATCCCGCGCTCCACCTTTCCAGCCTCACAGGCGGCCATCTTCGGCCCTTTGTGGGGCTTTTTTGG  >  NC_002947/2622189‑2622476

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