Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,520,464	A→G	8.5%	intergenic (+55/+17)	cspD → / ← icd	DNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,520,464	0	A	G	8.5%	123.6 / 2.6	47	intergenic (+55/+17)	cspD/icd	DNA replication inhibitor/NADP(+)‑specific isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (16/27);  new base G (2/2);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAGGTC  >  NC_002947/4520323‑4520612                                                                                                                                              |                                                                                                                                                     gCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGttt                                                                                                                                               >  3:159897/1‑149 (MQ=255)    gTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACCGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAACCGGTCGCTTCTTTTAAAGAAGCGGCCGAtttttt                                                                                                                                            <  8:294200/149‑3 (MQ=38)    gTTGAAATCAAAGATGCCCCCGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCt                                                                                                                                            <  8:284651/149‑1 (MQ=255)     ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTg                                                                                                                                           <  4:254691/149‑1 (MQ=255)     ttGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTg                                                                                                                                           <  7:258713/149‑1 (MQ=255)        aaaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGtt                                                                                                                                         <  6:294213/148‑1 (MQ=255)         aaTCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTg                                                                                                                                        <  4:159897/148‑1 (MQ=255)           tCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCt                                                                                                                                     >  4:293192/1‑149 (MQ=255)             cacgATGCCGCGGCCGCTGCGGCATAGGCTGACAGAGCCTTGGGCGAGCCGCAAAGCGAAACTCCAAACGCCTGACTGCCTGCGCTCTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGACGCGACCGGTTTTCTGTTGCCtt                                                                                                                                    <  6:189901/145‑1 (MQ=255)             aaaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTa                                                                                                                                   <  2:125256/149‑1 (MQ=255)              aaGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTAc                                                                                                                                  <  1:339970/149‑1 (MQ=255)               aGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACa                                                                                                                                 <  1:317841/149‑1 (MQ=255)                    ccGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATgtgt                                                                                                                            <  5:162860/149‑1 (MQ=255)                    ccGCGCCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATgtgt                                                                                                                            <  8:337367/149‑1 (MQ=255)                      gcgGCCGCTGCGGCATAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGttt                                                                                                                          <  6:6123/149‑1 (MQ=255)                          ccGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGAt                                                                                                                       >  1:213328/1‑148 (MQ=255)                          ccGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATc                                                                                                                      <  3:293192/149‑1 (MQ=255)                             cTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCAtt                                                                                                                   <  1:85594/149‑1 (MQ=255)                                            ccAGAGCCTTTGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCt                                                                                                    <  7:177532/149‑1 (MQ=255)                                            ccAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCt                                                                                                    <  2:213328/149‑1 (MQ=255)                                               gagCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAg                                                                                                 >  7:116128/1‑149 (MQ=255)                                                      gggCGAGCCGCAAAGCGAAAATGCAAACGCCTGACTGCCTCCGCTGTGCCTATAGACAAAAAAAGCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAgctgct                                                                                            <  6:209614/147‑1 (MQ=255)                                                               gCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt                                                                                                                           <  4:32063/107‑1 (MQ=255)                                                               gCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGtt                                                                                                                           >  3:32063/1‑107 (MQ=255)                                                                           tGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTAc                                                                                                                                  >  4:34180/1‑88 (MQ=255)                                                                           tGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTAc                                                                                                                                  <  3:34180/88‑1 (MQ=255)                                                                            gCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGt                                                                                                                                          <  4:54623/79‑1 (MQ=255)                                                                            gCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGt                                                                                                                                          >  3:54623/1‑79 (MQ=255)                                                                            gCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTcc                                                                         >  1:210905/1‑144 (MQ=255)                                                                            gCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTcc                                                                         <  2:210905/144‑1 (MQ=255)                                                                                  gCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTc                                                              <  4:248618/149‑1 (MQ=255)                                                                                     tGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGa                                                            >  3:246756/1‑148 (MQ=255)                                                                                           cctgcGGTGTGCCGATAGAAAAAAAAATCGGCCGATTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTa                                                     <  3:185910/149‑1 (MQ=255)                                                                                                  tgtgCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACg                                              >  8:149479/1‑149 (MQ=255)                                                                                                          tAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGc                                      <  8:116128/149‑1 (MQ=255)                                                                                                               aaaaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCaa                                 >  1:313810/1‑149 (MQ=255)                                                                                                                 aaaaaTCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTCCCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATg                               >  5:199030/1‑149 (MQ=255)                                                                                                                        ggtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAt                        >  7:208663/4‑149 (MQ=255)                                                                                                                         gtcGCTTCTTTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCAtt                       >  1:193391/3‑149 (MQ=255)                                                                                                                          ccGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTg                      <  7:149479/149‑1 (MQ=255)                                                                                                                            gCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAgg                                                   >  6:94269/1‑118 (MQ=255)                                                                                                                            gCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAgg                                                   <  5:94269/118‑1 (MQ=255)                                                                                                                              ttcttTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGc                  >  8:2457/1‑149 (MQ=255)                                                                                                                                cttTTAAAGAAGCGGCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGcc                 <  8:208663/148‑1 (MQ=255)                                                                                                                                 tttAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCt               >  2:218708/1‑149 (MQ=255)                                                                                                                                    aaaaGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCTCGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTtga            >  2:183411/1‑149 (MQ=255)                                                                                                                                           gCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAg     <  6:199030/149‑1 (MQ=255)                                                                                                                                              aCCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAGGTc  >  1:72291/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     GCCGTTGAAATCAAAGATGCCGCGGCCGCTGCGGCAGAGGCTGCCAGAGCCTTGGGCGAGCCGCAAAGCGAAACTGCAAAAGCCTGACTGCCTGCGCTGTGCCGATAGACAAAAAAATCGGCCGCTTCTTTAAAAGAAGCGACCGGTTTTCTGTTGCCTTACATGTGTTTGATCATTTCATCGCCAAAGCCTGAGCTGCTTACCAGCTTGGCGCCCTCCATCAATCGCTCGAAATCGTAGGTCACGGTCTTGGCCGCAATGGCGCCATTGGTGCCCTTGATGATCAGGTC  >  NC_002947/4520323‑4520612

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: ‑