Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,376,328 G→C 12.8% intergenic (+31/+18) kup → / ← PP_1201 potassium transport system protein Kup/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,376,3280GC12.8% 112.0 / 7.5 39intergenic (+31/+18)kup/PP_1201potassium transport system protein Kup/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/18);  new base C (3/2);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GTCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTAC  >  NC_002947/1376180‑1376471
                                                                                                                                                    |                                                                                                                                               
gTCGCGCTGGCGGGGCAGCCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGTTCGTGCTCTGAATGGGATTATACAGCAGGCTCAGGGCCGCagg                                                                                                                                                >  2:366508/1‑149 (MQ=255)
gTCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGTCCGCaga                                                                                                                                                <  1:997423/150‑1 (MQ=255)
 tCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCagag                                                                                                                                               >  2:771594/1‑150 (MQ=255)
 tCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCagag                                                                                                                                               >  1:389379/1‑150 (MQ=255)
 tCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCagag                                                                                                                                               >  1:389735/1‑150 (MQ=255)
    cgcTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATACCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCgg                                                                                                                                            <  2:1014279/150‑1 (MQ=255)
        ggCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                        <  1:889798/150‑1 (MQ=255)
                  ccTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTc                                                                                                                              >  1:1108517/1‑150 (MQ=255)
                            ttCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCATGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTc                                                                                                                    >  1:901711/1‑150 (MQ=255)
                               ttGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGc                                                                                                                 >  1:86990/1‑150 (MQ=255)
                                            ccAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGACTGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTg                                                                                                    >  2:306575/1‑150 (MQ=255)
                                                               cttcAACCTGCCGTTGAACCGGGTGATCGAGCGAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCGGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtct                                                                                 >  2:329429/1‑150 (MQ=255)
                                                                cgcaagttGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGAtctc                                                                                <  2:671389/143‑1 (MQ=255)
                                                                          cGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGTAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAAc                                                                      <  1:636051/150‑1 (MQ=255)
                                                                           gTTGAACCGGGTGATCGAGCTAGCGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATAAAcc                                                                     >  2:552064/1‑150 (MQ=255)
                                                                               aaCCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGGGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCtt                                                                 <  2:389735/150‑1 (MQ=255)
                                                                               aaCCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCtt                                                                 <  2:389379/150‑1 (MQ=255)
                                                                                aCCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCATTGTCAGTTTTGCTAAGCTACGCTCAAAGCGGCTTGGCGCGTTTCGATCTCTCCAAGCAACCGCTTg                                                                >  1:119721/1‑150 (MQ=255)
                                                                                  cGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGc                                                              <  1:783051/150‑1 (MQ=255)
                                                                                    ggTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGa                                                            >  1:593906/1‑150 (MQ=255)
                                                                                           gAGCTCGGGACCGAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                     <  2:923184/150‑1 (MQ=255)
                                                                                           gAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCgg                                                     <  2:1108517/150‑1 (MQ=255)
                                                                                               tAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCCGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAg                                                 >  1:529117/1‑150 (MQ=255)
                                                                                               tAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                        <  1:611865/63‑1 (MQ=255)
                                                                                               tAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCt                                                                                                                                        >  2:611865/1‑63 (MQ=255)
                                                                                               tAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAg                                                 >  1:1107643/1‑150 (MQ=255)
                                                                                                 gggACCCAGGTCGAGATCTGAATGGTGTTATACATCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGtt                                               <  2:529117/150‑1 (MQ=255)
                                                                                                  ggaacgaGGGCGAGATCTGAATGGGGTTATACAGCAGGGTGAGGGCCGCAGAGCGGCACTGAGCCTTGTCAGTGCTGGTCAGCTACGCTGGAGTTGCCTGGGCGGGTTTCGATCTCGCCAATCAACCGCTTGGTGAGTGCCGGGTAGTtc                                              <  1:366508/144‑1 (MQ=255)
                                                                                                    aCCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                      <  1:565418/60‑1 (MQ=25)
                                                                                                    aCCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCTCTGCGGCCCTGa                                                                                                                                      >  2:565418/1‑60 (MQ=25)
                                                                                                      ccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgt                                          <  1:739538/150‑1 (MQ=255)
                                                                                                      ccAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCtcgt                                          >  2:184518/1‑150 (MQ=255)
                                                                                                             gagaTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTc                                                                                      <  2:253614/99‑1 (MQ=255)
                                                                                                             gagaTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTc                                                                                      >  1:253614/1‑99 (MQ=255)
                                                                                                               gaTCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGTGAGTGCCGGGTAGTTCTCGTCGAAgtggt                                 <  2:947114/150‑1 (MQ=255)
                                                                                                                 tCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGc                               >  1:939012/1‑150 (MQ=255)
                                                                                                                   tGAATGGGGTTATACAGCAGGCTCAGGGCCGCaga                                                                                                                                                <  2:991575/35‑1 (MQ=255)
                                                                                                                   tGAATGGGGTTATACAGCAGGCTCAGGGCCGCaga                                                                                                                                                >  1:991575/1‑35 (MQ=255)
                                                                                                                   tGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGccc                             >  2:314415/1‑150 (MQ=255)
                                                                                                                        ggggTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAg                        <  1:184518/150‑1 (MQ=255)
                                                                                                                                  agcagGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAg              <  2:593906/150‑1 (MQ=255)
                                                                                                                                        cTCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAAGCAACCCCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGGGGTGCCCGCCAGGCAGGTTCAGGCCCTc        >  1:1011827/1‑150 (MQ=255)
                                                                                                                                         tCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGcc                              <  2:716925/127‑1 (MQ=255)
                                                                                                                                         tCAGGGCCGCTCTGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGcc                              >  1:716925/1‑127 (MQ=255)
                                                                                                                                         tCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCa       >  1:760030/1‑150 (MQ=255)
                                                                                                                                              gCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTAc  <  1:552064/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTCGCGCTGGCGGGGCAACCTGTTCGCGTTCTTGCTGAAGAATGCCAACGGCAACTTGCGCTTCTTCAACCTGCCGTTGAACCGGGTGATCGAGCTAGGGACCCAGGTCGAGATCTGAATGGGGTTATACAGCAGGCTCAGGGCCGCAGAGCGGCCCTGAGCCTTGTCAGTTCTGCTCAGCTACGCTGGACTTGCCTTGGCGCGTTTCGATCTCGCCAATCAACCGCTTGGCGAGTGCCGGGTAGTTCTCGTCGAAGTGGTGCCCGCCAGGCAGTTTCAGGCGCTCACCTAC  >  NC_002947/1376180‑1376471

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