Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 101,821 T→A 9.1% L214* (TTG→TAG)  prlC → oligopeptidase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947101,8210TA 9.1% 120.8 / 3.2 44L214* (TTG→TAG) prlColigopeptidase A
Reads supporting (aligned to +/- strand):  ref base T (20/20);  new base A (2/2);  total (22/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 = 9.30e-01

CCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGACAATGGCCCGGT  >  NC_002947/101681‑101962
                                                                                                                                            |                                                                                                                                             
ccGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCAccc                                                                                                                                       <  4:75863/149‑1 (MQ=255)
     tctcCAATCAGCTGCTGGACGCCCCCCAGGCCGGGACCAAGCAGGTTACCGACGAAGCTGCGCTGGCCGGCGTGACCGTTTCGGCCATGGCGCATATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGa                                                                                                                                   <  8:198028/148‑1 (MQ=255)
         cAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTAGATCACCCTTGAATTcc                                                                                                                              >  5:27810/1‑149 (MQ=255)
         cAATCAGCTGCTGGACGACAACCACCCCTGGACCAAGCATGTTACAGACGAAGCTGCGCTGGACGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGATGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTAGATCACCCTTGAATTcc                                                                                                                              >  5:33110/1‑149 (MQ=255)
                         gCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGACGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGt                                                                                                              >  8:287530/1‑149 (MQ=255)
                          ccACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGAGGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGtg                                                                                                             <  8:149954/149‑1 (MQ=255)
                                  tccTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                                           >  6:78932/2‑127 (MQ=255)
                                  gCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCa                                                                                                                           <  5:78932/127‑1 (MQ=255)
                                       tcccAAGAATTTTACCGACTACGCGGCGCTGGCCGGCGTGACCGTTTCGGCAAAGGCGCAGTTGGCCGCTGCAGCCCAGGCAAAGGGCCTGGACGGCTGGTAGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCa                                                                                                <  6:27810/147‑1 (MQ=255)
                                       gaacAAGCTTGTTACCGTCTACGCTGCGCTGGCCGGCGTGACCGATGCGTCCAAGTGGGTGTTGGCCGCGGCACCCCAGTCAATGGGGGTGGACGGCTGGTAGATCACCCTTGAATTCCCCAGCTACTACGCGGTTATGACCTACGCCa                                                                                                <  6:33110/146‑1 (MQ=255)
                                       gACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGc                                                                                                  <  7:185901/147‑1 (MQ=255)
                                               aTGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACcgc                                                                                        >  3:102541/1‑149 (MQ=255)
                                                    aCCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCGGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGccc                                                                                    <  7:268105/148‑1 (MQ=255)
                                                      cgacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgc                                                                                 >  1:214302/1‑149 (MQ=255)
                                                       gacgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCACGTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcg                                                                                >  5:311199/1‑149 (MQ=255)
                                                        acgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGAc                                                                                                        <  6:86072/124‑1 (MQ=255)
                                                        acgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGAc                                                                                                        >  5:86072/1‑124 (MQ=255)
                                                        acgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcgc                                                                               >  8:291812/1‑149 (MQ=255)
                                                        acgaAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTgcgc                                                                               >  7:274235/1‑149 (MQ=255)
                                                           aaGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGaa                                                                            <  5:153957/149‑1 (MQ=255)
                                                                         ggccTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATAcgccgc                                                              >  3:41923/1‑149 (MQ=255)
                                                                           ccTGACCGATTCGGCCAAGGCGCAGACGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGGGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCt                                                            <  7:287530/149‑1 (MQ=255)
                                                                              gACCGATCCGCCCAAGGCGCAGATGGCCGCGGCAGCCCCGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACt                                                         <  1:354642/149‑1 (MQ=255)
                                                                               aCCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAaga                                                                          >  4:349830/1‑131 (MQ=255)
                                                                               aCCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAaga                                                                          <  3:349830/131‑1 (MQ=255)
                                                                                   aTTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCAcc                                                    >  7:141408/1‑149 (MQ=255)
                                                                                           aaGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTc                                            <  4:41923/149‑1 (MQ=255)
                                                                                               cgcAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGAc                                        <  2:214302/149‑1 (MQ=255)
                                                                                                     tGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAggg                                   <  4:102541/148‑1 (MQ=255)
                                                                                                         cGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCg                               >  8:287637/1‑148 (MQ=255)
                                                                                                           cTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATTACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAAt                            <  7:18294/149‑1 (MQ=255)
                                                                                                               aGCCCAGGCAAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCg                        <  8:141408/149‑1 (MQ=255)
                                                                                                                gtccAGGCCCTTGGCCTGGGCGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCGCCCGCGCCTCCGACCAGGGCCCGAATGCCgg                       >  4:12409/3‑149 (MQ=255)
                                                                                                                   cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTcc                                           <  7:215669/126‑1 (MQ=255)
                                                                                                                   cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTcc                                           >  8:215669/1‑126 (MQ=255)
                                                                                                                   cAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGc                      >  5:117049/1‑147 (MQ=255)
                                                                                                                       ccAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTc                >  3:318396/1‑149 (MQ=255)
                                                                                                                         aaGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGC‑‑CCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGt                  <  2:184680/143‑1 (MQ=255)
                                                                                                                         aaGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGC‑‑CCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGt                  >  1:184680/1‑143 (MQ=255)
                                                                                                                            ggCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGACaa           >  4:21702/1‑149 (MQ=255)
                                                                                                                                  gACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCg                        >  6:107364/1‑130 (MQ=255)
                                                                                                                                  gACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCg                        <  5:107364/130‑1 (MQ=255)
                                                                                                                                    cGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGACAATGGCCCgg   >  4:243778/1‑149 (MQ=255)
                                                                                                                                     ggCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGACAATGGCCCGGt  <  8:274235/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
CCGTTTCTCCAATCAGCTGCTGGACGCCACCCAGGCCTGGACCAAGCATGTTACCGACGAAGCTGCGCTGGCCGGCCTGACCGATTCGGCCAAGGCGCAGATGGCCGCTGCAGCCCAGGCCAAGGGCCTGGACGGCTGGTTGATCACCCTTGAATTCCCCAGCTACTACGCGGTGATGACCTACGCCAGCGACCGCGCCCTGCGCGAAGAGCTATACGCCGCCTACTGCACCCGCGCCTCCGACCAGGGCCCGAATGCCGGCCAGTTCGACAATGGCCCGGT  >  NC_002947/101681‑101962

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