Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 336,288 T→A 19.7% intergenic (+136/+257) PP_0278 → / ← PP_0279 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947336,2880TA19.7% 45.5 / 3.2 24intergenic (+136/+257)PP_0278/PP_0279hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/6);  new base A (3/2);  total (16/8)
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

CCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCTGGTCTTTGCGCAAACGCCTGT  >  NC_002947/336141‑336432
                                                                                                                                                   |                                                                                                                                                
acAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGCCCGCtt                                                                                                                                                 >  3:6006/2‑149 (MQ=255)
      cgcTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCATTCGTAGCGCAAGGGCGGCAAAGCAGCCCCCCCCATTTAAAAGAATGCGCTGAAAAACGGGGCCCCATTGCGGc                                                                                                                                           >  2:42826/1‑149 (MQ=255)
         tgatgaACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCAAAGCGGcccc                                                                                                                                        <  6:101896/149‑2 (MQ=255)
              aaCGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTcg                                                                                                                                    <  3:119489/148‑1 (MQ=255)
                cGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGGTGCAAATGAGCCCCGGCGCTTTGAGAGATTGCGCGGAAAACGGGGGCCGCTTTTCGGGCCGTTtgcgg                                                                                                                                 >  1:50574/1‑148 (MQ=255)
                                          cgcATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCATTCGTAGCGAAAGGGCTGCAAAACACCCCCCGCGAATTTAAAGAATTCGCAGAAAAACGGGGCCCCATTTAGGCCCATTAGCCAACCAAAGCCCCACCCACaaaaaacc                                                                                                        >  2:34962/1‑144 (MQ=255)
                                                                ttttGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTg                                                                                                                                               <  8:174743/86‑1 (MQ=255)
                                                                ttttGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTg                                                                                                                                               >  7:174743/2‑87 (MQ=255)
                                                                                      cgcgAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTgg                                                            >  2:212094/1‑148 (MQ=255)
                                                                                                 tGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGa                                                 >  5:45650/1‑148 (MQ=255)
                                                                                                  gCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGa                                                 >  2:254505/1‑147 (MQ=255)
                                                                                                       gcagcCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAggag                                           >  1:24754/1‑146 (MQ=255)
                                                                                                       gcagcCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAg                                              <  2:24754/145‑1 (MQ=255)
                                                                                                         agcCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGc                                        <  1:254505/149‑1 (MQ=255)
                                                                                                           ccccGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCg                                      >  8:27323/1‑149 (MQ=255)
                                                                                                              cGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGttt                                                                               <  8:131108/105‑1 (MQ=255)
                                                                                                                          agaTTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCt                       >  5:230466/1‑149 (MQ=255)
                                                                                                                          agaTTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCt                       >  4:297656/1‑149 (MQ=255)
                                                                                                                             ttGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCTGGt                    <  7:27323/149‑1 (MQ=255)
                                                                                                                                  cTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCTGGTCTTTg               >  6:280663/1‑149 (MQ=255)
                                                                                                                                   tGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCTGGTCTTTgc              >  2:69247/1‑149 (MQ=255)
                                                                                                                                     aaaaCGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCTGGTCTTTgcg             >  7:184792/1‑148 (MQ=255)
                                                                                                                                        gaaagggCCGCAAAGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTAcg                            <  5:160967/126‑1 (MQ=39)
                                                                                                                                        gaaagggCCGCAAAGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTAcg                            >  6:160967/5‑130 (MQ=39)
                                                                                                                                               ccGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCTGGTCTTTGCGCAAACGCCTGt  >  8:42231/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTGCGCTGAAAACGGGGGCCGCTTTGCGGCCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCCGCTTGCGCTAACCACAAGAAATGTTTATGCACAAGTGGCATTTGGGACCTGAGCGAAAGGCATGCCGGATCGCTACGCCGCTGGTCTTTGCGCAAACGCCTGT  >  NC_002947/336141‑336432

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: