Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,867	+A	coding (660/1173 nt)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,866	1	.	A	97.1%	113.2 / ‑3.5	34	coding (661/1173 nt)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (20/13);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑AC‑A‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCA  >  NC_000913/4558730‑4558994                                                                                                                                           |                                                                                                                                  ccTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGc‑ac‑aat                                                                                                                                <  2:620869/140‑3 (MQ=255)         aGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGttt                                                                                                                        <  2:675632/140‑1 (MQ=255)               gACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCcagc                                                                                                                  <  2:603196/140‑1 (MQ=255)               gACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCcagc                                                                                                                  <  1:669526/140‑1 (MQ=255)               gACCGGTTCGTCAATCCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCcagc                                                                                                                  <  1:317875/140‑1 (MQ=255)                       cGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATa                                                                                                          >  2:837749/1‑140 (MQ=255)                       cGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTAGCTGATCGGCAAA‑A‑CCGCAGTGTTCCAGCAGGTCATa                                                                                                          >  1:248979/1‑140 (MQ=255)                        gTCAATGCTGCTGGTGATATCGTTGGTGCCGACTTAGCGAGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATaa                                                                                                         <  2:248979/140‑1 (MQ=255)                                   tgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTaa                                                                                              <  2:6011/140‑1 (MQ=255)                                   tgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTaa                                                                                              >  1:790542/1‑140 (MQ=255)                                     gTGATATCGATGGTGCCGCCTTCGCGCGCGAACTCTAGCGCCTGCTCATACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACg                                                                                            <  2:382912/140‑1 (MQ=255)                                          atCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCTGTTTTCCAGCACGTGATAAATGGGCTGTAACGCCttt                                                                                       >  1:352134/1‑140 (MQ=255)                                            cGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCttttt                                                                                     >  1:708053/1‑140 (MQ=255)                                                    ccgccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCg                                                                              <  2:856866/139‑1 (MQ=255)                                                    ccgccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCg                                                                              >  1:856866/1‑139 (MQ=255)                                                    ccgccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGc                                                                             >  2:743718/1‑140 (MQ=255)                                                    ccgccTTTGCGCGCGAACTCCACCGCCTGCTCAAAAAACCGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGc                                                                             <  2:266108/140‑1 (MQ=255)                                                     cgccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCGGGTCATAAATAGGCTGTAACGTCTTTTTACTGTCGcc                                                                            >  1:307844/1‑140 (MQ=255)                                                       ccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCa                                                                          >  2:771904/1‑140 (MQ=255)                                                       ccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCa                                                                          >  1:472934/1‑140 (MQ=255)                                                           tGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCAtgtg                                                                      >  2:654782/1‑140 (MQ=255)                                                            gcgcgcGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATCAATAGGCGGTAGCGACTTTTTACTGTCGCCCATGTgg                                                                     >  1:297133/1‑140 (MQ=255)                                                                          gcgcCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGcc                                                       <  1:837749/140‑1 (MQ=255)                                                                              cTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                                                   >  1:618014/1‑140 (MQ=255)                                                                              cTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGt                                                   >  2:707283/1‑140 (MQ=255)                                                                                                gTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCg                                 >  2:779694/1‑140 (MQ=255)                                                                                                  tGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGCACACGGTCACGTCATGTTTACTGCCGAGCAAACCGcc                               >  1:399855/1‑140 (MQ=255)                                                                                                    cGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc                               >  2:737555/1‑138 (MQ=255)                                                                                                    cGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGcc                               <  1:737555/138‑1 (MQ=255)                                                                                                         aaCGTGGGTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATATGGAACACCGTGACACCAGGTTTACCGCCGAGCAAACCGCCAACGCgg                        >  2:402245/1‑140 (MQ=255)                                                                                                               ggTCGGCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTTATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCt                  >  2:711301/1‑140 (MQ=255)                                                                                                                   ggCAGCAGCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgccg              <  2:618014/140‑1 (MQ=255)                                                                                                                     cagcagCTTGCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTgccgcc            >  2:184475/1‑140 (MQ=255)                                                                                                                        cagcTTGCTGATCGGC‑AG‑AGTCCCAGTTTTCTAGCCGGTTATAAATAGACTGTAACGCCTTTTTACTGTCGCCCATGGGGAACACGGTGACGCAGGGTTTACCGCCGAGCAACCCGCAAACGGGGGATTCTGCCGCCata         <  2:281397/140‑1 (MQ=255)                                                                                                                              gCTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGcc   <  1:654782/140‑1 (MQ=255)                                                                                                                               cTGATCGGC‑ACAA‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCa  >  1:102321/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                  CCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACAACGGTACGTTGCGGTTAACGTGGGTCGGCAGCAGCTTGCTGATCGGC‑AC‑A‑TCGCAGTTTTCCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTCGCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACCGCCAACGCGGGATTCTGCCGCCATATTGGCCA  >  NC_000913/4558730‑4558994

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