Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,793,752 G→T 17.9% intergenic (+61/+66) 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,7520GT17.9% 55.7 / 7.2 28intergenic (+61/+66)PP_2445/PP_2446membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (9/14);  new base T (2/3);  total (11/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

CGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTG  >  NC_002947/2793606‑2793889
                                                                                                                                                  |                                                                                                                                         
cGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTgg                                                                                                                                          <  7:156159/148‑1 (MQ=255)
 ggCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCa                                                                                                                                        <  1:295532/149‑1 (MQ=255)
     aTGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTTGccagtg                                                                                                                                    <  7:278224/149‑6 (MQ=255)
              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)
                                              ccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTg                                                                                                                                           <  5:320641/101‑1 (MQ=255)
                                              ccGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTg                                                                                                                                           >  6:320641/1‑101 (MQ=255)
                                                                  gTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCAGCCCAAGTATTAGCGTGGCACAGGCAAAGCGAATAACAGTGCCCCGGGCGTGTCAGGCACTGGGCGTAGGGAAATCCTGAGGTTa                                                                       >  2:318813/1‑149 (MQ=255)
                                                                                          aaaaGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTAATCGCGAGGCCCTGTCACCGCGAATACCtgtac                                                                                                                  >  7:207725/1‑80 (MQ=37)
                                                                                                            cccTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTgg                                                                                                                                          <  2:184549/40‑1 (MQ=255)
                                                                                                            cccTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTgg                                                                                                                                          >  1:184549/1‑40 (MQ=255)
                                                                                                                    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)
                                                                                                                                                  |                                                                                                                                         
CGGCCATGGCGTTCTCGACGGCTATCGAGGTCCTCAACATCCTGGCCCGTCGAGCGCGTTTGAAGCGTGAGGCGGCTGAAGGTTGATGTAAAAAGGGGCGCCTCGCGCCCCTCTTGCGGACAAGCCCGCTGCACAGGTATTCGCGTGGCACTGGCAACGCGAATACCTGTGCAGCGGGCTTGTCAGTCAGTGGGCTTCGGGAAATCCTGAGGTTAGCTCACCCCGGCATGCCGCCGGGGCCGTCTTACTGGCTGTGGCTTGGCGGGTGCCTGCGGTGGCAGGTG  >  NC_002947/2793606‑2793889

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: