Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,512 T→A 8.7% D220V (GAC→GTC)  fkl ← FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947243,5120TA 8.7% 133.7 / 4.2 46D220V (GAC→GTC) fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (21/21);  new base A (2/2);  total (23/23)
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

GTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCA  >  NC_002947/243364‑243660
                                                                                                                                                    |                                                                                                                                                    
gtgtAAGGAATCCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCTCTTCGATCTTGAATACTGCCTGATCGAAGTGCTGCACAACGTACACAAGGAAGAGGGTTCAGTGCAGGCCTGAGGTTGGCGTCAGTc                                                                                                                                                     <  2:505719/150‑1 (MQ=255)
   cAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTGGGCGTCAGTCgac                                                                                                                                                  >  1:1002377/1‑148 (MQ=255)
      ggACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCaa                                                                                                                                                >  2:86361/1‑149 (MQ=255)
      ggACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCaa                                                                                                                                                <  1:86361/149‑1 (MQ=255)
      ggACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGATGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGACTGACGCCAAc                                                                                                                                               <  2:26131/150‑1 (MQ=255)
                         gagaAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGACTGACGCCAAcctcag                                                                                                                                          >  1:1189335/1‑132 (MQ=255)
                         gagaAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGACTGACGCCAAcctcag                                                                                                                                          <  2:1189335/136‑5 (MQ=255)
                           gaAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGaa                                                                                                                          <  2:542418/150‑1 (MQ=255)
                                  gCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAg                                                                                                                   <  2:210644/150‑1 (MQ=255)
                                  gCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAg                                                                                                                   >  2:1349265/1‑150 (MQ=255)
                                                     tgCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCa                                                                                                <  2:459656/150‑1 (MQ=255)
                                                            gcgcTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGgtg                                                                                                             <  2:490145/130‑1 (MQ=255)
                                                            gcgcTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGgtg                                                                                                             >  1:490145/1‑130 (MQ=255)
                                                                ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGAtctc                                                                                                                             <  1:1018186/110‑1 (MQ=255)
                                                                ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGAtctc                                                                                                                             >  2:1018186/1‑110 (MQ=255)
                                                                ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgggg                                                                                                                >  1:496345/1‑123 (MQ=255)
                                                                ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAgggg                                                                                                                <  2:496345/123‑1 (MQ=255)
                                                                ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGcc                                                                                     >  1:889326/1‑150 (MQ=255)
                                                                               tGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAg                                                                                               <  2:1356212/125‑1 (MQ=255)
                                                                               tGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAg                                                                                               >  1:1356212/1‑125 (MQ=255)
                                                                                gCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                                     <  1:1349265/150‑1 (MQ=255)
                                                                                        gAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                                     <  2:1293069/142‑1 (MQ=255)
                                                                                        gAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                                     >  1:1293069/1‑142 (MQ=255)
                                                                                                  acaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAggg                                                                                                        <  1:298909/97‑1 (MQ=255)
                                                                                                  acaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAggg                                                                                                        >  2:298909/1‑97 (MQ=255)
                                                                                                  acaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTa                                                   >  1:1057280/1‑150 (MQ=255)
                                                                                                   caACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTg                                                            >  1:39560/1‑140 (MQ=255)
                                                                                                   caACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTg                                                            <  2:39560/140‑1 (MQ=255)
                                                                                                    aaCGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTAcc                                                 >  1:1018727/1‑150 (MQ=255)
                                                                                                      cGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGcccc                                                                                   >  1:765039/1‑114 (MQ=255)
                                                                                                      cGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGcccc                                                                                   <  2:765039/114‑1 (MQ=255)
                                                                                                                      ggTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCg                               >  1:370727/1‑150 (MQ=255)
                                                                                                                             tgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCa                        >  1:488181/1‑150 (MQ=255)
                                                                                                                                agGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGTCGCCATAGGCTTGCGCTGAAGGAATTTCCAGGCGCCATTTTGCGCCGGTGTGCAtgt                     <  2:1002377/150‑3 (MQ=255)
                                                                                                                                    tggcggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCgg                 >  1:756723/5‑150 (MQ=255)
                                                                                                                                         gTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGc                >  1:789281/1‑146 (MQ=255)
                                                                                                                                         gTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGc                <  2:789281/146‑1 (MQ=255)
                                                                                                                                             gCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAtt                      >  2:1296441/1‑136 (MQ=255)
                                                                                                                                             gCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAtt                      <  1:1296441/136‑1 (MQ=255)
                                                                                                                                             gCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGAc        <  2:889326/150‑1 (MQ=255)
                                                                                                                                               gTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCt      <  2:370727/150‑1 (MQ=255)
                                                                                                                                                  aGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCa                        <  2:403051/129‑1 (MQ=255)
                                                                                                                                                  aGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCa                        >  1:403051/1‑129 (MQ=255)
                                                                                                                                                   gTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTg                                                                                                      <  1:760165/50‑1 (MQ=255)
                                                                                                                                                   gTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTg                                                                                                      >  2:760165/1‑50 (MQ=255)
                                                                                                                                                   gTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCa  >  1:111194/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCTGCCA  >  NC_002947/243364‑243660

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