Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,624,601 G→A 7.8% G65R (GGG→AGG)  PP_2296 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,624,6010GA 7.8% 145.0 / 2.4 51G65R (GGG→AGG) PP_2296hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base A (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

CACGACCTGGACCCCAAAATGCCGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTCGAATCAGAGGACAAAAATGCCCA  >  NC_002947/2624457‑2624746
                                                                                                                                                |                                                                                                                                                 
cACGACCTGGACCCCAAAATGCCGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCg                                                                                                                                               >  4:16341/1‑149 (MQ=255)
  cGAGCTTGACCCCAAAATGCCGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCgtg                                                                                                                                             <  8:132072/149‑1 (MQ=255)
      cTGGACCCCAAAATGCCGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGtggtgg                                                                                                                                         >  4:103448/1‑149 (MQ=255)
           ccccAAAATGCCGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGCGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCa                                                                                                                                    >  6:54155/1‑149 (MQ=255)
              cAAAATGCCGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACg                                                                                                                                  >  4:212417/1‑148 (MQ=255)
              cAAAATGCCGATTCCGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACgg                                                                                                                                 >  3:145940/1‑149 (MQ=255)
                      cGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGt                                                                                                                         <  2:173567/149‑1 (MQ=255)
                          tCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCGGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTg                                                                                                                      <  4:145940/148‑1 (MQ=255)
                                aGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTcc                                                                                                               <  6:209586/149‑1 (MQ=255)
                                aGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTcc                                                                                                               <  3:84480/149‑1 (MQ=255)
                                aGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTcc                                                                                                               <  3:59446/149‑1 (MQ=255)
                                        aTTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGc                                                                                                        >  5:29506/1‑148 (MQ=255)
                                         tttGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTc                                                                                                          <  3:212417/145‑1 (MQ=255)
                                                aCTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGa                                                                                               <  6:29506/149‑1 (MQ=255)
                                                aCTCTTGCACAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTTGTTGGGGTAGGTGTAGCTGTCGCATCCCGATACGCTGCGACCTATGGGCGTGTTGGGCCCACGGACGCTGGTTTGGGCTTCATGGTCGATATGGAAGa                                                                                               >  8:276555/1‑149 (MQ=255)
                                                  tctTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGAcc                                                                                             >  6:310849/1‑149 (MQ=255)
                                                    ttGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGAccc                                                                                            >  3:101969/1‑148 (MQ=255)
                                                       cccAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAAc                                                                                        >  3:48816/1‑149 (MQ=255)
                                                         cAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCgtggt                                                                                                                                           >  5:211316/1‑96 (MQ=255)
                                                         cAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCgtggt                                                                                                                                           <  6:211316/96‑1 (MQ=255)
                                                           gCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGTGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGAtt                                                                                    >  4:32023/1‑149 (MQ=255)
                                                            cAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTCGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTg                                                                                   <  5:310849/149‑1 (MQ=255)
                                                            cAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTg                                                                                   >  4:156919/1‑149 (MQ=255)
                                                              aatagagCCTGTGGGCATTTAAGTAATGGCCGTCTTGGTTGGTGTAGGTGGAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTTGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGc                                                                                 <  2:156028/146‑1 (MQ=255)
                                                              aaaaGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGc                                                                                 <  1:89249/149‑1 (MQ=255)
                                                                 agagCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTa                                                                              <  3:16341/149‑1 (MQ=255)
                                                                 agagCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTa                                                                              <  4:48816/149‑1 (MQ=255)
                                                                  gagCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTaa                                                                             <  8:95831/149‑1 (MQ=255)
                                                                    gCCTGTGGGCATTTCAGAAATGGCCGACTTTGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTaag                                                                            <  7:276555/148‑1 (MQ=255)
                                                                        gtgGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATAGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGt                                                                                                           >  2:12735/1‑113 (MQ=255)
                                                                        gtgGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATAGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGt                                                                                                           <  1:12735/113‑1 (MQ=255)
                                                                            gCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGt                                                                                                           <  5:238131/109‑1 (MQ=255)
                                                                            gCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGt                                                                                                           >  6:238131/1‑109 (MQ=255)
                                                                             cATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATAGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTgg                                                                  >  3:179065/1‑149 (MQ=255)
                                                                               tttCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGgtg                                                                <  6:309371/149‑1 (MQ=255)
                                                                                ttCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGCTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGgtgt                                                               >  1:76386/1‑149 (MQ=255)
                                                                                  cAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGtt                                                              <  2:137103/148‑1 (MQ=255)
                                                                                   aGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCa                                                            >  4:53705/1‑149 (MQ=255)
                                                                                                       gggTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATAGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTc                                        <  4:179065/149‑1 (MQ=255)
                                                                                                          tAGGTGTAGCTGTCGCATCCCGCTACTCTGCGAGCTATGGGCGTGGTGGTCGCACGTACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCAtt                                     <  4:18341/149‑1 (MQ=255)
                                                                                                             gtgtAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGtt                                                              >  2:132821/1‑121 (MQ=255)
                                                                                                             gtgtAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGtt                                                              <  1:132821/121‑1 (MQ=255)
                                                                                                                tAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCa                                       >  7:116084/1‑141 (MQ=255)
                                                                                                                tAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCa                                       <  8:116084/141‑1 (MQ=255)
                                                                                                                     gTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGtt                          <  3:156919/149‑1 (MQ=255)
                                                                                                                      tCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTc                         >  7:43890/1‑149 (MQ=255)
                                                                                                                          aTCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTCGAAt                     <  3:103448/149‑1 (MQ=255)
                                                                                                                                   cGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTCGAATCAGAGGACa            <  8:74207/149‑1 (MQ=255)
                                                                                                                                          aGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCTACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTCGAATCAGAGGACAAAAATGc     >  3:14566/1‑149 (MQ=255)
                                                                                                                                            cTATGGGCGTGGTGGGCGGACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTCGAATCAGAGGACAAAAATGccc   <  3:32023/149‑1 (MQ=255)
                                                                                                                                             tATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTCGAATCAGAGGACAAAAATGCCCa  >  7:227572/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CACGACCTGGACCCCAAAATGCCGATTCGGCAAGTCATCTATTTGCTCACTCTTGCCCAGCAAAAAGAGCCTGTGGGCATTTCAGAAATGGCCGACTTGGTTGGGGTAGGTGTAGCTGTCGCATCCCGCTACGCTGCGAGCTATGGGCGTGGTGGGCGCACGGACGCTGGTTTGGGCTTCCTGGTCGCTAAGGAAGACCCGAACGATTGGCGTAAGAAGCAACTGGTGTTCACTGACGCCGGGCGCAAATTCATTGAAAGCTTGTTCGAATCAGAGGACAAAAATGCCCA  >  NC_002947/2624457‑2624746

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