Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,498,893 G→A 8.9% intergenic (‑40/+1160) PP_5563 ← / ← PP_5565 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,498,8930GA 8.9% 122.5 / 2.2 45intergenic (‑40/+1160)PP_5563/PP_5565membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base A (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

GGCGTACACCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGCTTCGACATCGGCAATGTCACCCGCGCAGCCTTG  >  NC_002947/3498746‑3499037
                                                                                                                                                   |                                                                                                                                                
ggCGTACACCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGc                                                                                                                                                 <  3:20003/149‑1 (MQ=255)
   gTACACCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATg                                                                                                                                              >  7:14297/1‑149 (MQ=255)
   gTACACCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATg                                                                                                                                              >  4:247740/1‑149 (MQ=255)
      cacCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATTGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCAt                                                                                                                                               <  6:12727/145‑1 (MQ=255)
      cacCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCAt                                                                                                                                               >  5:12727/1‑145 (MQ=255)
       acCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCt                                                                                                                                           >  1:325949/1‑148 (MQ=255)
       acCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACAGGCATGGCt                                                                                                                                           >  4:9218/1‑148 (MQ=255)
        ccGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTg                                                                                                                                          >  6:78015/1‑148 (MQ=255)
              gTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAacaaca                                                                                                                                   <  4:67781/149‑1 (MQ=255)
                  tCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGtt                                                                                                                               >  7:52997/1‑149 (MQ=255)
                  tCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGtt                                                                                                                               >  6:279545/1‑149 (MQ=255)
                    ccATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGtt                                                                                                                               >  8:286060/1‑147 (MQ=255)
                    ccATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGtt                                                                                                                               >  6:172304/1‑147 (MQ=255)
                    ccATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGtt                                                                                                                               <  5:52859/147‑1 (MQ=255)
                    ccATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTc                                                                                                                              <  6:328985/148‑1 (MQ=255)
                     cATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAg                                                                                                                            <  5:78015/149‑1 (MQ=255)
                         gCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACg                                                                                                                        <  5:172304/149‑1 (MQ=255)
                           gAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGt                                                                                                                       >  5:123314/1‑148 (MQ=255)
                            aGCATGGACATCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTa                                                                                                                     <  3:268674/149‑1 (MQ=255)
                            aGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGACATGGCTGAACAACAGGTTCAGCACGTTa                                                                                                                     >  8:194370/1‑149 (MQ=255)
                            aGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGACATGGCTGAACAACAGGTTCAGCACGTTa                                                                                                                     >  8:189428/1‑149 (MQ=255)
                               aTGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATca                                                                                                                  <  8:285328/149‑1 (MQ=255)
                                  gACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAg                                                                                                               >  7:236462/1‑149 (MQ=255)
                                           aGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCa                                                                                                      >  4:249665/1‑149 (MQ=255)
                                              ccGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCTGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTg                                                                                                   <  8:120499/149‑1 (MQ=255)
                                               cGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGACATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTgg                                                                                                  <  7:189428/149‑1 (MQ=255)
                                               cGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGACATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTgg                                                                                                  <  7:194370/149‑1 (MQ=255)
                                                             cGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCTCAGGACCAGGCACTGGACCTTTGCCGATAc                                                                                    <  2:325949/149‑1 (MQ=255)
                                                                 aaaGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATg                                                                                <  6:123314/149‑1 (MQ=255)
                                                                       tAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTgcgc                                                                          <  5:279545/149‑1 (MQ=255)
                                                                       tAGGGGACCTTACGTAGAGATCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGCTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTgcgc                                                                          >  6:83664/1‑149 (MQ=255)
                                                                         ggggACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCt                                                                         <  7:286060/148‑1 (MQ=255)
                                                                               cTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCa                                                                                                      >  2:227599/1‑113 (MQ=255)
                                                                                   cgcagagGTCGGGATGGTCCTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCCGCATGGCTGAACAACAGGTTCAGCACGTTATCACCGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGc                                                                <  8:52997/144‑1 (MQ=255)
                                                                                                          aTGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGcc                                        <  8:14297/148‑1 (MQ=255)
                                                                                                             tgtgCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGtt                                     <  3:9218/148‑1 (MQ=255)
                                                                                                              gtgCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGc                                   <  8:148794/149‑1 (MQ=255)
                                                                                                                 ccTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGCtt                                 <  8:236462/148‑1 (MQ=255)
                                                                                                                          tCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGACGTTGCTTCGAAATCgg                        >  4:190324/1‑148 (MQ=255)
                                                                                                                            ccTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACcgcg                                                 >  1:325548/1‑121 (MQ=255)
                                                                                                                            ccTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACcgcg                                                 <  2:325548/121‑1 (MQ=255)
                                                                                                                                             aaaCCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGCTTCGACATCGGCAATGTCACCCGCGCAGCCt    >  2:260381/1‑149 (MQ=255)
                                                                                                                                              aaCCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGCTTCGACATCGGCAATGTCACCCGCGCAGCCtt   <  4:133696/149‑1 (MQ=255)
                                                                                                                                               aCCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGCTTCGACATCGGCAATGTCACCCGCGCAGCCTTg  >  6:13742/1‑149 (MQ=255)
                                                                                                                                                ccGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGCTTCGACATCGGCAATGTCACCCGCGCAGCCTTg  >  4:113930/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGCGTACACCGCCAGTGCTCCCATGGCGAGCATGGACAGCATCAGACCGATGGCCGTAAGCCGGAAAAGAATAGGGGACCTTACGTAGAGTTCGGGATGGTCTTTCATGTGTGCCTGAGCATTCCCTGGCTGCCGCGGATCAAACCGGCATGGCTGAACAACAGGTTCAGCACGTTATCACAGGACCAGGCACTGGACCTTTGCCGATACCATGGTGCGCTGGTACCGGCAGGGTGAAAACCGCGTCAACGGCCGTTGCTTCGACATCGGCAATGTCACCCGCGCAGCCTTG  >  NC_002947/3498746‑3499037

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: ‑