Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,087,485 C→G 13.3% intergenic (+134/‑292) 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,4850CG13.3% 73.3 / 5.7 30intergenic (+134/‑292)PP_3597/PP_3598amino acid ABC transporter ATP‑binding protein/peptidase C26
Reads supporting (aligned to +/- strand):  ref base C (12/14);  new base G (2/2);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

ACTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACA  >  NC_002947/4087347‑4087625
                                                                                                                                          |                                                                                                                                              
aCTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATg                                                                                                                                      >  4:86647/1‑149 (MQ=255)
         aaGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTg                                                                                                                            <  1:174246/149‑1 (MQ=255)
          aGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGc                                                                                                                           <  4:209112/149‑1 (MQ=255)
             cTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCAGCGGCGCAATGat                                                                                                                                    <  5:83982/138‑2 (MQ=255)
             cTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCAGCGGCGCAATGat                                                                                                                                    >  6:83982/1‑137 (MQ=255)
                             tttACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGAc                                                                                                        >  1:220771/1‑149 (MQ=255)
                                aCCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCg                                                                                                      <  2:236677/148‑1 (MQ=255)
                                         gaagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGc                                                                                            <  4:85508/149‑1 (MQ=255)
                                          aagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCAGCGGCGCAATGatt                                                                                                                                   <  7:15196/110‑3 (MQ=25)
                                          aagaCACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCAGCGGCGCAATGatt                                                                                                                                   >  8:15196/1‑108 (MQ=25)
                                                ccgccgTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAAcc                                                                                      >  2:77128/1‑148 (MQ=255)
                                                 cgccgTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATg                                                                                                                                      >  6:8658/1‑100 (MQ=255)
                                                 cgccgTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATg                                                                                                                                      <  5:8658/100‑1 (MQ=255)
                                                      tGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCAcgccg                                                                               <  3:39960/149‑1 (MQ=255)
                                                         aTGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAAc                                                                                       <  2:223822/138‑1 (MQ=255)
                                                         aTGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAAc                                                                                       >  1:223822/1‑138 (MQ=255)
                                                            gCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAg                                                                         <  3:86647/149‑1 (MQ=255)
                                                                gggCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGccc                                                                     >  1:67875/1‑149 (MQ=255)
                                                                         ccGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGtt                                                           <  2:220771/149‑1 (MQ=255)
                                                                           gtgtTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCTCGCATCGGTTTc                                                         >  7:225114/1‑149 (MQ=255)
                                                                            tgtTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCa                                                        <  2:67875/149‑1 (MQ=255)
                                                                            tgtTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCa                                                        <  1:235395/149‑1 (MQ=255)
                                                                                         gCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACCTGCATCGTTAAATATTATTGACCGCAGACCCTGc                                                                                            <  1:36528/102‑1 (MQ=255)
                                                                                         gCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGc                                                                                            >  2:36528/1‑101 (MQ=255)
                                                                                                 ccACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATc                                   >  7:20140/1‑149 (MQ=255)
                                                                                                                      cAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCtt               <  8:20140/148‑1 (MQ=255)
                                                                                                                              cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTTACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCtttt      >  2:74832/1‑149 (MQ=255)
                                                                                                                              cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCACACGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGtt       <  1:77128/149‑1 (MQ=255)
                                                                                                                              cTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGttt      >  5:142315/1‑149 (MQ=255)
                                                                                                                                  tCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACa  >  5:32567/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                              
ACTGATTGAAAGCCTGTCACTGTGGGGGGTTTACCCGCGAAGAAGACACCGCCGTGGATGGCACGGGCTTCGCCCGTGTTCGCGGGCAAGCCCACTCCCACAGTGACCGAGGTTACTTCAAGTTACCTAATCATTGCGCCGCTGCAATGACCCAACTGCATCGTTAAATATTATTGACCGCAGACCCTGCCAAACCCACGCCGCCACAGGCCCGCATCGGTTTCAGCCCCCAGAAAACGGCGTATCCCCCGCGTTAAAGCCCGCTTCCGTCGTTTGACA  >  NC_002947/4087347‑4087625

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