Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,902,230 G→A 8.0% R135H (CGT→CAT)  PP_3445 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,902,2300GA 8.0% 139.9 / 2.1 50R135H (CGT→CAT) PP_3445hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base A (2/2);  total (23/27)
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

TCAAGAACAACCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGACCTGACCACACAGTTTGC  >  NC_002947/3902083‑3902374
                                                                                                                                                   |                                                                                                                                                
tCAAGAACAACCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGt                                                                                                                                                 >  8:49739/1‑149 (MQ=255)
tCAAGAACAACCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCACCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGt                                                                                                                                                 >  8:274792/1‑149 (MQ=255)
tCAAGAACAACCCGGGCACCGCCGAGCACACGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGt                                                                                                                                                 >  5:98540/1‑149 (MQ=255)
  aaGAACAACCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTg                                                                                                                                                <  2:173462/148‑1 (MQ=255)
     aacaacCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCGGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTgg                                                                                                                                            >  1:187182/1‑149 (MQ=255)
     aacaacCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTgg                                                                                                                                            <  1:45430/149‑1 (MQ=255)
     aacaacCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGAGAAATCGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTgg                                                                                                                                            <  1:96320/149‑1 (MQ=255)
           ccGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGc                                                                                                                                      <  7:179901/149‑1 (MQ=255)
               gCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGt                                                                                                                                  >  5:48916/1‑149 (MQ=255)
                cACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTg                                                                                                                                 >  4:210680/1‑149 (MQ=255)
                      cgAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCCGCGCCTTCCGTGCTGGCCCGGCCGGTGGCGCCg                                                                                                                           >  6:43257/1‑149 (MQ=255)
                            cgcgCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATg                                                                                                                     <  8:1400/149‑1 (MQ=255)
                                    gCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGa                                                                                                              >  8:209140/1‑148 (MQ=255)
                                      tGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAg                                                                                                           >  8:51283/1‑149 (MQ=255)
                                        cccGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCACCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGAc                                                                                                         <  7:274792/149‑1 (MQ=255)
                                           gTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCCCCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGAcgc                                                                                                       >  5:101978/1‑148 (MQ=255)
                                              cACCAAGGACGAGGGGGAATCCGCCGAACGTGAAAACGACGCCGCCCCACCCGCCGCCTGGCTGCTACCCGATACACCCTTGCGCCTGATCAGCGACTTCGGTGCTTGCACGGCCGGTGGCGCCGGCGATTGCAACGACAGACGCGaa                                                                                                    <  2:38521/148‑1 (MQ=255)
                                                 cAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGa                                                                                                          >  6:270740/1‑139 (MQ=255)
                                                 cAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGa                                                                                                          <  5:270740/139‑1 (MQ=255)
                                                          gCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCctgc                                                                                       <  3:2241/149‑1 (MQ=255)
                                                                aGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCAc                                                                                 >  3:241172/1‑149 (MQ=255)
                                                                  ccgccgAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCActg                                                                               <  7:209140/149‑1 (MQ=255)
                                                                  ccgccgAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCActg                                                                               <  3:210680/149‑1 (MQ=255)
                                                                    gccgAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCActgct                                                                             <  3:204570/149‑1 (MQ=255)
                                                                         aCGTGAAACCGACCCCGCCCCACCCGGCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGa                                                                        <  6:101978/149‑1 (MQ=255)
                                                                            tGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCActgct                                                                             >  2:42257/1‑141 (MQ=255)
                                                                            tGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCActgct                                                                             <  1:42257/141‑1 (MQ=255)
                                                                             gAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTg                                                                    >  3:173166/1‑149 (MQ=255)
                                                                                 aCGACCCCGCCCCACCCGCCGCCTGGCGGCCACCCGATCCACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGAtt                                                                <  3:176703/149‑1 (MQ=255)
                                                                                        ctccacACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCTACTTCCATGCTGGCACGGCCGGTGGAGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCATGCTGTAACTGCTGCTGACCTGGATTGaacc                                                           >  7:54783/6‑145 (MQ=255)
                                                                                        cGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGCCTTCCATGCTGGCACGGCCGGTGGCGCCGGCGATGGCACCGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGCAAc                                                           >  7:50667/1‑147 (MQ=255)
                                                                                            ccACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCgg                                                     <  6:48916/149‑1 (MQ=255)
                                                                                             cACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGa                                                    >  5:211110/1‑149 (MQ=255)
                                                                                                    cgccTGGCTGCCACCCGCTACCCCCTTGCGCCTGATCAGCGACTTCCATGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAAcc                                             <  8:50667/149‑1 (MQ=255)
                                                                                                    cgccTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAAcc                                             <  1:33174/149‑1 (MQ=255)
                                                                                                    cgccTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAAcc                                             <  6:211110/149‑1 (MQ=255)
                                                                                                    cgccTGGCTGACACCCTATACACCATTGCGCCTGATCAGCGACTTCCATGCTTGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACTCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCTTACTCaaac                                             <  8:54783/149‑3 (MQ=255)
                                                                                                         ggCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTAc                                        <  7:49739/149‑1 (MQ=255)
                                                                                                          gCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCGGGATTGAAACCAGCGGACTCAACCTTTa                                         <  4:241172/147‑1 (MQ=255)
                                                                                                           cTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGc                                      <  2:187182/149‑1 (MQ=255)
                                                                                                                  ccGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGccacccac                               >  8:118097/1‑149 (MQ=255)
                                                                                                                        cacCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTgaag                         <  4:173166/149‑1 (MQ=255)
                                                                                                                        cacCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTg                            >  5:297577/1‑146 (MQ=255)
                                                                                                                         acCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTgaag                         >  8:159358/1‑148 (MQ=255)
                                                                                                                                  ccTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGACCTGa               <  7:51283/149‑1 (MQ=255)
                                                                                                                                  ccTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGACCTGa               <  5:43257/149‑1 (MQ=255)
                                                                                                                                         aGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGACCTGACcacac         <  7:47137/148‑1 (MQ=255)
                                                                                                                                               ttCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGACCTGACCACACAGTTTGc  >  7:118623/1‑149 (MQ=255)
                                                                                                                                               ttCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGACCTGACCACACAGTTTGc  >  3:129228/1‑149 (MQ=255)
                                                                                                                                                  cGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGAc                   >  4:8820/1‑129 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCAAGAACAACCCGGGCACCGCCGAGCACGCGCCTGGCTGCCCGTTCACCAAGGACGAGCGGGAAGCCGCCGAACGTGAAAACGACCCCGCCCCACCCGCCGCCTGGCTGCCACCCGATACACCCTTGCGCCTGATCAGCGACTTCCGTGCTGGCACGGCCGGTGGCGCCGGCGATGGCAACGACAGACGCGAACAGCAGCGCCTGCTGTCACTGCTGCTGACCTGGATTGAAACCAGCGGACTCAACCTTTACGCCACCCACCTGAAGAAAGACCTGACCACACAGTTTGC  >  NC_002947/3902083‑3902374

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