Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,210,906 G→A 11.7% intergenic (‑181/‑58) PP_2814 ← / → PP_2815 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,210,9060GA11.7% 84.0 / 3.4 34intergenic (‑181/‑58)PP_2814/PP_2815hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/14);  new base A (2/2);  total (18/16)
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.32e-01

GCAGGCGGGTTGACGTTTGATCAACGGGCATCAGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGCTTATGCGCAAAGGCGCCGAACTCGCTCAGGCGCTGC  >  NC_002947/3210758‑3211051
                                                                                                                                                    |                                                                                                                                                   
gcAGGCGGGTTGACGTTTGATCAACGGGCATCAGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTg                                                                                                                                                     >  4:114388/1‑149 (MQ=255)
      gggTTGACGTTTGATCAACGGGCATCAGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGcc                                                                                                                                                 >  7:210154/1‑147 (MQ=255)
       ggTTGACGTTTGATCAACGGTCATCAGTAAACAGGCCTGCTTGACACTGCGCTAAACCCAAAAAAGTGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGAAGTGAGTTCATTCTAGGTTTCCACTGGCCCTGCGCCAAg                                                                                                                                              <  8:210154/149‑1 (MQ=255)
              gTTTGATCAACGGGCATCAGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGaa                                                                                                                                       >  5:111392/1‑149 (MQ=255)
                              tCAGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTACGCCAAGGTACGAACTGCACCATTTTCCAg                                                                                                                      >  2:192787/1‑149 (MQ=255)
                              tCAGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTACGCCAAGGTACGAACTGCACCATTTTCCAg                                                                                                                      >  2:197109/1‑149 (MQ=255)
                                aGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGcc                                                                                                                    >  1:122728/1‑149 (MQ=255)
                                aGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGc                                                                                                                     <  2:244701/148‑1 (MQ=255)
                                          ccTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAg                                                                                                          >  5:104651/1‑149 (MQ=255)
                                                 gCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGAtt                                                                                                     <  6:111392/147‑1 (MQ=255)
                                                       gcgcTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGc                                                                                            >  1:137216/1‑149 (MQ=255)
                                                           tAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAAc                                                                                                                                      <  8:244349/105‑1 (MQ=255)
                                                           tAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAAc                                                                                                                                      >  7:244349/1‑105 (MQ=255)
                                                           tAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGTACCGCCCAt                                                                                        >  6:158252/1‑149 (MQ=255)
                                                             aaCCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAg                                                                                                                                              <  8:61969/95‑1 (MQ=255)
                                                             aaCCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAg                                                                                                                                              >  7:61969/1‑95 (MQ=255)
                                                             aaCCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATg                                                                                       <  3:114388/148‑1 (MQ=255)
                                                              aCCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGAc                                                                                     <  3:182498/149‑1 (MQ=255)
                                                              aCCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGAc                                                                                     <  5:245981/149‑1 (MQ=255)
                                                                       aGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCCCGcc                                                                            >  3:248740/1‑149 (MQ=255)
                                                                                  gCGTAAGGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTgg                                                                  <  7:187141/148‑1 (MQ=255)
                                                                                       aCGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGc                                                            <  4:125325/149‑1 (MQ=255)
                                                                                       aCGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGc                                                            >  3:122636/1‑149 (MQ=255)
                                                                                       aCGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGc                                                            >  3:125325/1‑149 (MQ=255)
                                                                                        cGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTACGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCCATGACCAGCACGCCGTCGCCCTGGCCACGc                                                            <  1:197109/149‑1 (MQ=255)
                                                                                        aggTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTACGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCCATGACCAGCACGCCGTCGCCCTGGCCACGc                                                            <  1:192787/148‑1 (MQ=255)
                                                                                                  ttGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACAGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAAc                                                 >  2:89543/1‑149 (MQ=255)
                                                                                                              ttCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGc                                      <  4:122636/148‑1 (MQ=255)
                                                                                                                           ttCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGcccc                                                         >  2:292887/1‑116 (MQ=255)
                                                                                                                           ttCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGcccc                                                         <  1:292887/116‑1 (MQ=255)
                                                                                                                              aTGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGTACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGCTTATGCGCAAAGgcgc                      <  5:158252/148‑1 (MQ=255)
                                                                                                                                gCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTATTCCAGCCGATTTTTGAGTGATTTGCCCCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGCTTATGCGCAAAGGCGcc                     >  1:62732/1‑148 (MQ=255)
                                                                                                                                   aGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGCTTATGCGCAAAGGCGCCGAACt                >  3:289904/1‑149 (MQ=255)
                                                                                                                                    ggTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGCTTATGCGCAAAGGCGCCGAACTc               <  2:137216/149‑1 (MQ=255)
                                                                                                                                                 ccTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGCTTATGCGCAAAGGCGCCGAACTCGCTCAGGCGCTGc  >  2:286703/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAGGCGGGTTGACGTTTGATCAACGGGCATCAGGAAACAGGCCTGCTTGCCACTGCGCTAAACCCAAAAAAGGGACATGAAGCGTCACGGTGCCTCCTTGTATTTTTTGTTCAGGCAGTGAGTTCATGCTAGGTTTCCACTGGCCCTGCGCCAAGGTACGAACTGCACCATTTTCCAGCCGATTTTTGAGTGATTTGCACCGCCCATGACCAGCACGCCGTCGCCCTGGCCACGCCCCTCGGAAACCGTTCGCGAGCTTATGCGCAAAGGCGCCGAACTCGCTCAGGCGCTGC  >  NC_002947/3210758‑3211051

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: