Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,087,541 C→A 8.1% intergenic (+190/‑236) PP_3597 → / → PP_3598 amino acid ABC transporter ATP‑binding protein/peptidase C26

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,087,5410CA 8.1% 151.4 / 2.8 49intergenic (+190/‑236)PP_3597/PP_3598amino acid ABC transporter ATP‑binding protein/peptidase C26
Reads supporting (aligned to +/- strand):  ref base C (21/24);  new base A (2/2);  total (23/26)
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.96e-01

CGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGA  >  NC_002947/4087399‑4087684
                                                                                                                                              |                                                                                                                                               
cgTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCAcgc                                                                                                                                           <  8:207127/149‑1 (MQ=255)
    gcgaacACGGGCGAAGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCAcgccgc                                                                                                                                        >  1:305271/6‑148 (MQ=255)
          acaggCTTCGCCCGTGTTCGCGGGCAAGCCCACACCCAAAGTCACCGACGTAACTACAAGTTACCTAAACAATGCGCCGCTGCAATAACCCAACTGCATCGATAAATATTATTGACCTCATACACTCCCAAACCCACGCCGCCACACGc                                                                                                                                 >  6:201391/4‑149 (MQ=255)
           cGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGcc                                                                                                                                <  5:82392/149‑1 (MQ=255)
              gCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCg                                                                                                                              >  7:268228/1‑148 (MQ=255)
                    cccGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGt                                                                                                                       <  3:230859/149‑1 (MQ=255)
                     ccGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATAATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGtt                                                                                                                      >  6:21284/1‑149 (MQ=255)
                      cGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGttt                                                                                                                     <  2:159929/149‑1 (MQ=255)
                       gtgtTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGTCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTc                                                                                                                    <  3:8059/149‑1 (MQ=255)
                               gggCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGcccc                                                                                                              >  1:326389/1‑147 (MQ=255)
                                        cACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTc                                                                                                                    >  3:3770/1‑132 (MQ=255)
                                        cACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATAATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTc                                                                                                                    <  4:3770/132‑1 (MQ=255)
                                          cTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATGATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGt                                                                                                 <  5:21284/149‑1 (MQ=255)
                                            cccACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCACGCATCGGTTTCAGCCCCCAGAAAACGGCGTAt                                                                                               >  1:237603/1‑149 (MQ=255)
                                            cccACAGTGACCGAGGTTACTTCAAGTTACCTAATAATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTAt                                                                                               >  8:108596/1‑149 (MQ=255)
                                             ccACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATc                                                                                              <  6:92851/149‑1 (MQ=255)
                                                     aCCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATcc                                                                                             >  3:20532/1‑142 (MQ=255)
                                                     aCCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATcc                                                                                             <  4:20532/142‑1 (MQ=255)
                                                     aCCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCcgcg                                                                                       <  2:305271/148‑1 (MQ=255)
                                                     aCCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCcgcg                                                                                       <  2:281107/148‑1 (MQ=255)
                                                        gAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTaa                                                                                   <  2:326389/149‑1 (MQ=255)
                                                                         ccTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCATACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTaaa                                                                                  <  1:305446/133‑1 (MQ=255)
                                                                         ccTAATCATTGCGCCGCTGCAATGACCCAACAGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCAAGAAAACTGCGTATCACCCGCGTTaaa                                                                                  >  2:305446/1‑133 (MQ=255)
                                                                          cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGttt                                                                 >  3:196875/1‑149 (MQ=255)
                                                                           tAATCATTGCGCCGCTGCAATGACACAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGttt                                                                 >  6:315564/1‑148 (MQ=255)
                                                                            aaTCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGa                                                               >  1:213730/1‑149 (MQ=255)
                                                                                   gcgcCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGc                                                                                                   <  8:322958/106‑1 (MQ=255)
                                                                                     gccgcTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGaa                                                      >  4:81300/1‑149 (MQ=255)
                                                                                         tctcAATGACCCAACGTCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACg                                                    <  7:149759/144‑1 (MQ=255)
                                                                                            cAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAAACCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGctct                                               >  3:11373/1‑149 (MQ=255)
                                                                                            cAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAAACCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGctct                                               >  3:16039/1‑149 (MQ=255)
                                                                                             aaTGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCTAACg                                                    >  8:149759/1‑143 (MQ=255)
                                                                                                  cccAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCaa                                          <  7:108596/148‑1 (MQ=255)
                                                                                                           aTCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTGAGCCCCCAGAAAACGGCGTATCCCCCGCGTTaaa                                                                                  >  8:149106/1‑99 (MQ=255)
                                                                                                           aTCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTGAGCCCCCAGAAAACGGCGTATCCCCCGCGTTaaa                                                                                  <  7:149106/99‑1 (MQ=255)
                                                                                                            tCGTTAAATATTATTGACCGCAGACGCTGCCATACCCACTCCGCCACAGGCCCGCATAGGTTTCATCCCCCAGAAAAAGGCGTATCCACCGCGTTAATGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACaa                               <  5:201391/149‑1 (MQ=255)
                                                                                                               ttAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGACCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAAAGGCTCTGGCAAGCTTTCTACACGcc                            >  5:180999/1‑149 (MQ=255)
                                                                                                                  aaTATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCGCAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCg                         >  3:19281/1‑149 (MQ=255)
                                                                                                                    tattatTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGAc                       >  5:220975/1‑149 (MQ=255)
                                                                                                                      ttattGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCg                     <  6:220975/149‑1 (MQ=255)
                                                                                                                        attGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGa                   <  2:213730/149‑1 (MQ=255)
                                                                                                                              cGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAAtcgtc             <  3:81300/149‑1 (MQ=255)
                                                                                                                              cGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAAtcgtc             >  8:181954/1‑149 (MQ=255)
                                                                                                                              cGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAAtcgtc             <  8:268228/149‑1 (MQ=255)
                                                                                                                              agCAGACCCTGCCAAACCCACGCCGCCAGAGGCCCGCATCGGTTTCAGCCCACAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAAtcgtc             <  6:180999/148‑1 (MQ=255)
                                                                                                                               gCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTg                                              >  7:86303/1‑115 (MQ=255)
                                                                                                                               gCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTg                                              <  8:86303/115‑1 (MQ=255)
                                                                                                                                         ccAAAACCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGa  <  4:11373/149‑1 (MQ=255)
                                                                                                                                         ccAAAACCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGa  <  4:16039/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
CGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACATATCGAACGGCTCTGGCAAGCTATCTACAAGCCCCGACCGGAATCGTCGCCTCTGACGA  >  NC_002947/4087399‑4087684

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: ‑