Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,604,250 C→G 12.2% P12A (CCG→GCG) ‡ PP_4076 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,604,2500CG12.2% 118.2 / 8.7 41P12A (CCG→GCG) ‡PP_4076hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/16);  new base G (3/2);  total (23/18)
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.41e-01

ACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAT  >  NC_002947/4604108‑4604396
                                                                                                                                              |                                                                                                                                                  
aCGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTGGATGGCGGACGGGTCCGc                                                                                                                                             >  1:260369/1‑150 (MQ=255)
          tGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAg                                                                                                                                   <  1:307144/150‑1 (MQ=255)
                      gTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAg                                                                                                                       >  1:911432/1‑150 (MQ=255)
                            ggTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAgc                                                                                                                 >  1:623741/1‑150 (MQ=255)
                                   ctcgcCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGc                                                                                                          >  2:198894/1‑150 (MQ=255)
                                                     ttCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCa                                                                                        >  1:525312/1‑150 (MQ=255)
                                                       cAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTc                                                                                                                                                >  1:115253/1‑92 (MQ=255)
                                                       cAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTc                                                                                                                                                <  2:115253/92‑1 (MQ=255)
                                                        aaGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGt                                                                                     <  1:198894/150‑1 (MQ=255)
                                                              aTGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTGGATGGCGGACGGGTCCGCCATc                                                                                                                                         >  2:493931/1‑92 (MQ=255)
                                                              aTGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTGGATGGCGGACGGGTCCGCCATc                                                                                                                                         <  1:493931/92‑1 (MQ=255)
                                                                  gcgcGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGa                                                                           >  1:13817/1‑150 (MQ=255)
                                                                      gAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACg                                                                       >  1:317159/1‑150 (MQ=255)
                                                                          ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAgg                                                                                                                                  <  2:965754/87‑1 (MQ=255)
                                                                          ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAgg                                                                                                                                  >  1:965754/1‑87 (MQ=255)
                                                                          ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGc                                                                                                                                 >  2:426289/1‑88 (MQ=255)
                                                                          ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGc                                                                                                                                 <  1:426289/88‑1 (MQ=255)
                                                                          ttCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGc                                                                   >  2:708032/1‑150 (MQ=255)
                                                                                  ttCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTg                                                                 <  1:551820/144‑1 (MQ=255)
                                                                                  ttCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTg                                                                 >  2:551820/1‑144 (MQ=255)
                                                                                    cccTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGtt                                                                                                   >  2:748408/1‑108 (MQ=255)
                                                                                    cccTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGtt                                                                                                   <  1:748408/108‑1 (MQ=255)
                                                                                     ccTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtctc                                                                                                                            >  1:581380/1‑82 (MQ=255)
                                                                                     ccTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGtctc                                                                                                                            <  2:581380/82‑1 (MQ=255)
                                                                                         gTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCa                                                                                                                        >  1:622940/1‑82 (MQ=255)
                                                                                         gTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCa                                                                                                                        <  2:622940/82‑1 (MQ=255)
                                                                                         gTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAg                                                    >  2:737498/1‑150 (MQ=255)
                                                                                                      acacGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCagag                                       <  2:1178509/150‑1 (MQ=255)
                                                                                                        acGCAATGAGTGGTTTGTACGAACCCTTGATGGCGGACCCGTCCGCCCTCCAAGAGGCGCATCGCCAGGAAAGCGCTTGTCACAAGTTAGCTGCGATCAAGTGCGACCACGAAAAGATGCTTAACGCTCGCGAAGGTTGCTGCCAGagcc                                     >  2:588168/1‑149 (MQ=255)
                                                                                                                   ggTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTg                          <  2:623741/150‑1 (MQ=255)
                                                                                                                     tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGa                              <  2:372489/144‑1 (MQ=255)
                                                                                                                     tttCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGa                              >  1:372489/1‑144 (MQ=255)
                                                                                                                         ttcGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTg                    >  2:1015435/1‑150 (MQ=255)
                                                                                                                             aTCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGcgcg                >  1:372310/1‑150 (MQ=255)
                                                                                                                                  cggATGGCGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTc           >  1:851488/3‑150 (MQ=255)
                                                                                                                                         cGGACGGGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGc    <  2:851488/150‑1 (MQ=255)
                                                                                                                                         cGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGc    <  2:13817/150‑1 (MQ=255)
                                                                                                                                         cGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGc    <  2:911432/150‑1 (MQ=255)
                                                                                                                                          ggACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCa   <  1:708032/150‑1 (MQ=255)
                                                                                                                                            aCCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAt  <  2:489887/149‑1 (MQ=255)
                                                                                                                                            aCCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAt  >  1:489887/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
ACGTGGCGGGTGATCCGGGCTGGTACCAGGTCTCGCTCGCCGTGCGCCCGCACTTCAAGTACATGGGCGCGAACTTCGAAATTTCCCTGGTAGGCCGGCTCGACACGCAATGAGTGGTTTCTTCGATCGCTTGATGGCGGACCCGTCCGCCATCCGGGAGGCGTCTCGCCAGGAAAGCGCTTCGCACAAGTTCGCGGCCATCAAGTGCCACCTCGAAACGTTGCTGAACGCTCGCCAAGGTTGCTCCCAGAGCAGCCCCGAGCTGGGCCTGCGCGACTTCAACGGGCAT  >  NC_002947/4604108‑4604396

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