Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,171,608 T→G 9.5% Q21P (CAG→CCG)  PP_1026 ← M23/M37 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,171,6080TG 9.5% 121.0 / 4.9 42Q21P (CAG→CCG) PP_1026M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base T (19/19);  new base G (2/2);  total (21/21)
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.56e-01

TTGTCCTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAGGCG  >  NC_002947/1171473‑1171750
                                                                                                                                       |                                                                                                                                               
ttGTCCTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGAGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGc                                                                                                                                   <  1:179469/149‑1 (MQ=255)
     cTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCcag                                                                                                                                >  2:110128/1‑147 (MQ=255)
        cTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCcagcagca                                                                                                                           <  5:176853/149‑1 (MQ=255)
                         ccGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAg                                                                                                           >  3:322362/1‑148 (MQ=255)
                         ccGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAg                                                                                                           >  6:114639/1‑148 (MQ=255)
                                 ccGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAgg                                                                                                                >  2:339409/1‑135 (MQ=255)
                                 ccGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAgg                                                                                                                <  1:339409/135‑1 (MQ=255)
                                   gTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAg                                                                                                                 >  2:24453/1‑132 (MQ=255)
                                      gAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGACTG‑GGCGCCGCCGGCCAGAAGCAGAAGGGAATGGGCAAGCATGGGAGCGAGCa                                                                                             >  6:122777/1‑149 (MQ=255)
                                          cGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCTGGTGATGTAGCTGGCCTT‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCAcgcg                                                                                         <  5:122777/149‑1 (MQ=255)
                                            ggCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCggg                                                                                       <  5:114639/149‑1 (MQ=255)
                                            ggCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCggg                                                                                       <  4:322362/149‑1 (MQ=255)
                                               gcTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCa                                                                                    <  7:119348/149‑1 (MQ=255)
                                                            gcggcgGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCAt                                                                                   <  1:41504/137‑1 (MQ=255)
                                                            gcggcgGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCAt                                                                                   >  2:41504/1‑137 (MQ=255)
                                                             cggcggGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCagt                                                                      <  1:264655/149‑1 (MQ=255)
                                                                  gggCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑CGCGCCGCCGGCCAGCAGCAGCAGGGAAAGTGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTga                                                                 <  1:110128/149‑1 (MQ=255)
                                                                      ccAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                          >  1:7517/1‑136 (MQ=255)
                                                                      ccAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                          <  2:7517/136‑1 (MQ=255)
                                                                      ccAGGTCGACGACTGCCACACCGCCGGGCACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAgg                                                             >  7:142727/1‑149 (MQ=255)
                                                                                    gCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGa                                                                                                >  4:193930/1‑100 (MQ=255)
                                                                                           cgccgGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCaa                                                                                                            <  1:221018/81‑1 (MQ=255)
                                                                                           cgccgGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCaa                                                                                                            >  2:221018/1‑81 (MQ=255)
                                                                                           cgccgGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTg                                        >  3:259924/1‑149 (MQ=255)
                                                                                                ggcaCCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTc                                   <  8:142727/146‑1 (MQ=255)
                                                                                                     cGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGAccc                              >  6:47993/1‑149 (MQ=255)
                                                                                                            ttCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAg                       >  5:270109/1‑149 (MQ=255)
                                                                                                              cAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGc                      <  6:270109/148‑1 (MQ=255)
                                                                                                               aGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCaa                    >  2:315970/1‑149 (MQ=255)
                                                                                                                 tgtgCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                          >  5:13375/1‑93 (MQ=255)
                                                                                                                 tgtgCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                          <  6:13375/93‑1 (MQ=255)
                                                                                                                      gggTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCt             <  5:47993/149‑1 (MQ=255)
                                                                                                                             gTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAggg                                                                                                      <  4:101531/53‑1 (MQ=255)
                                                                                                                             gTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAggg                                                                                                      >  3:101531/1‑53 (MQ=255)
                                                                                                                              ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc           <  8:239645/141‑1 (MQ=255)
                                                                                                                              ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc           >  7:239645/3‑143 (MQ=255)
                                                                                                                              ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc           <  4:8897/141‑1 (MQ=255)
                                                                                                                              ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGCGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc           >  3:8897/3‑143 (MQ=255)
                                                                                                                                gCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCACGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAggc   >  4:205833/1‑149 (MQ=255)
                                                                                                                                 cTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAggcg  <  4:259924/149‑1 (MQ=255)
                                                                                                                                  tGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCggg                                                    >  2:296831/1‑98 (MQ=255)
                                                                                                                                  tGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCggg                                                    <  1:296831/98‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                               
TTGTCCTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAGGCG  >  NC_002947/1171473‑1171750

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: ‑