Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,228,383 G→A 13.3% intergenic (‑127/+409) PP_5624 ← / ← PP_3707 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,228,3830GA13.3% 67.8 / 3.3 30intergenic (‑127/+409)PP_5624/PP_3707hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base A (2/2);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-02

CCCACATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACATTTCCGT  >  NC_002947/4228250‑4228517
                                                                                                                                     |                                                                                                                                      
cccacATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATg                                                                                                                         >  3:272191/1‑149 (MQ=255)
 ccacATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGc                                                                                                                        >  3:173426/1‑149 (MQ=255)
     aTGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGAc                                                                                                                     >  3:296417/1‑148 (MQ=255)
      tGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACaa                                                                                                                   <  5:161477/149‑1 (MQ=255)
              tCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTcga                                                                                                           >  5:100647/1‑149 (MQ=255)
                aGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTcgac                                                                                                          >  2:5963/1‑148 (MQ=255)
                          gAGCCAGGGATACTTTCCTCCCCTGTCCTTGGAGTGGCGAAGCTAGAAAGGTACGGTGATTTGCCCAACGATTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATCCCATTAg                                                                                               <  7:121068/149‑1 (MQ=255)
                             ccAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGc                                                                                             >  4:215664/1‑148 (MQ=255)
                                  gATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCAc                                                                                       <  4:173426/149‑1 (MQ=255)
                                               ccTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGa                                                                          <  6:100647/149‑1 (MQ=255)
                                                   tCCTTGGAGTTGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGGCATAACCTTGTTAAATGGAGTAt                                                                      <  2:268031/149‑1 (MQ=255)
                                                   tCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTAt                                                                      <  4:296417/149‑1 (MQ=255)
                                                              gCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGAtt                                                           >  3:176126/1‑149 (MQ=255)
                                                              gCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGAtt                                                           >  1:202477/1‑149 (MQ=255)
                                                              gCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGAt                                                            >  1:160592/1‑148 (MQ=255)
                                                                      aGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGcc                                                   >  6:136835/1‑149 (MQ=255)
                                                                      aGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGcc                                                   >  6:141291/1‑149 (MQ=255)
                                                                                ggTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCtt                                                                                    <  4:205423/106‑1 (MQ=255)
                                                                                ggTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTc                                         <  4:272191/149‑1 (MQ=255)
                                                                                ggTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATAACCtt                                                                                    >  3:205423/1‑106 (MQ=255)
                                                                                         cccAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGt                                 <  3:215664/148‑1 (MQ=255)
                                                                                                  ttttCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGa                        >  7:42601/1‑148 (MQ=255)
                                                                                                         gTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGc                   <  1:5963/146‑1 (MQ=255)
                                                                                                         gTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCAc                >  2:278678/1‑149 (MQ=255)
                                                                                                            cTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTa              <  1:278678/148‑1 (MQ=255)
                                                                                                               ccAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGa           >  1:258656/1‑148 (MQ=255)
                                                                                                               ccAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGa           <  5:136835/148‑1 (MQ=255)
                                                                                                               ccAAGCTGGCGTCTCCTATGCCAGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGa           <  5:141291/148‑1 (MQ=255)
                                                                                                                   gCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACAttt      >  5:297400/1‑149 (MQ=255)
                                                                                                                       gCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACATTTCCGt  <  6:297400/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
CCCACATGTTGAGCTCAGATGAATAGGAGCCAGGGATACTATCCTCCCCTGTCCTTGGAGTGGCGAACCTAGAAAGGTACGGTGATTTGCCCAACGAGTTTTCTAGTGCTCCCAAGCTGGCGTCTCCTATGCCGGCACGTCGCGTTATGCGACAAAGCCTCGACGATACCATTAGGCCATCACCTTGGTAAATGGAGTATGGGATCCGATTCTCCTGCCTTGGAGTGTCTGATCAGTGATACCAGATAGGCCACTACGACATTTCCGT  >  NC_002947/4228250‑4228517

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: