Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 243,510 A→T 9.3% *221R (TGA→AGA) ‡ 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,5100AT 9.3% 126.8 / 5.1 43*221R (TGA→AGA) ‡fklFKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base T (2/2);  total (21/22)
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.61e-01

GGAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCT  >  NC_002947/243361‑243656
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                                                   >  2:1356795/1‑51 (MQ=255)
                                                                                                    cacaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAg                                                                                                                                                   <  1:1356795/51‑1 (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                                                                               <  2:765039/114‑1 (MQ=255)
                                                                                                         cGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGcccc                                                                               >  1:765039/1‑114 (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            <  2:789281/146‑1 (MQ=255)
                                                                                                                                            gTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGc            >  1:789281/1‑146 (MQ=255)
                                                                                                                                                gCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAtt                  <  1:1296441/136‑1 (MQ=255)
                                                                                                                                                gCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCAtt                  >  2:1296441/1‑136 (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)
                                                                                                                                                     |                                                                                                                                                  
GGAGTGCAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGGCGACCT  >  NC_002947/243361‑243656

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: