Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,938,882 C→T 7.7% E73K (GAG→AAG)  PP_1738 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,938,8820CT 7.7% 141.8 / 3.0 52E73K (GAG→AAG) PP_1738hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base T (2/2);  total (25/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 = 4.40e-01

ACAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTC  >  NC_002947/1938735‑1939023
                                                                                                                                                   |                                                                                                                                              
aCAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCGACGTGAACTcgc                                                                                                                                              >  2:1036985/1‑150 (MQ=255)
  aGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCgcgc                                                                                                                                            <  1:1036985/150‑1 (MQ=255)
                   aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAg                                                                                                                          >  1:613396/1‑150 (MQ=255)
                   aTTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATAAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAg                                                                                                                          >  1:615297/1‑150 (MQ=255)
                    ttGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTAGGGCAACGTGAACTTGCGCCCGGA‑GCGCGAGTTCac                                                                                                                          >  1:874951/1‑148 (MQ=255)
                    ttGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTAGGGCAACGTGAACTTGCGC‑CCGGAGCGCGAGTTCac                                                                                                                          <  2:874951/149‑2 (MQ=255)
                      gATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCt                                                                                                                       >  1:983513/1‑150 (MQ=255)
                        tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGt                                                                                                                     <  1:1112470/150‑1 (MQ=255)
                        tgctgctgcTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGt                                                                                                                     >  2:1094535/1‑150 (MQ=255)
                          ctgctgctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAg                                                                                                                   <  1:980403/150‑1 (MQ=255)
                             ctgctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCg                                                                                                                >  1:1021316/1‑150 (MQ=255)
                               gctgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGa                                                                                                              <  1:212397/150‑1 (MQ=255)
                                 tgTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCACAGAAGTTCTCGCGACACGTTCATTCGCGCCACGC‑GGTCAATTGCAGAGGGTAGTCGGACt                                                                                                            >  2:140405/1‑150 (MQ=255)
                                      tcatcgcCGGCACCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGc                                                                                                                        <  2:498336/129‑1 (MQ=255)
                                      ggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGC‑CCGGAGCGCGAGTTCacg                                                                                                                         <  1:535276/132‑3 (MQ=255)
                                      ggCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTTGCGC‑CCGGAGCGCGAGTTCacg                                                                                                                         >  2:535276/1‑130 (MQ=255)
                                           gcCGGTGCCGGCGATGAAGAAGATCAGAGACCCTTCGTTATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTc                                                                                                  <  2:1104719/150‑1 (MQ=255)
                                              ggTGCCGGCGATGATTAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCATAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTTTAGTCGGACTGAGTGATTTCAtc                                                                                               <  2:562803/150‑1 (MQ=255)
                                                           atgaAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGc                                                                                  >  1:723254/1‑150 (MQ=255)
                                                             gaAGATCAGCGGCGCTTCTTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTc                                                                                >  2:1022496/1‑150 (MQ=255)
                                                                gATCAGCGGCGCTTCGTGATCCTGGCGTGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAgg                                                                             <  2:965765/150‑1 (MQ=255)
                                                                           cTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCcg                                                                  <  2:1021316/150‑1 (MQ=255)
                                                                                 gATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTACGgcggcg                                                            <  2:983513/150‑1 (MQ=255)
                                                                                   tCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtg                                                          >  2:486659/1‑150 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtc                                                                                               <  1:766135/110‑1 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCAtc                                                                                               >  2:766135/1‑110 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTc                                                                    <  2:1040129/137‑1 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTc                                                                    >  1:1040129/1‑137 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggtg                                                       <  2:615297/150‑1 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggtg                                                       <  2:613396/150‑1 (MQ=255)
                                                                                      tGGCGGGCCAGGCGGTACTTGAGCTAGGTGACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGACTGTAGTCTGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCgtggtg                                                       >  1:174684/1‑150 (MQ=255)
                                                                                                 gCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAt                                                  <  2:498551/144‑1 (MQ=255)
                                                                                                 gCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAt                                                  >  1:498551/1‑144 (MQ=255)
                                                                                                        ttGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAAc                                     >  1:218381/1‑150 (MQ=255)
                                                                                                           aGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACggg                                  <  1:786080/150‑1 (MQ=255)
                                                                                                           aGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACggg                                  <  2:551956/150‑1 (MQ=255)
                                                                                                             cttcttCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGtt                                >  1:673280/1‑150 (MQ=255)
                                                                                                                cttcACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGt                             >  2:427231/1‑150 (MQ=255)
                                                                                                                    aCTGCCGAGAAGTTGTCGGGCAACGTGAACTCGCTCTGCGG‑GCGCAAGTTCATGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCTCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAaga                         <  2:178860/150‑1 (MQ=255)
                                                                                                                        tccAGGAGTTGGCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTAAGGCTATAGTCGCACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGt                     >  2:343313/2‑150 (MQ=255)
                                                                                                                         ccAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGc                                                                                     <  1:948034/85‑1 (MQ=255)
                                                                                                                         ccAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGc                                                                                     >  2:948034/1‑85 (MQ=255)
                                                                                                                         ccAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTa                    <  2:723254/150‑1 (MQ=255)
                                                                                                                                     gggCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGc        <  1:427231/150‑1 (MQ=255)
                                                                                                                                       gCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGCTTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg      >  2:15236/1‑150 (MQ=255)
                                                                                                                                       gCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAg      >  2:335329/1‑150 (MQ=255)
                                                                                                                                         aaCGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        <  1:673120/130‑1 (MQ=255)
                                                                                                                                         aaCGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        >  2:673120/1‑130 (MQ=255)
                                                                                                                                           cGTGACCTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTc  <  1:1094535/150‑1 (MQ=255)
                                                                                                                                              gAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAtggtg                                              >  1:320999/1‑103 (MQ=255)
                                                                                                                                              gAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGAtggtg                                              <  2:320999/103‑1 (MQ=255)
                                                                                                                                                  tCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        <  2:567485/121‑1 (MQ=255)
                                                                                                                                                  tCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAagag                        >  1:567485/1‑121 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
ACAGCTTTTTCAGGTATTCATTGATGCTGCTGCTGTAGGGCGCGCCGGTGCCGGCGATGATGAAGATCAGCGGCGCTTCGTGATCCTGGCGGGCCAGGCGGTACTTGAGCTTCTTCACTGCCCAGAAGTTGTCGGGCAACGTGAACTCGCGCTCCGG‑GCGCAAGTTCAGGCTGTAGTCGGACTGAGTGATTTCATCGTCGCTGGGCAGCTCAGGGCGCTTGTCCGGCGGCGTGGTGGCGATGGTGGCCTCGAACGGGTTGGTAAGAGGGTAGCCATAGCTGGCAGCGTC  >  NC_002947/1938735‑1939023

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