Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	539,260	C→A	G38G (GGC→GGA)	allB →	allantoinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	539,260	0	C	A	100.0%	107.8 / NA	31	G38G (GGC→GGA)	allB	allantoinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (15/16);  total (15/16)

TTTACAAGGAGTTTGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGCGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCATATTTCTGAACCGGGTCGTAGCCACTGGGAAG  >  NC_000913/539130‑539360                                                                                                                                   |                                                                                                     tttACAAGGAGTTTGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCaaaa                                                                                             <  2:17865/140‑1 (MQ=255)   tACAAGGAGTTTGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAaaga                                                                                           <  1:54172/140‑1 (MQ=255)     cAAGGAGTTTGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCaaaa                                                                                             <  2:280560/136‑1 (MQ=255)     cAAGGAGTTTGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCaaaa                                                                                             >  1:280560/1‑136 (MQ=255)        ggAGTTTGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGTAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTa                                                                                      >  2:305192/1‑140 (MQ=255)            tttGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGa                                                                                  <  2:194726/140‑1 (MQ=255)                 tATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgt                                                                             <  1:182279/140‑1 (MQ=255)                    gTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgtctg                                                                          >  1:346632/1‑140 (MQ=255)                    gTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgtctg                                                                          >  2:244196/1‑140 (MQ=255)                    gTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgtctg                                                                          >  1:158825/1‑140 (MQ=255)                    gTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgtctg                                                                          >  2:191133/1‑140 (MQ=255)                      cTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgtctggt                                                                        >  2:63169/1‑140 (MQ=255)                      cTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgtctggt                                                                        >  1:58905/1‑140 (MQ=255)                      cTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCgtctggt                                                                        >  2:16853/1‑140 (MQ=255)                          tGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGtctggtctgg                                                                    >  1:253080/1‑140 (MQ=255)                            aTTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCtggtg                                                                  >  2:159671/1‑140 (MQ=255)                               tAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGtt                                                               <  1:305192/140‑1 (MQ=255)                                 aTCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTc                                                             >  2:158202/1‑140 (MQ=255)                                 aTCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTc                                                             >  2:209900/1‑140 (MQ=255)                                  tCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCg                                                            <  1:244196/140‑1 (MQ=255)                                         aaaCGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCa                                                     <  2:40953/140‑1 (MQ=255)                                            cGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATgg                                                  <  1:16853/140‑1 (MQ=255)                                               cACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTg                                               <  1:191133/140‑1 (MQ=255)                                                    ttATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATgcg                                          <  2:158825/140‑1 (MQ=255)                                                     tATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATgcgc                                         <  1:63169/140‑1 (MQ=255)                                                      aTTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGca                                        <  2:346632/140‑1 (MQ=255)                                                            gAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCa                                  <  1:36473/140‑1 (MQ=255)                                                            gAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCa                                  >  2:36473/1‑140 (MQ=255)                                                                                   tatCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCATATTTCTGAACCGGGTCGTAGcc           <  2:58905/140‑1 (MQ=255)                                                                                         cGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCATATTTCTGAACCGGGTCGTAGCCACTggg     <  2:253080/140‑1 (MQ=255)                                                                                         cGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCATATTTCTGAACCGGGTCGTAGCCACTggg     >  1:286190/1‑140 (MQ=255)                                                                                            tAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGAGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCATATTTCTGAACCGGGTCGTAGCCACTGGGAAg  >  1:4213/1‑140 (MQ=255)                                                                                                                                   |                                                                                                     TTTACAAGGAGTTTGTTATGTCTTTTGATTTAATCATTAAAAACGGCACCGTTATTTTAGAAAACGAAGCTCGCGTTGTAGATATCGCCGTTAAAGGCGGAAAAATTGCTGCTATCGGTCAGGATCTGGGCGATGCAAAAGAAGTTATGGATGCGTCTGGTCTGGTGGTTTCGCCGGGCATGGTTGATGCGCACACCCATATTTCTGAACCGGGTCGTAGCCACTGGGAAG  >  NC_000913/539130‑539360

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