Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,125,675 C→T 9.4% E53K (GAA→AAA) ‡ PP_3629 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,125,6750CT 9.4% 147.1 / 4.7 53E53K (GAA→AAA) ‡PP_3629hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/20);  new base T (2/3);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

TCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTT  >  NC_002947/4125530‑4125819
                                                                                                                                                 |                                                                                                                                                
tCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCAc                                                                                                                                               >  3:216296/1‑149 (MQ=255)
tCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCAc                                                                                                                                               >  2:381104/1‑149 (MQ=255)
      tcgtTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGGAATTCACGGTAAt                                                                                                                                         <  8:119720/149‑1 (MQ=255)
        gtTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGGAATTCACGGTAAtgg                                                                                                                                       <  7:333327/149‑3 (MQ=255)
         ttACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATcac                                                                                                                                      <  4:216296/149‑1 (MQ=255)
         ttACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATca                                                                                                                                       >  1:198383/1‑148 (MQ=255)
                cAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGccc                                                                                                                               <  3:28549/149‑1 (MQ=255)
                   tGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGc                                                                                                                            >  5:25361/1‑149 (MQ=255)
                           cacaAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGGAATTCACGGTAAt                                                                                                                                         >  1:328528/1‑128 (MQ=255)
                           cacaAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGGAATTCACGGTAAt                                                                                                                                         <  2:328528/128‑1 (MQ=255)
                           cacaAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACg                                                                                                                                              <  8:135521/123‑1 (MQ=255)
                                        tctcTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATcggcg                                                                                                       <  6:212356/149‑1 (MQ=255)
                                            tGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGGCTGAAGGGTAATCGGCGGGAg                                                                                                   >  8:49397/1‑149 (MQ=255)
                                              ccTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGcc                                                                                                 >  1:378739/1‑149 (MQ=255)
                                                     cgcgGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGca                                                                                          >  1:142659/1‑149 (MQ=255)
                                                           aCTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGa                                                                                    >  1:11853/1‑149 (MQ=255)
                                                               aCCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAAc                                                                                >  7:190826/1‑149 (MQ=255)
                                                                   ggAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCgg                                                                            >  4:333269/1‑149 (MQ=255)
                                                                      aCGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTa                                                                         >  3:247448/1‑149 (MQ=255)
                                                                         aaaaTGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCg                                                                      >  7:9356/1‑149 (MQ=255)
                                                                               ccTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGaaa                                                                <  8:190826/149‑1 (MQ=255)
                                                                                cTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAg                                                               >  8:138226/1‑149 (MQ=255)
                                                                                 tCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGt                                                              >  6:137096/1‑149 (MQ=255)
                                                                                   acaGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGtt                                                             <  3:333269/148‑1 (MQ=255)
                                                                                    caGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCg                                                           >  6:6276/1‑149 (MQ=255)
                                                                                     aGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCgg                                                          <  2:198383/149‑1 (MQ=255)
                                                                                          ttGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCg                                                       >  6:192014/1‑147 (MQ=255)
                                                                                             aaGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGaa                                                  <  5:6276/149‑1 (MQ=255)
                                                                                                      ccaccaTTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTcgccgc                                         >  3:129393/1‑149 (MQ=255)
                                                                                                       caccaTTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCgccgcc                                        >  8:225860/1‑149 (MQ=255)
                                                                                                          caTTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCGTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCctct                                     >  1:181539/1‑149 (MQ=255)
                                                                                                                   ccACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCg                            <  7:49397/149‑1 (MQ=255)
                                                                                                                     aCGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTg                          <  8:9356/149‑1 (MQ=255)
                                                                                                                     aCGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTg                          >  6:137124/1‑149 (MQ=255)
                                                                                                                     aCGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATAGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTg                          >  7:202884/1‑149 (MQ=255)
                                                                                                                      cGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTgg                         <  7:138226/149‑1 (MQ=255)
                                                                                                                       gCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGaaa                                                                <  3:44236/109‑1 (MQ=255)
                                                                                                                       gCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGaaa                                                                >  4:44236/1‑109 (MQ=255)
                                                                                                                         aaGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCg                                  >  8:239459/1‑137 (MQ=255)
                                                                                                                         aaGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCg                                  <  7:239459/137‑1 (MQ=255)
                                                                                                                          aGGGATACCCATTACCGTGAATTCCAGGGTGATCACATCGCCCTGCTTGCCTGAATGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGt                     <  2:378739/149‑1 (MQ=255)
                                                                                                                             gATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTAcc                  <  1:381104/149‑1 (MQ=255)
                                                                                                                                aCCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTc                       >  6:122552/1‑141 (MQ=255)
                                                                                                                                aCCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTc                       <  5:122552/141‑1 (MQ=255)
                                                                                                                                    atcaCCGTGGAATTCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAgg             >  2:231954/4‑147 (MQ=255)
                                                                                                                                    aTTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAgg             >  4:302046/1‑147 (MQ=255)
                                                                                                                                        ccGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGAt        <  2:11853/148‑1 (MQ=255)
                                                                                                                                        ccGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGAt        <  5:192014/148‑1 (MQ=255)
                                                                                                                                        ccGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATg       >  6:52365/1‑149 (MQ=255)
                                                                                                                                         cGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATgg      >  2:300871/1‑149 (MQ=255)
                                                                                                                                             aaTTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGtt  <  5:100505/149‑1 (MQ=255)
                                                                                                                                             aaTTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGATAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGc            <  4:193748/139‑1 (MQ=255)
                                                                                                                                             aaTTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGATAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGc            >  3:193748/1‑139 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATTCCACGGTGATCACATCGCCCTGCTTGCCTGAAGGGTAATCGGCGGGAGCCTGGTGCACAGCGACAACCCGGCTATCGGGGAAAGTTCGGGCGTAGAAGTTCGCCGCCTCTTCGGCATCGTGGTCGTACCACAGGCAGATGGTGTT  >  NC_002947/4125530‑4125819

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