Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,209,062 G→C 10.0% A96P (GCG→CCG)  PP_2813 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,209,0620GC10.0% 116.4 / 4.2 40A96P (GCG→CCG) PP_2813BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (20/16);  new base C (2/2);  total (22/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.51e-01

CCGGTGCGCGCGTGGTTGCCGTGGGCGACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAGGCCTCTGGTCGT  >  NC_002947/3208919‑3209209
                                                                                                                                                |                                                                                                                                                    
ccGGTGCGCGCGTGGTTGCCGTGGGCGACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGtt                                                                                                                                                 <  8:102074/149‑1 (MQ=255)
     gacaccgTGGTCGCCCACGGCAACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCa                                                                                                                                            >  2:54227/6‑149 (MQ=255)
        cgcgTGGTTGCCGTGGGCGACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTg                                                                                                                                         <  4:209462/149‑1 (MQ=255)
                    gTGGGCGACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCt                                                                                                                             >  2:234937/1‑149 (MQ=255)
                        gCGACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGAt                                                                                                                          >  8:115389/1‑148 (MQ=255)
                        gCGACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATc                                                                                                                         <  7:113412/149‑1 (MQ=255)
                                         ttgttgtCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCggg                                                                                                          >  3:219819/1‑147 (MQ=255)
                                           gttgtCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGCTCCTGATCACCACCGATGGCGGGGAGt                                                                                                      >  6:312459/1‑149 (MQ=255)
                                           gttgtCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGc                                                                                                             >  4:188367/1‑142 (MQ=255)
                                           gttgtCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGc                                                                                                             <  3:188367/142‑1 (MQ=255)
                                              gtCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGGCTGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTc                                                                                                     <  1:234937/148‑1 (MQ=255)
                                              gtCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTcc                                                                                                    <  4:219819/148‑1 (MQ=255)
                                                   acgacCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGCCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAg                                                                                              >  1:180761/1‑149 (MQ=255)
                                                   acgacCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGCCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAg                                                                                              >  1:180759/1‑149 (MQ=255)
                                                            gCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTgg                                                                                                 <  2:68044/137‑1 (MQ=255)
                                                            gCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTgg                                                                                                 >  1:68044/1‑137 (MQ=255)
                                                             ccTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGcc                                                                                    >  6:222516/1‑149 (MQ=255)
                                                                 cAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCggg                                                                                                          <  3:29190/123‑1 (MQ=255)
                                                                 cAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCggg                                                                                                          >  4:29190/1‑123 (MQ=255)
                                                                                 gTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCgg                                                               >  6:188013/1‑149 (MQ=255)
                                                                                   tCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCt                                                                                                   >  6:263777/1‑112 (MQ=255)
                                                                                   tCGGAACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCt                                                                                                   <  5:263777/112‑1 (MQ=255)
                                                                                      gTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACTCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACt                                                          <  5:291914/149‑1 (MQ=255)
                                                                                                    acCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCAct                                            >  7:132402/1‑149 (MQ=255)
                                                                                                        ttGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCActtctt                                        <  5:11216/149‑1 (MQ=255)
                                                                                                          gACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACttcttc                                       >  4:180668/1‑148 (MQ=255)
                                                                                                           aCGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCa                                      <  7:115389/148‑1 (MQ=255)
                                                                                                            cGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAAt                                    <  5:222516/149‑1 (MQ=255)
                                                                                                            cGGGCGTGTCTATCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAAt                                    <  5:312459/149‑1 (MQ=255)
                                                                                                               gCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAGCAAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGc                                  <  1:54227/149‑1 (MQ=255)
                                                                                                                  tgtCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCg                               >  7:106126/1‑148 (MQ=255)
                                                                                                                      tttCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAgg                          >  2:109729/1‑149 (MQ=255)
                                                                                                                                  cAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACTGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGAGGTAg              >  8:256239/1‑149 (MQ=255)
                                                                                                                                  cAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAg              >  4:204049/1‑149 (MQ=255)
                                                                                                                                  cAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAg              >  1:232153/1‑149 (MQ=255)
                                                                                                                                  cAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAg              <  2:232153/149‑1 (MQ=255)
                                                                                                                                  cAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGCCCACCGAAG‑AAGACAGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAg              >  6:195580/1‑149 (MQ=255)
                                                                                                                                   aaaCGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGACACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAgg             >  4:92566/1‑149 (MQ=255)
                                                                                                                                    aaCGCGGG‑TGGCCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAGGc            <  2:180759/149‑1 (MQ=255)
                                                                                                                                    aaCGCGGG‑TGGCCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAGGc            <  2:180761/149‑1 (MQ=255)
                                                                                                                                               ggCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAGGCCTCTGGTCGt  <  8:132402/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCGGTGCGCGCGTGGTTGCCGTGGGCGACCACGGTGTCGTGTTGTTGTCCGACGACCAGGGCCTGCAGTGGCGCCAGGCGAGTTCGGTACCGCTGTCCACACCATTGACGGGCGTGTCTTTCGTGGATGCCAAACGCGGG‑TGGGCGGTTGGCCACTGGGGCGCGATCCTGATCACCACCGATGGCGGGGAGTCCTGGCAGGTCCAGCGCCTGGCCACCGAAG‑AAGACCGGCCACTGTTCGCGGTGCACTTCTTCAATGCCCGTCAGGGTGTTGCGGTAGGCCTCTGGTCGT  >  NC_002947/3208919‑3209209

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