Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,758 G→T 14.9% intergenic (+67/+60) PP_2445 → / ← PP_2446 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,793,7580GT14.9% 57.7 / 4.9 27intergenic (+67/+60)PP_2445/PP_2446membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/13);  new base T (2/2);  total (12/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

TCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGCCGC  >  NC_002947/2793620‑2793905
                                                                                                                                          |                                                                                                                                                    
ccGACGGCTATCCAGCGCCTCAACCTCCTGGCCCGGCGAGCGCGTTCGCAGCGTGAGGCGGCTGAAGGATGATGTAAAAAGGCGCGCCGCGCGCCCCTCTTGCGGACAAGCCCGCTGCCCAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                                                                            <  8:295456/148‑1 (MQ=255)
 cGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGaa                                                                                                                                            >  8:206410/1‑148 (MQ=255)
  gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAAt                                                                                                                                           <  3:108018/148‑1 (MQ=255)
  gACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                                                                                          <  4:75099/149‑1 (MQ=255)
     ggCTATTGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGTCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCt                                                                                                                                       <  3:172399/149‑1 (MQ=255)
     ggCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATAcca                                                                                                                                       >  5:211336/1‑148 (MQ=38)
     ggCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCt                                                                                                                                       <  6:272494/149‑1 (MQ=38)
       cTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAACAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                                                                                     <  5:55577/149‑1 (MQ=255)
       cTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                                                                                     <  5:156223/149‑1 (MQ=255)
       cTATCGAGGTCCTCAACATCCTGGCCCATCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgt                                                                                                                                     <  8:227219/149‑1 (MQ=255)
          tCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCtgtg                                                                                                                                    >  7:303059/1‑147 (MQ=255)
                  cTCACCACCCTCGCCCGTCGTGCGCGTGTGTAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCtt                                                                                                                          <  4:236396/149‑1 (MQ=38)
                  cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAAt                                                                                                                                           <  1:128943/132‑1 (MQ=37)
                  cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                                                                                          >  2:301992/1‑133 (MQ=255)
                  cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATa                                                                                                                                          <  1:301992/133‑1 (MQ=255)
                  cTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCtt                                                                                                                          >  3:236396/1‑149 (MQ=255)
                     aaCATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAACGTTGATGTAAAAAGGGGCGCCTAGCGCCCCTCTTGAGGACAAGCCAGCTGCACAGGTATTCCCGTGGCACTg                                                                                                                                                      >  7:152397/1‑118 (MQ=255)
                               cccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAACCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTggg                                                                                                             <  7:206410/149‑1 (MQ=255)
                                                    gTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCAGCCCAAGTATTAGCGTGGCACAGGCAAAGCGAATAACAGTGCCCCGGGCGTGTCAGGCACTGGGCGTAGGGAAATCCTGAGGTTa                                                                                       >  2:318813/1‑149 (MQ=255)
                                                                            aaaaGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTAATCGCGAGGCCCTGTCACCGCGAATACCtgtac                                                                                                                                   >  7:207725/1‑80 (MQ=37)
                                                                                                      ctgACACGCCCGCTGCACAGGTATTCGCGGGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCgg                                     <  8:303059/147‑1 (MQ=255)
                                                                                                           aaGCCCGCTGCACAGGTATTCGCGTTGCCAGTGCCACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGcc                                >  3:218259/1‑149 (MQ=25)
                                                                                                                cGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTTGCTTGGCGGGTGCCTGCgg                           >  8:134286/1‑149 (MQ=255)
                                                                                                                        atgTAGTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAATTCCTGAGGTTAGCTCACAGCGGCATTCCGCGGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGt                                   <  5:198108/131‑1 (MQ=255)
                                                                                                                                gcgTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAggtg                  >  4:102536/1‑142 (MQ=255)
                                                                                                                                gcgTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAggtg                  <  3:102536/142‑1 (MQ=255)
                                                                                                                                       aCTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGGTTAGCTAACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGAcgc     >  1:256588/1‑149 (MQ=255)
                                                                                                                                         tGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGAcgccgc  >  4:71694/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
TCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGCCGC  >  NC_002947/2793620‑2793905

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: