Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,506,680 G→A 10.5% intergenic (+91/‑108) PP_4839 → / → cycA iron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,506,6800GA10.5% 101.8 / 3.7 38intergenic (+91/‑108)PP_4839/cycAiron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease
Reads supporting (aligned to +/- strand):  ref base G (20/14);  new base A (2/2);  total (22/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 = 6.73e-01

GTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACGCAAACCTCACCCATAGCGACAGATGAACCGCAC  >  NC_002947/5506542‑5506826
                                                                                                                                            |                                                                                                                                                  
gTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGCGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                                                             >  1:198488/1‑146 (MQ=11)
gTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGCGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                                                             <  2:198488/146‑1 (MQ=11)
     tAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCt                                                                                                                                    <  1:553586/150‑1 (MQ=255)
       gTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa                                                                                                                                  <  1:647654/150‑1 (MQ=255)
                    gcgcCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCag                                                                                                                     >  1:289590/1‑150 (MQ=255)
                                            ttGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCt                                                                                             >  1:555193/1‑150 (MQ=255)
                                               aGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAg                                                                                                                                              <  2:32904/98‑1 (MQ=255)
                                               aGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAg                                                                                                                                              >  1:32904/1‑98 (MQ=255)
                                                 ccTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAc                                                                                        >  2:575058/1‑150 (MQ=255)
                                                  cTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCa                                                                                                                      >  2:264382/1‑119 (MQ=255)
                                                  cTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCa                                                                                                                      <  1:264382/119‑1 (MQ=255)
                                                          gCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCaga                                                                                                                    >  2:627728/1‑113 (MQ=255)
                                                          gCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCaga                                                                                                                    <  1:627728/113‑1 (MQ=255)
                                                            cTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTtag                                                                             >  1:616005/1‑150 (MQ=255)
                                                                gCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCGGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTtagccta                                                                         <  1:365940/150‑1 (MQ=255)
                                                                  tGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa                                                                                                                                  >  2:278679/1‑91 (MQ=255)
                                                                  tGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa                                                                                                                                  <  1:278679/91‑1 (MQ=255)
                                                                  tGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTg                                                                                            <  2:393507/129‑1 (MQ=255)
                                                                  tGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTg                                                                                            >  1:393507/1‑129 (MQ=255)
                                                                    gCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGATGTAGGAGCGGGTTGACCCGCGAAGAGGGCGGTAAAGGCGAAATGAGGAGCCGAAAAGGCAGCGTCGATCTGTCGCGCCGCGATCTGGAACGCGCGTGTtcgactagcag                                                                     >  2:202000/1‑148 (MQ=37)
                                                                       tctTCGCGGGTAAACCCGCTGCTACGAGGGACCGCGATGAAAGAGCGGGGTTACCCGCGATAGGGCCGGTACAGGCGAAATTCTGAACCGAAAAACCAGAGTCCATCTTTCGCGCCCCGCGCTGGAGCGCTCGTGTTAGCTTAGACGACt                                                                  >  1:256629/1‑150 (MQ=255)
                                                                       tctTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCCCGTGTTAGCCTAGCCGGCt                                                                  >  2:116671/1‑150 (MQ=255)
                                                                                    aCCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTcg                                                                                                        <  2:451280/99‑1 (MQ=255)
                                                                                    aCCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTcg                                                                                                        >  1:451280/1‑99 (MQ=255)
                                                                                     cccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaa                                                    >  1:364089/1‑150 (MQ=255)
                                                                                     cccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaagaa                                                    >  1:585778/1‑150 (MQ=255)
                                                                                       cGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAagaaga                                                  <  2:364089/150‑1 (MQ=255)
                                                                                       cGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGcaa                                                                                                                                          >  2:654949/1‑60 (MQ=17)
                                                                                       cGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGcaa                                                                                                                                          <  1:654949/62‑3 (MQ=17)
                                                                                                              tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCATAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACGCAAAcc                             <  2:289590/150‑1 (MQ=255)
                                                                                                                     gCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACg                                   >  1:515601/1‑137 (MQ=255)
                                                                                                                     gCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACg                                   <  2:515601/137‑1 (MQ=255)
                                                                                                                     gCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACGCAAACCTCACCCa                      >  1:113292/1‑150 (MQ=255)
                                                                                                                               ccGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGa                                     <  2:187611/125‑1 (MQ=255)
                                                                                                                               ccGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGa                                     >  1:187611/1‑125 (MQ=255)
                                                                                                                                cgcgAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACGCAAACCTCACCCATAGCGACAGAt           >  2:532696/1‑150 (MQ=255)
                                                                                                                                           cGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACGCAAACCTCACCCATAGCGACAGATGAACCGCAc  >  2:498499/1‑148 (MQ=255)
                                                                                                                                           cGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACGCAAACCTCACCCATAGCGACAGATGAACCGCAc  <  1:498499/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTATGACGCAAACCTCACCCATAGCGACAGATGAACCGCAC  >  NC_002947/5506542‑5506826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: