Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,173,551 G→A 8.5% R76R (CGC→CGT)  PP_3672 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,173,5510GA 8.5% 130.5 / 2.4 47R76R (CGC→CGT) PP_3672hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base A (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

GCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTGCACCGCCGACAGGTTGACC  >  NC_002947/4173404‑4173698
                                                                                                                                                   |                                                                                                                                                   
gcgGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTg                                                                                                                                                     <  1:354441/148‑1 (MQ=255)
  ggCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTgcgc                                                                                                                                                  <  4:231844/149‑1 (MQ=255)
           acaTGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTctcc                                                                                                                                         >  7:21189/1‑149 (MQ=255)
           acaTGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTACGCACGTTctcc                                                                                                                                         >  5:324572/1‑149 (MQ=255)
             aTGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTctcctc                                                                                                                                       <  1:196257/149‑1 (MQ=255)
                 cATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGtt                                                                                                                                             <  8:281184/139‑1 (MQ=255)
                 cATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGtt                                                                                                                                             >  7:281184/1‑139 (MQ=255)
                      acgaacgaCTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGcc                                                                                                                              >  2:260464/1‑149 (MQ=255)
                        gaacgaCTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGGAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGCAATGGACACGCCCAGGTCTTTCAGCGCGAGGAAGGTACGCACGTTGTCCTCGCTGTTGcc                                                                                                                              >  8:30537/1‑147 (MQ=255)
                          acgaCTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCcagc                                                                                                                          <  3:110176/149‑1 (MQ=255)
                           cgaCTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCcagca                                                                                                                         <  3:217499/149‑1 (MQ=255)
                                 gTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTAcc                                                                                                                    >  7:246768/1‑148 (MQ=255)
                                    gATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTgc                                                                                                                                                    <  8:207464/113‑1 (MQ=255)
                                    gATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTgc                                                                                                                                                    >  7:207464/1‑113 (MQ=255)
                                    gATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCt                                                                                                                                    <  3:194477/129‑1 (MQ=255)
                                    gATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCt                                                                                                                                    >  4:194477/1‑129 (MQ=255)
                                        ttGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCgg                                                                                                             >  7:22150/1‑148 (MQ=255)
                                          gATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTaa                                                                                                          >  7:84188/1‑149 (MQ=255)
                                          gATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTaa                                                                                                          >  7:208011/1‑149 (MQ=255)
                                            tCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTACGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAAtt                                                                                                        <  6:324572/149‑1 (MQ=255)
                                               ggTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTtc                                                                                                      >  5:55552/1‑148 (MQ=255)
                                                  cGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTtccagc                                                                                                  >  5:184613/1‑149 (MQ=255)
                                                             gaacgCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTg                                                                                       >  6:354869/1‑149 (MQ=255)
                                                               acgCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTgcg                                                                                     >  5:73540/1‑149 (MQ=255)
                                                                cgCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTgcgg                                                                                    <  6:97901/149‑1 (MQ=255)
                                                                   aGGTAGCCCAACGATTATTAGCCGGTACCAAAGTCATCCAGTCAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGcggcgg                                                                                 <  2:279210/149‑1 (MQ=255)
                                                                   aGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGcggcgg                                                                                 >  3:125542/1‑149 (MQ=255)
                                                                        gCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAgg                                                                             <  8:21189/148‑1 (MQ=255)
                                                                                        ccGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAgcgc                                                            >  1:359075/1‑149 (MQ=255)
                                                                                           gTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAg                                                         <  8:208011/149‑1 (MQ=255)
                                                                                                  aGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGGACGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACg                                                  <  7:30537/149‑1 (MQ=255)
                                                                                                       tCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATg                                             >  8:111227/1‑149 (MQ=255)
                                                                                                         cAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGc                                            <  8:261560/148‑1 (MQ=255)
                                                                                                         cAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGc                                            >  7:261560/1‑148 (MQ=255)
                                                                                                         cAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCa                                           >  2:206981/1‑149 (MQ=255)
                                                                                                          aGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCAc                                          <  6:184613/149‑1 (MQ=255)
                                                                                                            tGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCTGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCa                                        <  1:260464/149‑1 (MQ=255)
                                                                                                             gAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCa                                        <  8:246768/148‑1 (MQ=255)
                                                                                                                 tGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCg                                   <  4:125542/149‑1 (MQ=255)
                                                                                                                   gACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGcc                                 <  2:359075/149‑1 (MQ=255)
                                                                                                                    acacGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCg                                >  3:178542/1‑149 (MQ=255)
                                                                                                                          ccAGGGCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTg                          <  8:310743/149‑1 (MQ=255)
                                                                                                                            aGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGaa                        >  2:248117/1‑149 (MQ=255)
                                                                                                                             ggTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAAc                       >  8:180450/1‑149 (MQ=255)
                                                                                                                                cTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTg                     >  4:111024/1‑148 (MQ=255)
                                                                                                                                cTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTGc                    >  8:154179/1‑149 (MQ=255)
                                                                                                                                  ttCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTGCAc                  >  7:96261/1‑149 (MQ=255)
                                                                                                                                     aGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGGACCAGATCGCCGTTCTTGAACTGCAccg                <  1:206981/148‑1 (MQ=255)
                                                                                                                                     aGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTGCAccg                <  8:22150/148‑1 (MQ=255)
                                                                                                                                         cgcgcAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTGCACCGCCGAc           <  1:248117/149‑1 (MQ=255)
                                                                                                                                                  tGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTGCACCGCCGACAGGTTGAcc  <  5:354869/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGGCTGTCGCACATGTCATGCACGAACGACTTGTCGATCTTGATCCGGTCGAACGGGAACGAACGCAGGTAGCCCAACGATGAGTAGCCGGTACCAAAGTCATCCAGTGAAATGGACACGCCCAGGTCTTTCAGCGCGCGCAAGGTGCGCACGTTCTCCTCGCTGTTGCCCAGCAGTACCGACTCGGTAATTTCCAGCTCCAGCCGCTGCGGCGGCAGGCCCGAGTCGGCCAGCGCGAGGGCGACGGTATGCACCAGATCGCCGTTCTTGAACTGCACCGCCGACAGGTTGACC  >  NC_002947/4173404‑4173698

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