Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,928,936 G→A 7.2% R386* (CGA→TGA)  PP_1728 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,928,9360GA 7.2% 159.0 / 2.3 55R386* (CGA→TGA) PP_1728membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (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 = 9.89e-01

ATCATCGACACGCCCCACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAAC  >  NC_002947/1928788‑1929075
                                                                                                                                                    |                                                                                                                                           
atcatcGACACGCCCCACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCg                                                                                                                                             >  6:47821/1‑149 (MQ=255)
   atcGACACGCCCCACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCa                                                                                                                                           <  2:214584/148‑1 (MQ=255)
             cccACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGtt                                                                                                                                 >  3:222798/1‑148 (MQ=255)
             cccACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTg                                                                                                                                <  2:205910/149‑1 (MQ=255)
             cccACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCACACCAGGCCGtt                                                                                                                                 >  7:214703/1‑148 (MQ=255)
             cccACACGCGCAGGCCCAACTCGATGAAAACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTCCGCAGGCGTTCGGTGTACAACGCCTGGCGCTTGCTGGCCCGCTAGGGGTTGGCGTCATGGCGTCACACCAGGCCGtt                                                                                                                                 >  5:13804/1‑148 (MQ=255)
              ccACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGtt                                                                                                                                 >  5:39747/1‑147 (MQ=255)
              ccACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGa                                                                                                                               <  6:119359/149‑1 (MQ=255)
                acacGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGc                                                                                                                             >  5:35984/1‑149 (MQ=255)
                    gcgcAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGc                                                                                                                                     >  3:122259/1‑137 (MQ=255)
                    gcgcAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGc                                                                                                                                     <  4:122259/137‑1 (MQ=255)
                          gCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGcaca                                                                                                                   <  2:106339/149‑1 (MQ=255)
                           cccAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCacac                                                                                                                  >  4:165158/1‑149 (MQ=255)
                              aaCTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCAcaccac                                                                                                               >  3:153751/1‑149 (MQ=255)
                              aaCTCGATGAACACCAGCCAGAGAAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAAGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCAaaccac                                                                                                               >  7:53967/1‑149 (MQ=255)
                                   gATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACAccaccacca                                                                                                          <  6:39747/149‑1 (MQ=255)
                                      gAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGtt                                                                                                                                 >  3:186838/1‑123 (MQ=255)
                                      gAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGtt                                                                                                                                 <  4:186838/123‑1 (MQ=255)
                                       aaCACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCAc                                                                                                      >  7:310955/1‑149 (MQ=255)
                                         cacCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACgg                                                                                                    >  5:30123/1‑149 (MQ=255)
                                              gccagATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACgg                                                                                                    >  1:166142/1‑144 (MQ=255)
                                              gccagATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACgg                                                                                                    <  2:166142/144‑1 (MQ=255)
                                              gccagATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGca                                                                                                >  3:158077/1‑148 (MQ=255)
                                                   aTGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAg                                                                                          <  6:274953/149‑1 (MQ=255)
                                                             gCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGgcg                                                                                 >  1:75851/1‑148 (MQ=255)
                                                              ccccAGCAGGTAGCCGAAGTTCCGCACGAGTTCCGTCTAAACCGCCTCGCGCTTGCTGACCCGCTTGGGGTTGGCGGCATCGCGTCGCACCAGGCCGTAGAGCACCATCCACACCACCACCAGCACGGTGCACATCAGCGACTGgcg                                                                                 >  7:22896/3‑147 (MQ=255)
                                                                ccAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAgg                                                                             <  6:270987/149‑1 (MQ=255)
                                                                  agcagGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGc                                                                           >  5:50798/1‑149 (MQ=255)
                                                                  agcagGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGc                                                                           >  5:56097/1‑149 (MQ=255)
                                                                       gTAGCCCAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGATTGCTGGCCCGCTTGGGGTTGGCGGAATGGCCTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCATAACGGTGCACATCAGCGACTGGCGAAGGGctgagc                                                                      >  8:233359/1‑144 (MQ=255)
                                                                          gCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTg                                                                    <  4:158077/148‑1 (MQ=255)
                                                                          gCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTg                                                                    <  6:56097/148‑1 (MQ=255)
                                                                          gCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTg                                                                    <  4:153751/148‑1 (MQ=255)
                                                                          gCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTg                                                                    <  6:50798/148‑1 (MQ=255)
                                                                           ccGCAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGGTGGCTGGCCGGCTTGGTGTTGGCGGCATGGCGGAGCACGTGGCGGTTGAGCACCATCCACACCACCACCAGCTCGTTGCACATGAGCGTGTGGCGCAGGTGGGTTCTGgt                                                                  <  8:22896/149‑1 (MQ=255)
                                                                           ccGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACTGTGCACAGCACCGACTGGCGCAGGGCGGGGCTGgt                                                                  <  6:35984/149‑1 (MQ=255)
                                                                              aaGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTc                                                               <  6:30123/149‑1 (MQ=255)
                                                                                gTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGc                                                             <  8:310955/149‑1 (MQ=255)
                                                                                       aGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCGCGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGga                                                      <  5:47821/149‑1 (MQ=255)
                                                                                         gCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCTGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGgaga                                                    <  3:181353/149‑1 (MQ=255)
                                                                                            ttCGGTGTACACCTCATGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCACAGCTGTCGGTTGATCTCCTTGCACACCACCAGGATCACGTTGCACATCTGCGACTTGCGCAGGGCGGTGCTGGTGTCGCCGGGGGAGActa                                                 <  6:13804/149‑3 (MQ=255)
                                                                                            ttCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCACACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGa                                                 <  8:214703/149‑1 (MQ=255)
                                                                                                   tACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGTGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGgaga                                                    >  3:288289/1‑139 (MQ=255)
                                                                                                   tACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGTGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGgaga                                                    <  4:288289/139‑1 (MQ=255)
                                                                                                    acacCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCg                                         >  5:300607/1‑149 (MQ=255)
                                                                                                      acCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGaa                                       <  6:300607/149‑1 (MQ=255)
                                                                                                          ccTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGc                                   >  6:88084/1‑149 (MQ=255)
                                                                                                           cTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGc                                   >  5:209428/1‑148 (MQ=255)
                                                                                                             ggCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGCCGAAGGTGGCGATCAGCgca                                >  5:208442/1‑147 (MQ=255)
                                                                                                                gcTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATg                             <  3:165158/149‑1 (MQ=255)
                                                                                                                      tGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCTCCATGCCGTTGAGCACCCTGCACACCACCACCAGCACGGTGCCCAGCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGGACAGCGAAATGGccgc                        >  5:367835/1‑146 (MQ=255)
                                                                                                                        gCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGccac                        >  6:234661/1‑146 (MQ=255)
                                                                                                                        gCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGccac                        >  6:300809/1‑146 (MQ=255)
                                                                                                                           cGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCTTAGAGCACCATGCACACCACCACCTGCACTGTTCACATCAGCGACTGGCTCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAAt                  <  7:233359/149‑1 (MQ=255)
                                                                                                                                ggggTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAg               <  6:209428/147‑1 (MQ=255)
                                                                                                                                            aTGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAAc  >  6:166419/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
ATCATCGACACGCCCCACACGCGCAGGCCCAACTCGATGAACACCAGCCAGATGAACAGGTGCACCAGCAGGTAGCCGAAGTTGCGCAGGCGTTCGGTGTACACCGCCTGGCGCTTGCTGGCCCGCTTGGGGTTGGCGGCATGGCGTCGCACCAGGCCGTTGAGCACCATGCACACCACCACCAGCACGGTGCACATCAGCGACTGGCGCAGGGCGGTGCTGGTGTCGCCGGCGGAGACGAAGGTGGCGAACAGCGAAATGGCCACCAGAATCAGCGCCGGAATGAAC  >  NC_002947/1928788‑1929075

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