Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,907,691 T→A 10.7% Q335L (CAG→CTG)  hasE ← secretion metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,907,6910TA10.7% 98.9 / 2.5 37Q335L (CAG→CTG) hasEsecretion metalloprotease
Reads supporting (aligned to +/- strand):  ref base T (15/18);  new base A (2/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCAGCCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTC  >  NC_002947/2907547‑2907827
                                                                                                                                                 |                                                                                                                                         
tcgAGCAGCCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCAGccg                                                                                                                                       >  4:217142/1‑147 (MQ=255)
              ggcggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCAGCCGCTGCAGGTCGatgc                                                                                                                         <  3:217142/149‑4 (MQ=255)
               gcggAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCgg                                                                                                                       >  2:223740/1‑149 (MQ=255)
                         tttttACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTc                                                                                                              >  1:211691/1‑148 (MQ=255)
                             tACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATg                                                                                                          <  4:17244/148‑1 (MQ=255)
                               cccGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATc                                                                                                       >  4:21569/1‑149 (MQ=255)
                               cccGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATc                                                                                                       <  2:88414/149‑1 (MQ=255)
                                ccggccggGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCTCATCGACCTGCAGCGGCTTGTCGA‑GCGGTACCACTTCCATGATCt                                                                                                      <  6:107202/149‑1 (MQ=255)
                                 cggccggGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCtt                                                                                                     <  6:143093/149‑1 (MQ=255)
                                 cggccggGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCtt                                                                                                     <  4:9998/149‑1 (MQ=255)
                                               cGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTc                                                                                                                              <  1:211289/111‑1 (MQ=255)
                                               cGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTc                                                                                                                              >  2:211289/1‑111 (MQ=255)
                                                ggCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGga                                                                                       <  2:211691/148‑1 (MQ=255)
                                                ggCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGga                                                                                       <  2:225499/148‑1 (MQ=255)
                                                    tCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGatgatg                                                                                  <  1:223740/149‑1 (MQ=255)
                                                                tGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTg                                                                      >  2:138488/1‑149 (MQ=255)
                                                                  aaGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTgg                                                                     <  1:138488/148‑1 (MQ=255)
                                                                              gCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCGAATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGc                                                         >  1:55059/1‑149 (MQ=255)
                                                                              gCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGcc                                                        >  4:302265/1‑149 (MQ=255)
                                                                                 aTGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTg                                                     <  3:21569/149‑1 (MQ=255)
                                                                                    tCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGAc                                                  <  6:4191/149‑1 (MQ=255)
                                                                                      cACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCa                                                >  1:198402/1‑149 (MQ=255)
                                                                                       aCGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCAt                                               >  1:155216/1‑149 (MQ=255)
                                                                                        cGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATg                                              <  2:55059/149‑1 (MQ=255)
                                                                                               gcaagggACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGatgatg                                                                                  <  1:256423/102‑1 (MQ=255)
                                                                                               gCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGatgatg                                                                                  >  2:256423/1‑106 (MQ=255)
                                                                                                        cATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGc                              <  1:233340/149‑1 (MQ=255)
                                                                                                                 g‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCg                     >  6:243736/1‑149 (MQ=255)
                                                                                                                                ggATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTc        >  8:86625/1‑148 (MQ=255)
                                                                                                                                ggATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCa       >  4:299268/1‑149 (MQ=255)
                                                                                                                                  atcatcGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATa     <  3:299268/149‑1 (MQ=255)
                                                                                                                                   tcatcGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCa       >  6:208681/1‑146 (MQ=255)
                                                                                                                                   tcatcGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTAGATCGGCGAGCGGATCACCGTGTGGGTCACTTCa       <  5:208681/146‑1 (MQ=255)
                                                                                                                                     atcgcatcgACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCAATGATCTTGGCCCCGGGCTGGATGATGCCTCCGACCGTTGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCTGATCACCGTGTGGGTCACTTCATAGTc  <  3:193098/149‑1 (MQ=255)
                                                                                                                                      tcgcatcgACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTc  >  5:304398/1‑148 (MQ=255)
                                                                                                                                         gctcgACCAGCCGCTGCAGGTCGATGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATa     >  7:242121/3‑143 (MQ=255)
                                                                                                                                         gctcgACCAGCCGCTGCAGGTCGATGCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATa     <  8:242121/141‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
TCGAGCAGCCGGTCGGCGGAAATGGTTTTTACCCGGCCGGGAATGTTCGGCGTCTGCGCATGGTTGAAGGCCGGGAAGGCGATGTCCACGGCAAGGCCTGGGACCATCTTGTCG‑ATCGCCTGCACCGGGATCATCGCATCGACCTGCAGCGGCTGGTCGA‑GCGGTACCACTTCCATGATCTTGGCCCCGGGCTGGATGATGCCGCCGACCGTGGCGATGCTCAGGCCCTGGACCATGCCGTCGATCGGCGAGCGGATCACCGTGTGGGTCACTTCATAGTC  >  NC_002947/2907547‑2907827

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