Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,042,395 A→T 13.4% E92D (GAA→GAT) ‡ PP_5295 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,042,3950AT13.4% 72.4 / 5.0 30E92D (GAA→GAT) ‡PP_5295hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/10);  new base T (2/2);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GCCGGCAGCCAAGGCACCCTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCAT  >  NC_002947/6042246‑6042540
                                                                                                                                                     |                                                                                                                                                 
gcCGGCAGCCAAGGCACCCTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGaa                                                                                                                                                   <  1:446718/150‑1 (MQ=255)
    gCAGCCAAGGCACCCTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCATGCGCAAGGGCCTGTCACCCAGCAAGGTGGAATGCa                                                                                                                                               >  2:78781/1‑150 (MQ=255)
           aGGCACCCTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTgcg                                                                                                                                        >  1:513257/1‑150 (MQ=255)
             gCACCCTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGcc                                                                                                                                      >  2:22290/1‑150 (MQ=255)
                 ccTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTAGGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCaa                                                                                                                                  >  2:141302/1‑150 (MQ=255)
                 ccTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCaa                                                                                                                                  <  2:175726/150‑1 (MQ=255)
                               tGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTACGCGTGCCCGAAGCTCTGCGCGAACTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGGACCTTGCGCCTCAACGAAGACAGCAATg                                                                                                                    >  1:569229/1‑150 (MQ=255)
                                       gTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGCATTCCACCTTgc                                                                                                                                         <  1:252931/121‑1 (MQ=255)
                                       gTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGCATTCCACCTTgc                                                                                                                                         >  2:252931/1‑121 (MQ=255)
                                                       tGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATg                                                                                            >  2:74809/1‑150 (MQ=255)
                                                            ccGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGAc                                                                                                >  2:272811/1‑141 (MQ=255)
                                                            ccGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGAc                                                                                                <  1:272811/141‑1 (MQ=255)
                                                              gCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGCATTCCACCTTg                                                                                                                                          <  1:16705/97‑1 (MQ=255)
                                                              gCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGCATTCCACCTTg                                                                                                                                          >  2:16705/1‑97 (MQ=255)
                                                                 ccTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGcgcgccgc                                                                                  <  1:78781/150‑1 (MQ=255)
                                                                        cTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTg                                                                           >  2:345509/1‑150 (MQ=255)
                                                                         tGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTgg                                                                          >  1:513468/1‑150 (MQ=255)
                                                                          gCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTgcgc                                                                                                                                       >  1:423367/1‑88 (MQ=255)
                                                                          gCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTgcgc                                                                                                                                       <  2:423367/88‑1 (MQ=255)
                                                                                     tGCGCGACATCGAAGGCTGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGccgccgccg                                                              >  2:596699/1‑150 (MQ=255)
                                                                                         cgACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGa                                                          >  2:547809/1‑150 (MQ=255)
                                                                                                    gCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGcc                                               >  1:608378/1‑150 (MQ=255)
                                                                                                       gtgtACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGa                                            <  2:608378/150‑1 (MQ=255)
                                                                                                                aaGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTgg                                                                                                  <  1:263254/87‑1 (MQ=255)
                                                                                                                aaGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTgg                                                                                                  >  2:263254/1‑87 (MQ=255)
                                                                                                                  ggCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCg                                 >  2:581137/1‑150 (MQ=255)
                                                                                                                     cTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGcc                              >  1:545864/1‑150 (MQ=255)
                                                                                                                        cgcCAGGGCCTGTCCCGCGCCAAGGTGGAATGCACCTTGCGCCTCCACGAAGACAGCAATGGCAAACCCCTGAAAGCGGACCATGAGCGCGCCTCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGccccc                           <  1:188718/150‑1 (MQ=255)
                                                                                                                                        gcgGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAg           <  1:101480/150‑1 (MQ=255)
                                                                                                                                          ggCAAGGTGGCATGCACCTTGCGCCTCAACGAATACAGCAATGGCAAACCGCTGAAAGTGGACCCTGTATCCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTg         <  1:596699/150‑1 (MQ=255)
                                                                                                                                                 tGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCAt  <  1:404707/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCCGGCAGCCAAGGCACCCTGGTCTGGGAGCTGCGTTCGGTCAACCACCGTTATCTGGAGCCGCACCTGCGCCTGCCCGAAGCCCTGCGCGACCTCGAAGGCGGTGTACGTGAAGGCCTGCGCCAGGGCCTGTCACGCGGCAAGGTGGAATGCACCTTGCGCCTCAACGAAGACAGCAATGGCAAACCGCTGAAAGTGGACCATGAGCGCGCCGCGCAACTGGTTGCCGCCGCCGAGGAAGTGGCTGGCCTGATCAAGCAGCCGGCCCCCCTGAACCCGCTGGAAGTGCTGTCAT  >  NC_002947/6042246‑6042540

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: