Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,316,584 G→T 8.7% intergenic (+32/+4) PP_2033 → / ← PP_2034 hypothetical protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,316,5840GT 8.7% 156.8 / 2.5 46intergenic (+32/+4)PP_2033/PP_2034hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base T (2/2);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

GCCCTGAAGGGCTTGCGCAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGCGACGAAAATCCCGTTGGGCAGGCGCTG  >  NC_002947/2316457‑2316716
                                                                                                                               |                                                                                                                                    
gCCCTGAAGGGCTTGCGCAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACa                                                                                                                              <  1:372143/136‑1 (MQ=255)
    tGAAGGGCTTGCGCAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGgc                                                                                                                          >  1:66245/1‑136 (MQ=255)
     gAAGGGCTTGCGCAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGgcg                                                                                                                         <  1:410703/136‑1 (MQ=255)
              gtgcAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTGGTACCGGGCCTACAGGGCGCGCGCCAGc                                                                                                                >  2:428015/3‑136 (MQ=255)
              gcgcAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGc                                                                                                                >  1:146833/1‑136 (MQ=255)
              gcgcAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGc                                                                                                                >  1:187310/1‑136 (MQ=255)
              gcgcAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGc                                                                                                                >  2:557400/1‑136 (MQ=255)
              gcgcAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGc                                                                                                                >  1:366546/1‑136 (MQ=255)
              gcgcAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGc                                                                                                                >  1:414408/1‑136 (MQ=255)
               cgcAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCt                                                                                                               >  2:363342/1‑136 (MQ=255)
               cgcAACAGTCAGCTTGATGTGCTGTTCGTCAATCGGCCTTGGGCAGCCGGGCGCGTGGTGAAATCCAACCCGGCGCAATGAATGACTGATCATCGCCGGGCAATGTTTTACCGGGCCTACAGGGCGCGCCCCAGcg                                                                                                               >  2:617766/1‑135 (MQ=255)
                       tgtgCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTc                                                                                                       <  1:272351/136‑1 (MQ=255)
                       tgtgCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTc                                                                                                       >  1:70205/1‑136 (MQ=255)
                       ggtgCTGGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTc                                                                                                       <  2:668789/135‑1 (MQ=255)
                            ttGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTcgcgcg                                                                                                  >  2:117381/1‑136 (MQ=255)
                            ttGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTcgcgcg                                                                                                  >  1:290286/1‑136 (MQ=255)
                            ttGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTcgcgcg                                                                                                  >  1:468812/1‑136 (MQ=255)
                                 gtgCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTc                                                                                             <  1:142819/136‑1 (MQ=255)
                                 gtgCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTGGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTc                                                                                             <  1:428015/136‑1 (MQ=255)
                                    cTGTTCGTCAATCGGCCGGGGGCAGGCGGTCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCTGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTgg                                                                                          >  1:64975/1‑136 (MQ=255)
                                    cTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCTGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTgg                                                                                          >  1:228397/1‑136 (MQ=255)
                                    cTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGTGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTgg                                                                                          >  1:629913/1‑136 (MQ=255)
                                    cTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTgg                                                                                          <  1:637655/136‑1 (MQ=255)
                                    cTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTgg                                                                                          >  2:44766/1‑136 (MQ=255)
                                              aTCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACAt                                                                                <  1:654492/136‑1 (MQ=255)
                                                  gCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTcc                                                                            >  2:225509/1‑136 (MQ=255)
                                                  gCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTcc                                                                            >  2:32624/1‑136 (MQ=255)
                                                  gCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTcc                                                                            >  2:26523/1‑136 (MQ=255)
                                                      ggggCAGGCGGTCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCTGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACTTCTCCtgtg                                                                        <  2:64975/136‑1 (MQ=255)
                                                      ggggCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCTGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCtgtg                                                                        <  2:228397/136‑1 (MQ=255)
                                                      ggggCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCtgtg                                                                        <  1:363342/136‑1 (MQ=255)
                                                        ggCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGcc                                                                      <  1:117381/136‑1 (MQ=255)
                                                             gcgggcgCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGt                                                                 <  2:151661/136‑1 (MQ=255)
                                                               gggcgCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTca                                                               <  2:603143/136‑1 (MQ=255)
                                                                     gtggtgCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAgcagca                                                         >  2:508861/1‑136 (MQ=255)
                                                                             tGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTcagc                                                 <  1:295144/136‑1 (MQ=255)
                                                                                  aCCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCAt                                            <  1:44766/136‑1 (MQ=255)
                                                                                      ggCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCtgt                                                                         >  2:437412/1‑103 (MQ=255)
                                                                                      ggCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCtgt                                                                         <  1:437412/103‑1 (MQ=255)
                                                                                           aatgaatgACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAAcaccac                                   >  2:630802/1‑136 (MQ=255)
                                                                                                   aCTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGcga                           <  2:414408/136‑1 (MQ=255)
                                                                                                           tggcctggcAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGCGACGAAAATc                   <  1:617766/136‑1 (MQ=255)
                                                                                                                tggcAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGCGACGAAAATCCCGtt              >  1:349024/1‑136 (MQ=255)
                                                                                                                    aaTGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGCGACGAAAATCCCGTTGGGc          >  2:2629/1‑136 (MQ=255)
                                                                                                                    aaTGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGCGACGAAAATCCCGTTGGGc          >  2:70551/1‑136 (MQ=255)
                                                                                                                            aCCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGCGACGAAAATCCCGTTGGGCAGGCGCTg  <  2:468812/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                                    
GCCCTGAAGGGCTTGCGCAACAGTGTGCTTGATGTGCTGTTCGTCAATCGGCCGGGGGCAGGCGGGCGCGTGGTGCATGCCAACCCGGCGCAATGAATGACTGATCCTGGCCTGGCAATGTTGTACCGGGCCTACAGGGCGCGCGCCAGCTGTTCAGTCGCGCGCCTTCTGGCTGAATACATCTCCTGTGCCAGGGTCAGCAGCACCGTCAGCAGCATGAACACCACCGATGCGACGAAAATCCCGTTGGGCAGGCGCTG  >  NC_002947/2316457‑2316716

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