Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,508 T→A 9.1% *221C (TGA→TGT) ‡ 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,5080TA 9.1% 124.4 / 4.3 44*221C (TGA→TGT) ‡fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (19/21);  new base A (2/2);  total (21/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 = 5.20e-01

GGGGAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCT  >  NC_002947/243359‑243656
                                                                                                                                                     |                                                                                                                                                    
ggggAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGt                                                                                                                                                      >  1:1036133/1‑150 (MQ=255)
  ggAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCa                                                                                                                                                    <  2:394822/150‑1 (MQ=255)
     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                                                                                           >  1:1356212/1‑125 (MQ=255)
                                                                                    tGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAg                                                                                           <  2:1356212/125‑1 (MQ=255)
                                                                                     gCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                                 <  1:1349265/150‑1 (MQ=255)
                                                                                             gAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                                 <  2:1293069/142‑1 (MQ=255)
                                                                                             gAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                                 >  1:1293069/1‑142 (MQ=255)
                                                                                                      cacaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAg                                                                                                                                                   <  1:1356795/51‑1 (MQ=255)
                                                                                                      cacaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAg                                                                                                                                                   >  2:1356795/1‑51 (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)
                                                                                                                                                     |                                                                                                                                                    
GGGGAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCT  >  NC_002947/243359‑243656

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: