Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,761 A→C 14.5% intergenic (+70/+57) 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,7610AC14.5% 59.3 / 3.0 27intergenic (+70/+57)PP_2445/PP_2446membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base C (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 = 9.99e-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)
                               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: