Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,818,497 C→T 7.2% R213R (CGC→CGT)  PP_1622 → M23B subfamily metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,818,4970CT 7.2% 159.7 / 2.2 55R213R (CGC→CGT) PP_1622M23B subfamily metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base T (2/2);  total (28/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 = 2.84e-01

GCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCACGTCGTT  >  NC_002947/1818349‑1818645
                                                                                                                                                    |                                                                                                                                                    
gCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGc                                                                                                                                                      <  8:182929/149‑1 (MQ=255)
      atatCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTg                                                                                                                                                >  8:156566/1‑149 (MQ=255)
         tCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGtt                                                                                                                                              <  5:265022/148‑1 (MQ=255)
                  aTTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAg                                                                                                                                    <  5:209729/149‑1 (MQ=255)
                  aTTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGCTCATAACCGCAGGCTGTTGGTTCGGGAg                                                                                                                                    <  2:251579/149‑1 (MQ=255)
                      gggACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAgggg                                                                                                                                 >  3:231628/1‑148 (MQ=255)
                       ggACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAgg                                                                                                                                                   >  6:207138/1‑129 (MQ=255)
                         aCAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGgcagc                                                                                                                             <  4:231628/149‑1 (MQ=255)
                           aGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGcagcag                                                                                                                           <  1:137142/149‑1 (MQ=255)
                               tGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTc                                                                                                                        >  8:129456/1‑148 (MQ=255)
                                      gctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCTAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTTCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcag                                                                                                                 <  1:110429/148‑1 (MQ=255)
                                       ctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagg                                                                                                                <  3:168169/148‑1 (MQ=255)
                                           gTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGt                                                                                                           <  7:129456/149‑1 (MQ=255)
                                           gTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGt                                                                                                           >  4:126402/1‑149 (MQ=255)
                                                 tGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTg                                                                                                     <  3:126402/149‑1 (MQ=255)
                                                  ggTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGc                                                                                                    >  2:189507/1‑149 (MQ=255)
                                                         tggcTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGTGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTTGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATg                                                                                             <  8:116882/149‑1 (MQ=255)
                                                             ttATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGt                                                                                         >  4:14967/1‑149 (MQ=255)
                                                                 gCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTAc                                                                                     >  6:211208/1‑149 (MQ=255)
                                                                        gTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCt                                                                                                   <  8:150475/128‑1 (MQ=255)
                                                                                     cTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGc                                                                 >  4:179230/1‑149 (MQ=255)
                                                                                               cTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTc                                                         >  1:50755/1‑147 (MQ=255)
                                                                                                  atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCgaga                                                     <  3:14967/148‑1 (MQ=255)
                                                                                                  atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGCTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGAt                                                    >  7:27314/1‑149 (MQ=255)
                                                                                                         tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGt                                                                                         <  5:272700/105‑1 (MQ=255)
                                                                                                         tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGt                                                                                         >  6:272700/1‑105 (MQ=255)
                                                                                                           aaGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcgccgc                                             >  2:263816/1‑147 (MQ=255)
                                                                                                              cacaTCGAAACCTCCGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCTATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                                        <  7:230481/149‑1 (MQ=255)
                                                                                                                         cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGt                               <  5:144926/147‑1 (MQ=255)
                                                                                                                         cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGt                               >  6:144926/1‑147 (MQ=255)
                                                                                                                         cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCg                             <  1:263816/149‑1 (MQ=255)
                                                                                                                         cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCg                             <  5:211208/149‑1 (MQ=255)
                                                                                                                         cTACGTCAGTGCCTACGGTCATAACCGCAGGCTCTTGGTTCGGGACGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATAGGGTGAAGCTGCATTTCGAGATTAGCCGCCAGGGCAAACCCGTc                              >  4:126657/1‑148 (MQ=255)
                                                                                                                          tACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTc                                                                                                                        >  7:479/1‑57 (MQ=255)
                                                                                                                          tACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTc                                                                                                                        <  8:479/57‑1 (MQ=255)
                                                                                                                          tACGTCAGTGACTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGa                            >  6:110347/1‑149 (MQ=255)
                                                                                                                               cAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCAc                       <  7:156566/149‑1 (MQ=255)
                                                                                                                                 gTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTc                     <  2:50755/149‑1 (MQ=255)
                                                                                                                                  tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTc                     <  8:51148/148‑1 (MQ=255)
                                                                                                                                  tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTc                     >  4:115273/1‑148 (MQ=255)
                                                                                                                                  tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTc                     >  7:51148/1‑148 (MQ=255)
                                                                                                                                   gCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGAt                           >  3:237876/1‑141 (MQ=255)
                                                                                                                                   gCCTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGAt                           <  4:237876/141‑1 (MQ=255)
                                                                                                                                    ccTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGGAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                                        <  8:104270/127‑1 (MQ=255)
                                                                                                                                    ccTGCGGTTATGACCGTAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                                        >  7:104270/1‑127 (MQ=255)
                                                                                                                                    ccTACGGTCATAACCGCAGGGTGTTGGTTCGGGAGGGGCAGCAGGTCTTGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTTCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCa                   <  3:126657/148‑1 (MQ=255)
                                                                                                                                     cTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGt                 >  5:89806/1‑149 (MQ=255)
                                                                                                                                     cTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCAGCCAGGGCAAACCCGTCGATCCACTCCAGt                 >  5:150422/1‑149 (MQ=255)
                                                                                                                                      tACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAATTGGGGTCTACGGGCACAGAGCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCCCTCCAGtt                >  7:75722/1‑149 (MQ=255)
                                                                                                                                        cGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGtt                <  3:179230/147‑1 (MQ=255)
                                                                                                                                          gTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTg            >  3:47323/1‑149 (MQ=255)
                                                                                                                                          gTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTg            <  8:27314/149‑1 (MQ=255)
                                                                                                                                           tCATAACCGCAGGCTGTTGGTTCGGGAGTGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGCGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGc           <  8:47273/149‑1 (MQ=255)
                                                                                                                                             aTAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCa         >  4:95987/1‑149 (MQ=255)
                                                                                                                                             aTAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCa         >  8:78105/1‑149 (MQ=255)
                                                                                                                                                   gCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCAcgtcgt   >  5:20514/1‑149 (MQ=255)
                                                                                                                                                    cAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCACGTCGtt  <  5:110347/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGCCACGTCGTT  >  NC_002947/1818349‑1818645

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