Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 249,248 T→C 8.1% Q200Q (CAA→CAG qmcA ← membrane protease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947249,2480TC 8.1% 137.6 / 2.6 49Q200Q (CAA→CAGqmcAmembrane protease family protein
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base C (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

CACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGATCTCGACGTTGGCCACCTTGA  >  NC_002947/249099‑249396
                                                                                                                                                     |                                                                                                                                                      
cACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCGTATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCt                                                                                                                                                        <  2:497882/150‑1 (MQ=255)
         tcCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGt                                                                                                                                               >  1:750359/1‑150 (MQ=255)
           cTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGttt                                                                                                                                             >  2:684643/1‑150 (MQ=255)
                 caccaACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGAGGCTTGCATCACTTTA‑‑TCCGa                                                                                                                                     >  2:612794/1‑150 (MQ=255)
                    caACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGc                                                                                                                                  >  1:413462/1‑150 (MQ=255)
                              agcagATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCa                                                                                                                        >  1:1075717/1‑150 (MQ=255)
                               gcagATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCAc                                                                                                                       >  2:938444/1‑150 (MQ=255)
                                      tATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGc                                                                                                                >  2:847087/1‑150 (MQ=255)
                                        tCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCAt                                                                                                              <  2:750359/150‑1 (MQ=255)
                                                                  ctctTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCgcgc                                                                                    <  1:847087/150‑1 (MQ=255)
                                                                   tctTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGct                                                                                   <  2:600266/150‑1 (MQ=255)
                                                                       gTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCCTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGTAt                                                                                                              <  1:770800/119‑1 (MQ=255)
                                                                       gTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCCTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCAt                                                                                                              <  1:214748/119‑1 (MQ=255)
                                                                       gTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCCTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCAt                                                                                                              >  2:214748/1‑119 (MQ=255)
                                                                       gTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCAGCTCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCCTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCAt                                                                                                              >  2:770800/1‑119 (MQ=255)
                                                                        tCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGctcggcc                                                                              <  2:823351/150‑1 (MQ=255)
                                                                           ccGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGctcggcctcg                                                                           >  2:659268/1‑150 (MQ=255)
                                                                            cGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGctcggcctcgg                                                                          >  1:387871/1‑150 (MQ=255)
                                                                                caTCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTg                                                                      <  1:706509/149‑1 (MQ=255)
                                                                                    gTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTggcgg                                                                  <  1:612794/150‑1 (MQ=255)
                                                                                     tCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTggcggg                                                                 >  1:1053036/1‑150 (MQ=255)
                                                                                     tCCCCAGTGGCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTggcggg                                                                 <  2:1024076/150‑1 (MQ=255)
                                                                                     cccccAGTGTCTCCATATAGCGCAACTGCATCGCCCCCGATTCCTTACTC‑‑CTGACGCTCTCTTTGCTGATTTTTCCAACCGACGCCTGCAGTTCACCCTCGCCATGAATCACCCTGGCCCTACCTTCGCGCTCGGCCTCGGCCTggcggg                                                                 <  1:684643/149‑1 (MQ=255)
                                                                                        ccAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGa                                                              >  1:60964/1‑150 (MQ=255)
                                                                                                 gCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTa                                                     <  2:1075717/150‑1 (MQ=255)
                                                                                                   atatAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAAGCACCTTGGCCCGACGTTCGCGCTCGGCCCCGGCCTGGCGGGCGCTGGCGCGTAcc                                                   >  1:359307/1‑150 (MQ=255)
                                                                                                         cgcAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGAc                                             >  1:360088/1‑150 (MQ=255)
                                                                                                           cAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCCCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTc                                           <  1:541111/150‑1 (MQ=255)
                                                                                                            aaCTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCg                                          >  1:563173/1‑150 (MQ=255)
                                                                                                             aCTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATc                                                                                                                                                  <  2:1102707/47‑1 (MQ=255)
                                                                                                             aCTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATc                                                                                                                                                  >  1:1102707/1‑47 (MQ=255)
                                                                                                               tGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTc                                       <  1:678676/150‑1 (MQ=255)
                                                                                                                 cATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAg                                     <  1:51229/150‑1 (MQ=255)
                                                                                                                         ccGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGt                             >  1:548271/1‑150 (MQ=255)
                                                                                                                            gTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGGCGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCt                          <  1:237334/150‑1 (MQ=255)
                                                                                                                                 ttGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGAtc                     >  1:1100203/1‑150 (MQ=255)
                                                                                                                                 ttGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCCCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCa                                      >  2:333265/1‑133 (MQ=255)
                                                                                                                                 ttGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCCCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCa                                      <  1:333265/133‑1 (MQ=255)
                                                                                                                                     tCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGATCTCGa                 >  1:1061931/1‑150 (MQ=255)
                                                                                                                                      cAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGCTCTCGAc                >  1:491296/1‑150 (MQ=255)
                                                                                                                                       aGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGTGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAAc                               <  2:830246/134‑1 (MQ=255)
                                                                                                                                       aGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGTGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAAc                               >  1:830246/1‑134 (MQ=255)
                                                                                                                                       aGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCAGCCTCGGCCTggcggg                                                                 >  1:353496/1‑100 (MQ=255)
                                                                                                                                       aGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCAGCCTCGGCCTggcggg                                                                 <  2:353496/100‑1 (MQ=255)
                                                                                                                                          aTTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGATCTCGACGTTg            <  2:1053036/150‑1 (MQ=255)
                                                                                                                                           tttGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTc                                                                                        <  2:726143/73‑1 (MQ=255)
                                                                                                                                           tttGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTc                                                                                        >  1:726143/1‑73 (MQ=255)
                                                                                                                                             tGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCAGGAAGCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGATCTCGACGTTGGcc         <  2:413462/150‑1 (MQ=255)
                                                                                                                                                 cgGCTTGCATCAGTTGC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGATCTCGACGTTGGCCACCt     <  2:360088/150‑1 (MQ=255)
                                                                                                                                                    cTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGATCTCGACGTTGGCCACCTTGa  <  2:1100203/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CACTATTTCTCCTTGTCCACCAACCCCTTCAGCAGATCTATCGGCAGCGGGAAGACGATGGTCGAGCTCTTGTCCCCGGCAATCGTCCCCAGTGTCTGCATATAGCGCAACTGCATCGCCCCCGGTTCCTTGCTCAGCATTTGCGCGGCTTGCATCAGTTTC‑‑TCCGACGCCTGCAGTTCACCCTCGGCATGAATCACCTTGGCCCGACGTTCGCGCTCGGCCTCGGCCTGGCGGGCGATGGCGCGTACCATCGACTCGTTCAGGTCAACGTGCTTGATCTCGACGTTGGCCACCTTGA  >  NC_002947/249099‑249396

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