Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,357,497 T→A 10.4% V289E (GTG→GAG) ‡ PP_3831 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,357,4970TA10.4% 116.0 / 4.5 47V289E (GTG→GAG) ‡PP_3831hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/22);  new base A (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

CCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAG  >  NC_002947/4357353‑4357642
                                                                                                                                                |                                                                                                                                                 
ccAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGc                                                                                                                                               <  8:189597/149‑1 (MQ=255)
ccAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGc                                                                                                                                               >  5:63666/1‑149 (MQ=255)
 cAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGcc                                                                                                                                              <  2:216584/149‑1 (MQ=255)
  aaGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGACACGCGCCTGGGCAATACGCCGGCAGTGTGCCg                                                                                                                                             >  8:284213/1‑149 (MQ=255)
    ggACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgc                                                                                                                                           <  3:299580/149‑1 (MQ=255)
          tCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCgctgct                                                                                                                                     <  2:6178/149‑1 (MQ=255)
                              gcCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTg                                                                                                                 >  3:269739/1‑149 (MQ=255)
                                cTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTg                                                                                                                 >  2:304891/1‑147 (MQ=255)
                                     aTTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgcg                                                                                                                                          >  6:121715/1‑117 (MQ=255)
                                     aTTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCACACTGCCGgcg                                                                                                                                          <  5:121715/117‑1 (MQ=255)
                                        gAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTAt                                                                                                        <  2:415/148‑1 (MQ=255)
                                           ccTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCt                                                                                                    >  1:141767/1‑149 (MQ=255)
                                              gCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTCGAGCACTATGCCTAgg                                                                                                 >  2:101681/1‑149 (MQ=255)
                                                tAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTgggg                                                                                               >  2:107282/1‑149 (MQ=255)
                                                tAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTgggg                                                                                               >  3:152162/1‑149 (MQ=255)
                                                   aGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGgc                                                                                              <  1:137882/147‑1 (MQ=255)
                                                   aGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGgc                                                                                              >  2:137882/1‑147 (MQ=255)
                                                       cTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTgg                                                                                        <  1:107282/149‑1 (MQ=255)
                                                       cTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTgg                                                                                        <  5:268629/149‑1 (MQ=255)
                                                        ttGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGc                                                                                       >  3:28596/1‑149 (MQ=255)
                                                            ccTGGGAGCCGGAGAGTGAAGCCAATCGCCTGGTCGCCGGAATAGTCCGACAGGTAGCCACGCGCCTGGGCAATACTCGGGGAGTGTGCCGGCGCTGCTAAATCCTCCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAt                                                                                   <  1:219116/149‑1 (MQ=255)
                                                                ggAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAAttt                                                                                 <  1:304891/147‑1 (MQ=255)
                                                                    ccGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGa                                                                            >  3:156430/1‑148 (MQ=255)
                                                                       gagaGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTCGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCCGCAGTTCAGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGaaaaa                                                                        >  3:226610/1‑149 (MQ=255)
                                                                            tGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCgt                                                                    <  6:63666/148‑1 (MQ=255)
                                                                               aGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTgc                                                                >  4:163263/1‑149 (MQ=255)
                                                                                     aCGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAg                                                          <  4:152162/149‑1 (MQ=255)
                                                                                         cAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCt                                                      <  7:240250/149‑1 (MQ=255)
                                                                                         cAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCt                                                      >  7:107487/1‑149 (MQ=255)
                                                                                          aGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTg                                                     >  3:212258/1‑149 (MQ=255)
                                                                                            gTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGa                                                   >  6:93537/1‑149 (MQ=255)
                                                                                             tCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCgtgt                                                                  <  7:171023/133‑1 (MQ=255)
                                                                                             tCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCgtgt                                                                  >  8:171023/1‑133 (MQ=255)
                                                                                                  gcAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCgga                                             <  4:28596/149‑1 (MQ=255)
                                                                                                           cGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTgg                                     <  7:284213/148‑1 (MQ=255)
                                                                                                                 gTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCt                              <  2:141767/149‑1 (MQ=255)
                                                                                                                    gCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCtgtt                           <  8:107487/149‑1 (MQ=255)
                                                                                                                    gCCACGCGCCCGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGCGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCtgt                            <  4:269739/148‑1 (MQ=255)
                                                                                                                      cACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCg                                                                             <  1:284762/97‑1 (MQ=255)
                                                                                                                      cACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCg                                                                             >  2:284762/1‑97 (MQ=255)
                                                                                                                            ccTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTttcttc                   <  4:175327/149‑1 (MQ=255)
                                                                                                                                  cAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCATTTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCAc             <  3:163263/149‑1 (MQ=255)
                                                                                                                                  cAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCAc             >  2:273484/1‑149 (MQ=255)
                                                                                                                                     tACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTATTAAGGCACtt           >  7:162783/1‑148 (MQ=255)
                                                                                                                                     agcGCCGGCACCCTGCCGGCGGTGCTACCTCCACCCGGCAGTTCTGGAGCACTCTGGCTTGGGGCGCCTGGCGAATTTGCCGAATAACCGTGTGCGCAACGGCCTGGAGCCGGAGGCAGTGGCCTGGCGGGTGTTTCTTCCGGCACGTg          >  8:300261/3‑149 (MQ=255)
                                                                                                                                     agcGCCGGCACACTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTg          >  7:37093/3‑149 (MQ=255)
                                                                                                                                             cAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACag  <  5:138132/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCAAGGACTATCGCACCTGGGCCGGCAGCAGCCTGGCATTGAACCTGCTAAAGCCCTTGGCCTGGGAGCCGGAGAGTGAAGCCAAACGCCAGGTCGCCGCAATAGTCCGCCAGGTAGCCACGCGCCTGGGCAATACGCCGGCAGTGTGCCGGCGCTGCTACATCCACCCGGCAGTTCTGGAGCACTATGCCTTGGGGCGCCTGGCCAATTTGCCGAAAAACCGTGTGCGCAAAGGCCTGGACCCGGAGGAAGTGGCCTTGCTGTTGTTTCTTCAGGCACTTGAGGAACAG  >  NC_002947/4357353‑4357642

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