Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,159,143 C→T 14.2% intergenic (+91/‑118) PP_3661 → / → PP_3662 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,159,1430CT14.2% 98.7 / 7.5 43intergenic (+91/‑118)PP_3661/PP_3662membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/15);  new base T (3/3);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

GGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAA  >  NC_002947/4158997‑4159283
                                                                                                                                                  |                                                                                                                                             
ggCGGGTTGGCGGGGGGCTGGCGATTCCCCTGGGGGGGCAGAAGCTGTTGGGCTGACTGGAGCCCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTTAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCAc                                                                                                                                             <  6:257749/149‑1 (MQ=39)
    ggTTGGCGGGTGGCTGGCGAATACCCTGGGCGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCcgc                                                                                                                                         <  8:326989/149‑1 (MQ=255)
                            ccTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTAGCGGGTAAACCCCCACCTACAAGTGGGTGCTGGCCTGAAGACGGCAACGCCCCAGGAGGAGCGGGTTAACCCGCGAAGCAGGCAGCACGGTATCaacg                                                                                                                 >  6:368057/1‑147 (MQ=255)
                                   gTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTca                                                                                                          <  7:348748/149‑1 (MQ=255)
                                      cAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAgcaggcag                                                                                                                              <  1:41720/126‑1 (MQ=255)
                                                  ggCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTg                                                                                           <  8:76477/149‑1 (MQ=255)
                                                                       gctgctTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtg                                                                     >  6:300700/1‑149 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc                                                                           <  6:73580/138‑1 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTtctc                                                                           >  5:73580/1‑138 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATggtggt                                                                >  3:308235/1‑149 (MQ=255)
                                                                            ttCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACAGTTCTCGGTATggtggt                                                                >  7:173226/1‑149 (MQ=255)
                                                                               gcgGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGgtgt                                                              <  1:24590/148‑1 (MQ=255)
                                                                                 gggTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCa                                                            >  1:67446/1‑148 (MQ=255)
                                                                                  gctAAACCCGCTCCTACACGTGGGTGATGCCCGGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTAAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACAATTCTCGGTAAGGAGGTGTCa                                                            >  2:56897/3‑147 (MQ=255)
                                                                                   gTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCa                                                            >  6:270981/1‑146 (MQ=255)
                                                                                     aaaCCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaag                                                         >  2:284414/1‑147 (MQ=255)
                                                                                      AACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaa                                                       >  5:270981‑M1/1‑148 (MQ=255)
                                                                                        cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGa                                                    >  5:287231‑M2/59‑149 (MQ=255)
                                                                                        cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCCAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaag                                                      >  6:369784‑M2/59‑147 (MQ=255)
                                                                                        cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGCCGACTCGGTCTCCAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAaagaag                                                      <  5:369784‑M2/89‑1 (MQ=255)
                                                                                        cccgctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGa                                                    >  2:347975‑M2/59‑149 (MQ=255)
                                                                                         CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTc                                                                                                           >  6:98119‑M1/1‑94 (MQ=255)
                                                                                         CCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTc                                                                                                           <  5:98119‑M1/94‑1 (MQ=255)
                                                                                          CGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATg                                                  >  5:300700‑M1/1‑149 (MQ=255)
                                                                                           gctcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGt                                                 >  3:38019‑M2/56‑149 (MQ=255)
                                                                                            CTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGtt                                                >  4:108109‑M1/1‑149 (MQ=255)
                                                                                             tcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtggtg                                                                  >  5:238445‑M2/54‑130 (MQ=255)
                                                                                             tcctacaggggccttgccgtcttcagggcagcacccacctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTAtggtg                                                                  <  6:238445‑M2/77‑1 (MQ=255)
                                                                                               CTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGAAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGttt                                               >  3:368547‑M1/1‑147 (MQ=255)
                                                                                                     GTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGttttt                                                                                                            <  1:225369‑M1/81‑1 (MQ=255)
                                                                                                     GTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCGCGGTCTCAAGGttttt                                                                                                            >  2:225369‑M1/1‑81 (MQ=255)
                                                                                                      tgggtgCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGcac                                      <  5:368057/149‑1 (MQ=255)
                                                                                                                 cTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACtt                           <  8:173226/149‑1 (MQ=255)
                                                                                                                 cTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACt                            >  6:200480/1‑148 (MQ=255)
                                                                                                                    cagggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTa                          <  4:38019‑M2/117‑1 (MQ=255)
                                                                                                                      gggcagcacccacctgtaggagcgggttTACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCCCAGGAAAGACTTAGATg                      <  6:287231‑M2/121‑1 (MQ=255)
                                                                                                                      gACGGCAATGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGAt                       <  1:284414/149‑1 (MQ=255)
                                                                                                                                  cctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAt          >  5:19783‑M2/17‑149 (MQ=255)
                                                                                                                                  cctgtaggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCAt          >  8:31218‑M2/17‑149 (MQ=255)
                                                                                                                                      taggagcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGtc       >  2:376392‑M2/13‑148 (MQ=255)
                                                                                                                                          agcgggttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCaa  <  1:342132‑M2/141‑1 (MQ=255)
                                                                                                                                               gttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGa                              >  2:370471‑M2/4‑116 (MQ=255)
                                                                                                                                               gttCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGa                              <  1:370471‑M2/113‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GGCTGGTTGGCGGGGGGCTGGCGATTACCCTGGGGGTGCAGAAGCTGTTGGGCTGACTGGAGACCGAGTCGGCTGCTTCGCGGGTAAACCCGCTCCTACAGGTGGGTGCTGCCCTGAAGACGGCAAGGCCCCTGTAGGAGCGGGTTCACCCGCGAAGCAGGCAGCTCGGTCTCAAGGTTTTTCACAGCCCCCCTCGTTGCCGTCCACATTTCTCGGTATGGTGGTGTCAAAGAAGGATGTTTTACCCGCACAGGAAAGACTTAGATGCCTCAACGCCATGTCATCAA  >  NC_002947/4158997‑4159283

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