Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,811 G→C 12.1% R143P (CGG→CCG)  PP_0891 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,8110GC12.1% 88.3 / 5.3 33R143P (CGG→CCG) PP_0891hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/13);  new base C (2/2);  total (18/15)
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.11e-01

CACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATC  >  NC_002947/1033665‑1033950
                                                                                                                                                  |                                                                                                                                           
cacaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCgcggcg                                                                                                                                          >  1:433645/1‑150 (MQ=255)
 acaGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGgcgc                                                                                                                                         >  1:207899/1‑150 (MQ=255)
          gCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtg                                                                                                                                >  2:99566/1‑150 (MQ=255)
          gCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtg                                                                                                                                <  2:207899/150‑1 (MQ=255)
            acaGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAggtggt                                                                                                                              >  1:360444/1‑150 (MQ=255)
                 tCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGAccc                                                                                                                         <  2:360444/150‑1 (MQ=255)
                    aCCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTg                                                                                                                      >  1:482714/1‑150 (MQ=255)
                     ccTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGa                                                                                                                     <  1:377693/150‑1 (MQ=255)
                             cggtCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTgg                                                                                                             >  1:384730/4‑150 (MQ=255)
                                           gCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGaa                                                                                               <  1:77344/150‑1 (MQ=255)
                                             cGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGa                                                                                             >  2:18259/1‑150 (MQ=255)
                                             cGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGa                                                                                             <  2:265349/150‑1 (MQ=255)
                                                      acaGATCGAATTTCGCCTGCCGCTGACCGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTTGCCACCCGCCTGGATGTGCGACGCTTCCAGGTGGTGACCCTTGATCCTCTGGTGCTGATCGCCGAAGATTATGGCtt                                                                                    >  2:281397/1‑150 (MQ=255)
                                                      acaGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCtt                                                                                    <  1:99566/150‑1 (MQ=255)
                                                                            cTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTg                                                              <  2:11498/150‑1 (MQ=255)
                                                                                    aaCCCGGCATGGGCAGTCGCACAGCTACAATGTCCTGTTGCTGGCAACCCGACTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGTGGCTGGAAAGCCTGCGCAAGtt                                                      >  2:277446/2‑150 (MQ=255)
                                                                                       ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGATCCtt                                                                                                                       >  1:292130/1‑82 (MQ=255)
                                                                                       ccTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCtt                                                                                                                       <  2:292130/82‑1 (MQ=255)
                                                                                          gCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAg                                                >  1:61839/1‑150 (MQ=255)
                                                                                                 cAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCGGGATGCGCGGCGCTCCGAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCCTGCTGCCGGGGCTGTGAAGCCTGCGCAAGTTCGCCAGGTTGAAg                                         >  2:245818/1‑150 (MQ=255)
                                                                                                     cGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCa                                     <  1:548456/150‑1 (MQ=255)
                                                                                                     cGCACAGCTACAATGTCCTGTTGCTGGAATTTCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCa                                     <  1:463402/150‑1 (MQ=255)
                                                                                                      gCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCAt                                    <  2:433645/150‑1 (MQ=255)
                                                                                                                       tGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                         <  1:464966/112‑1 (MQ=255)
                                                                                                                       tGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCaa                                                         >  2:464966/1‑112 (MQ=255)
                                                                                                                        gTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGt                  >  1:597844/1‑150 (MQ=255)
                                                                                                                            cTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAgc              >  1:390445/1‑150 (MQ=255)
                                                                                                                                   cccGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGt                            >  2:162650/1‑129 (MQ=255)
                                                                                                                                   cccGCCTGGATGCGCGGCGCGTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGt                            <  1:162650/129‑1 (MQ=255)
                                                                                                                                     caccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGt                                        <  2:403383/113‑1 (MQ=255)
                                                                                                                                     caccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGt                                        >  1:403383/3‑115 (MQ=255)
                                                                                                                                       ccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGa    <  1:466488/149‑1 (MQ=255)
                                                                                                                                       ccTGGAAGCGCCGCGCATCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGa    >  2:466488/1‑149 (MQ=255)
                                                                                                                                        cTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATc  >  2:387452/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
CACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGTCACCCTGCATGGCCAGTCGCACAGCTACAATGTCCTGTTGCTGGCCACCCGCCTGGATGCGCGGCGCTTCCAGGTGGTGACCCTTGAGCCTCTGGTGCTGCGCGCCGAAGATTTTGGCTTGCTGCCGGGGCTGGAAAGCCTGCGCAAGTTCGCCAGGTTGAAGTCCATCAACCCGTCGGTACCGGTTAGCGCGGTGCTGATC  >  NC_002947/1033665‑1033950

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