Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,041,206 G→A G294D (GGC→GAC)  polA → DNA polymerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,041,2060GA98.0% 193.9 / ‑6.0 51G294D (GGC→GAC) polADNA polymerase I
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/27);  total (23/28)
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.99e-01

TGGCGACGATTAAAACCGACGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGGCAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTACGTCACCATCCTTGATGAAGAAACAC  >  NZ_CP009273/4041071‑4041340
                                                                                                                                       |                                                                                                                                      
tGGCGACGATTAAAACCGACGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACa                                                                                                                                      >  1:27089/1‑138 (MQ=255)
       gaTTAAAACCGACGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTa                                                                                                                              <  1:475472/139‑1 (MQ=255)
       gaTTAAAACCGACGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTa                                                                                                                              <  1:1210241/139‑1 (MQ=255)
       gaTTAAAACCGACGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTa                                                                                                                              <  2:1117514/139‑1 (MQ=255)
                cGACGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCaaa                                                                                                                     <  2:758841/139‑1 (MQ=255)
                  aCGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAgg                                                                                                                   >  1:489457/1‑139 (MQ=255)
                             gAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAg                                                                                                        >  2:521224/1‑139 (MQ=255)
                                   aCCTGTGAACAACTGGAAGTGCAGCAACCGGCCGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCcgcg                                                                                                   >  1:1130038/1‑138 (MQ=255)
                                   aCCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCcgcg                                                                                                   >  1:767745/1‑138 (MQ=255)
                                   aCCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCcgc                                                                                                    <  2:27089/137‑1 (MQ=255)
                                     cTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAg                                                                                                >  1:1204035/1‑139 (MQ=255)
                                        tgAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCca                                                                                             <  1:257290/139‑1 (MQ=255)
                                        tgAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCca                                                                                             >  2:137363/1‑139 (MQ=255)
                                           acaacTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAg                                                                                          <  2:780530/139‑1 (MQ=255)
                                             aacTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTa                                                                                                                              >  1:304551/1‑101 (MQ=255)
                                             aacTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTa                                                                                                                              <  2:304551/101‑1 (MQ=255)
                                               cTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGaaa                                                                                      <  2:1045748/139‑1 (MQ=255)
                                               cTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGaaa                                                                                      <  2:489457/139‑1 (MQ=255)
                                                  gAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAAc                                                                                     >  1:359336/1‑137 (MQ=255)
                                                  gAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCa                                                                                   <  1:624210/139‑1 (MQ=255)
                                                     gTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAgtg                                                                                <  2:767745/139‑1 (MQ=255)
                                                         agcaACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGCAGCCACAGGAAACCAGTGTTg                                                                             >  2:1153332/1‑138 (MQ=255)
                                                            aaCCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGa                                                                         <  2:1094192/139‑1 (MQ=255)
                                                               cGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGa                                                                      <  2:931421/139‑1 (MQ=255)
                                                               cGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGa                                                                      <  1:419746/139‑1 (MQ=255)
                                                                                     gCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACg                                                >  2:536394/1‑139 (MQ=255)
                                                                                         ttCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGCAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCcac                                                                                            >  2:217973/1‑91 (MQ=255)
                                                                                         ttCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCcac                                                                                            <  1:217973/91‑1 (MQ=255)
                                                                                         ttCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGa                                            >  1:258/1‑139 (MQ=255)
                                                                                            aaaaaGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGAttt                                         <  1:713955/139‑1 (MQ=255)
                                                                                                aGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTa                                     <  1:536394/139‑1 (MQ=255)
                                                                                                  tATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATg                                   <  1:147213/139‑1 (MQ=255)
                                                                                                   aTGAGTTCAAACGCTGGACTGCTGATGTCTAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGa                                                                      >  1:62051/1‑103 (MQ=255)
                                                                                                   aTGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGa                                                                      <  2:62051/103‑1 (MQ=255)
                                                                                                   aTGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGa                                            >  1:1107558/1‑129 (MQ=255)
                                                                                                   aTGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGa                                            <  2:1107558/129‑1 (MQ=255)
                                                                                                   aTGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATg                                   >  2:898737/1‑138 (MQ=255)
                                                                                                   aTGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATg                                   <  1:898737/138‑1 (MQ=255)
                                                                                                        ttCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACt                             <  2:1204035/139‑1 (MQ=255)
                                                                                                         tCAAACGCTGGACTGCTGATGTCGAAGCGGACAAATGGTTCCGGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAATCAGTGTGGCAGACGAAGCCCCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTa                            <  2:65877/139‑1 (MQ=255)
                                                                                                          cAAACGCTGGACTGCTGATGTCGAAGCGGGCAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTa                            <  1:830629/138‑1 (MQ=255)
                                                                                                                cTGGACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGGGTTGCAGACGAAGCACCAGAAGTGACGGGAACGGTGATTTCTTATGACAACTACGTCAcc                     >  2:58839/1‑139 (MQ=255)
                                                                                                                  ggACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGc                                                    >  1:690458/1‑106 (MQ=255)
                                                                                                                  ggACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGc                                                    <  2:690458/106‑1 (MQ=255)
                                                                                                                  ggACTGCTGATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTACGTCACCAt                   >  2:523259/1‑139 (MQ=255)
                                                                                                                        ctgATGTCGAAGCGGACAAATGGTTACAGGCCGAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAg                                                                    >  1:782405/1‑84 (MQ=255)
                                                                                                                        ctgATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAg                                                                    <  2:782405/84‑1 (MQ=255)
                                                                                                                        ctgATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTACGTCACCATCCTtga             >  1:837131/1‑139 (MQ=255)
                                                                                                                        ctgATGTCGAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTACGTCACCATCCTtga             >  1:6990/1‑139 (MQ=255)
                                                                                                                                gAAGCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTACGTCACCATCCTTGATGAAGaaa     >  1:795189/1‑139 (MQ=255)
                                                                                                                                   gCGGACAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTACGTCACCATCCTTGATGAAGAAacac  <  2:1130038/139‑1 (MQ=255)
                                                                                                                                      ggCAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCca                                                                                             >  1:178320/1‑45 (MQ=255)
                                                                                                                                      ggCAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCca                                                                                             <  2:178320/45‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
TGGCGACGATTAAAACCGACGTTGAACTGGAGCTGACCTGTGAACAACTGGAAGTGCAGCAACCGGCAGCGGAAGAGTTGTTGGGGCTGTTCAAAAAGTATGAGTTCAAACGCTGGACTGCTGATGTCGAAGCGGGCAAATGGTTACAGGCCAAAGGGGCAAAACCAGCCGCGAAGCCACAGGAAACCAGTGTTGCAGACGAAGCACCAGAAGTGACGGCAACGGTGATTTCTTATGACAACTACGTCACCATCCTTGATGAAGAAACAC  >  NZ_CP009273/4041071‑4041340

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: