Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,756 C→A 14.2% intergenic (+65/+62) 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,7560CA14.2% 65.8 / 2.3 27intergenic (+65/+62)PP_2445/PP_2446membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (10/13);  new base A (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 = 1.00e+00

TCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGC  >  NC_002947/2793620‑2793902
                                                                                                                                        |                                                                                                                                                   
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)
                                                                                                                                        |                                                                                                                                                   
TCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTGGTGTGCATGACGC  >  NC_002947/2793620‑2793902

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: