Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,388,807 C→T 12.0% R291H (CGC→CAC) ‡ PP_2991 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,388,8070CT12.0% 126.5 / 7.4 50R291H (CGC→CAC) ‡PP_2991hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (3/3);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

GGCCGGCGAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAG  >  NC_002947/3388664‑3388941
                                                                                                                                               |                                                                                                                                      
ggCCGGCGAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTc                                                                                                                                   >  3:44530/1‑149 (MQ=255)
ggCCGGCGAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTc                                                                                                                                   >  3:48223/1‑149 (MQ=255)
ggCCGGCGAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTc                                                                                                                                   <  1:99298/149‑1 (MQ=255)
  ccGGCGAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGa                                                                                                                                 <  4:108868/149‑1 (MQ=255)
   cGGCGAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAg                                                                                                                                <  8:241915/149‑1 (MQ=255)
          tGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTcgcg                                                                                                                                       >  3:175616/1‑135 (MQ=255)
          tGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTcgcg                                                                                                                                       <  4:175616/135‑1 (MQ=255)
             tttGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACg                                                                                                                        >  4:230626/1‑147 (MQ=255)
             tttGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCa                                                                                                                      >  5:90907/1‑149 (MQ=255)
                  gTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGc                                                                                                                       <  6:11334/143‑1 (MQ=255)
                  gTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGc                                                                                                                       >  5:11334/1‑143 (MQ=255)
                   tGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGc                                                                                                                >  8:114261/1‑149 (MQ=255)
                   tGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCAga                                                                                                                    <  2:148492/145‑2 (MQ=255)
                   tGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCAga                                                                                                                    >  1:148492/1‑144 (MQ=255)
                         gagCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGaaa                                                                                                          <  6:90907/149‑1 (MQ=255)
                               aCTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCa                                                                                                     >  4:9121/1‑148 (MQ=255)
                                   tGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAAcc                                                                                                >  2:130506/1‑149 (MQ=255)
                                   tGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCAga                                                                                                                    <  2:240451/129‑2 (MQ=255)
                                   tGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCATGGCTCGAGCCGCGACGCAga                                                                                                                    >  1:240451/1‑128 (MQ=255)
                                     aaaCGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCAt                                                                                                    <  1:135051/143‑1 (MQ=255)
                                     aaaCGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCAt                                                                                                    >  2:135051/1‑143 (MQ=255)
                                         gCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCg                                                                                                               >  2:186550/1‑128 (MQ=255)
                                              catgTATTCGTACTTCATCATTTTCCCCAACTTCATTGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGATGGTGTCTGAATCTGCGTCGCGGCTCGAga                                                                                                                               <  8:197799/104‑2 (MQ=255)
                                              caagTATTAGTACTTCATAATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTAGCGGCTAGAGc                                                                                                                               >  7:197799/4‑107 (MQ=255)
                                              cACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACt                                                                                     <  1:130506/149‑1 (MQ=255)
                                                 gTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGc                                                                                  <  4:48223/149‑1 (MQ=255)
                                                 gTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGc                                                                                  <  4:44530/149‑1 (MQ=255)
                                                  tATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGGGGGTGAATCTGCGTCGCGGCTCTAGCCATGACGCATTTGGCTGGAAACCGCATAACCGTATTCCAACTGGcc                                                                                 <  2:43448/149‑1 (MQ=255)
                                                   aTTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCg                                                                                                       <  3:27476/126‑1 (MQ=255)
                                                   aTTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCg                                                                                                       >  4:27476/1‑126 (MQ=255)
                                                       gTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCtga                                                                             <  2:118767/148‑1 (MQ=255)
                                                          cTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTgatgat                                                                         >  3:261687/1‑149 (MQ=255)
                                                              atcatTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGTCCTGATGATc                                                                        >  1:228041/1‑146 (MQ=255)
                                                              atcatTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATc                                                                        <  2:228041/146‑1 (MQ=255)
                                                                       cccAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGc                                                                                          <  5:1566/119‑1 (MQ=255)
                                                                       cccAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGc                                                                                          >  6:1566/1‑119 (MQ=255)
                                                                        ccAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAg                                                           <  5:25105/149‑1 (MQ=255)
                                                                        ccAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAg                                                           >  6:25105/1‑149 (MQ=255)
                                                                                       ggAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGCGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGat                                             >  6:170285/1‑148 (MQ=255)
                                                                                       ggAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTgg                                                                <  4:14692/129‑1 (MQ=255)
                                                                                       ggAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTgg                                                                >  3:14692/1‑129 (MQ=255)
                                                                                               ttCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGGTGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTc                                    >  5:231944/1‑149 (MQ=255)
                                                                                                   aCAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGTga                                <  4:261687/149‑1 (MQ=255)
                                                                                                             gTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGt                      >  8:60043/1‑149 (MQ=255)
                                                                                                             gTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGt                      <  3:230626/149‑1 (MQ=255)
                                                                                                                 ccTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCga                                <  4:173140/135‑1 (MQ=255)
                                                                                                                 ccTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCga                                >  3:173140/1‑135 (MQ=255)
                                                                                                                        tgAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCttcttc           <  7:60043/149‑1 (MQ=255)
                                                                                                                            ggTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCtt       <  7:114261/149‑1 (MQ=255)
                                                                                                                               ggCTGAAGCTGCGGCGCGGCTCTAGCCATGACGCAGTTGGCGGGAACCCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGa    <  4:65541/149‑1 (MQ=255)
                                                                                                                                gCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGaa   <  3:9121/149‑1 (MQ=255)
                                                                                                                                  tGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAg  >  3:131635/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                      
GGCCGGCGAGTGGTTTGGGTGCATAGAGCGCACTCTGAAACGCTGGCACGTATTCGTACTTCATCATTTTCCCCAACTTCAGGGTGCGGAGGTAATTCGACAGAGCGCCGTCGCCTAGTGTGAGGGTGGCTGAATCTGCGTCGCGGCTCGAGCCATGACGCAGTTGGCGGGAAACCGCATAACCGTATGCCAACTGGCCCTGATGATCGAGATTGGTAAAGCTGCTGGTGACGATATCGACTTCGCGAGAGAGCTTGTCGAGCTTCTTCGCTTTGAAG  >  NC_002947/3388664‑3388941

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