Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,102,237 T→A 11.5% L9M (TTG→ATG) ‡ PP_2712 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,102,2370TA11.5% 117.8 / 7.0 44L9M (TTG→ATG) ‡PP_2712hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/20);  new base A (2/3);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AAAAGCTCTTCGTGCGTATCAACGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGCCTTCGAAAAGCTGGCAGG  >  NC_002947/3102098‑3102376
                                                                                                                                           |                                                                                                                                           
aaaaGCTCTTCGTGCGTATCAACGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCAAGCTGGCaat                                                                                                                                   <  2:718164/150‑3 (MQ=255)
aaaaGCTCTTCGTGCGTATCAACGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCAAGCTGGCaat                                                                                                                                   <  2:168660/150‑3 (MQ=255)
 aaaGCTCTTCGTGCGTATCAACGGCCTCTGACCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCAttt                                                                                                                                  >  1:451148/1‑150 (MQ=255)
    gCTCTTCGTGCGTATCAAAGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGACGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCa                                                                                                                               <  1:683984/150‑1 (MQ=255)
     ctctTCGTGCGTATCAACGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGTAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCAc                                                                                                                              >  2:574582/1‑150 (MQ=255)
      tctTCGTGCGTATCAACGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCAcc                                                                                                                             <  1:375799/150‑1 (MQ=255)
                  tCAACGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGCCATCTGGCCCCTTGACCTGGGc                                                                                                                 >  2:370841/1‑150 (MQ=255)
                    aaCGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCAt                                                                                                               >  2:656317/1‑150 (MQ=255)
                          aTGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCAATGCCAGCAAGCTGGCaat                                                                                                                                   <  1:672856/124‑3 (MQ=255)
                          aTGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCAATGCCAGCAAGCTGGCaat                                                                                                                                   >  2:672856/1‑122 (MQ=255)
                                   ggCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCCCCGTGAGCTGGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAggcaggca                                                                                                <  1:370841/150‑1 (MQ=255)
                                   ggCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAggcaggca                                                                                                <  1:475627/150‑1 (MQ=255)
                                   ggCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAggcaggca                                                                                                <  1:796124/150‑1 (MQ=255)
                                   ggCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAggcaggca                                                                                                <  1:796233/150‑1 (MQ=255)
                                     cATGGTCGATCGCGCCCTACGCAAGAAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGcagc                                                                                              <  2:273306/150‑1 (MQ=255)
                                      aTGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCAtt                                                                                                              >  1:259896/1‑133 (MQ=255)
                                      aTGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCAtt                                                                                                              <  2:259896/133‑1 (MQ=255)
                                               cgcgCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGTAGTCGAACAATATAACCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGa                                                                                    <  1:574582/150‑1 (MQ=255)
                                                gcgcTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACcaggca                                                                                                    >  1:329266/1‑133 (MQ=255)
                                                gcgcTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACcaggca                                                                                                    <  2:329266/133‑1 (MQ=255)
                                                  gcTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTg                                                                                 >  1:815243/1‑150 (MQ=255)
                                                    tctACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCt                                                                                                                                         <  1:205807/92‑1 (MQ=255)
                                                    tctACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCt                                                                                                                                         >  2:205807/1‑92 (MQ=255)
                                                     ctACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCACTGAGGCAGCCAGCGCCTGAGTGATc                                                                              <  1:331418/150‑1 (MQ=255)
                                                                    ccAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTggg                                                               <  2:83122/150‑1 (MQ=255)
                                                                      aGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGgcg                                                             >  2:837329/1‑150 (MQ=255)
                                                                       gTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTg                                                                 >  1:467341/1‑145 (MQ=255)
                                                                       gTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTg                                                                 <  2:467341/145‑1 (MQ=255)
                                                                                 gcTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTa                                                  >  2:24493/1‑150 (MQ=255)
                                                                                     tatCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCaa                                              >  2:507784/1‑150 (MQ=255)
                                                                                      atCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTACCTGTAAACGGTCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCaaa                                             <  2:451148/150‑1 (MQ=255)
                                                                                         ggCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGc                                          >  2:483335/1‑150 (MQ=255)
                                                                                             gCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCaa                                      <  2:815243/150‑1 (MQ=255)
                                                                                                  tCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATa                                 >  1:685214/1‑150 (MQ=255)
                                                                                                               aCAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGc                    <  1:196705/150‑1 (MQ=255)
                                                                                                                     gAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGCCTTCGa              >  2:289236/1‑150 (MQ=255)
                                                                                                                       aGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGCCTTCGaaa            <  2:678889/150‑1 (MQ=255)
                                                                                                                        gCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGCCTTCGaaaa           <  1:24493/150‑1 (MQ=255)
                                                                                                                          gTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGa                                                                                    >  1:505687/1‑75 (MQ=255)
                                                                                                                          gTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGa                                                                                    <  2:505687/75‑1 (MQ=255)
                                                                                                                              gTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGCCTTCGAAAAGCTGGc     >  2:200809/1‑150 (MQ=255)
                                                                                                                               ttCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTg                                                                                                             <  1:523869/45‑1 (MQ=255)
                                                                                                                               ttCATTGCCAGCTTGCTGGCATTTGCACCCTTCACATGGGCATTg                                                                                                             >  2:523869/1‑45 (MQ=255)
                                                                                                                                 aaatGCCAGCAAGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGCCTTCGAAAAGCTGGCAgg  >  2:668778/4‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                           
AAAAGCTCTTCGTGCGTATCAACGGCATGTTGCCGGGCATGGTCGATCGCGCTCTACGCAAGCAACTGCCAGTCATTCGCCGCTATATCGGCAGCCACTCCAAGGAGTCGAACAAATGAAGCGTCTGTTCATTGCCAGCTTGCTGGCATTTGCACCCTTCACCTGGGCATTGGACCAGGCAGGCAGCCAGCGCCTGAGTGATCTGCAGCAGCGCTGGGCGCAGATTCAGTACCAAATGCCCAAGGATAAGCGTGCCGACGCCTTCGAAAAGCTGGCAGG  >  NC_002947/3102098‑3102376

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