Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,250,887	G→C	pseudogene (639/986 nt)	BW25113_RS21020 →	DUF2713 domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,250,887	0	G	C	100.0%	157.9 / NA	42	pseudogene (639/986 nt)	BW25113_RS21020	DUF2713 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (11/31);  total (11/31)

TTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGT  >  NZ_CP009273/4250800‑4250973                                                                                        |                                                                                       ttAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTccc                                                                                      <  1:690789/90‑1 (MQ=255)      gTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGc                                                                                 <  2:230323/90‑1 (MQ=255)           ttAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGAtt                                                                            <  2:66365/90‑1 (MQ=255)           ttAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGAtt                                                                            <  2:281266/90‑1 (MQ=255)           ttAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGAtt                                                                            <  1:480622/90‑1 (MQ=255)           ttAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGAtt                                                                            <  2:273607/90‑1 (MQ=255)                  ctgccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGa                                                                     <  1:331262/90‑1 (MQ=255)                   tgccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGaa                                                                    >  2:558974/1‑90 (MQ=255)                   tgccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGaa                                                                    >  1:211751/1‑90 (MQ=255)                   tgccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGaa                                                                    <  2:2098/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  2:170270/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  2:211751/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  1:743851/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  1:66409/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  1:64194/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  1:513191/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  2:520793/90‑1 (MQ=255)                     ccAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAAc                                                                  <  2:529437/90‑1 (MQ=255)                       aGGATCATGATGACAACAGCTTAAAATTATTTCCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTCCGCAGATTTATCAGAAACAg                                                                <  1:151887/90‑1 (MQ=255)                       aGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAg                                                                <  1:788686/90‑1 (MQ=255)                           aCATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGta                                                            >  1:630875/1‑90 (MQ=255)                                atgaCAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTa                                                       <  2:156884/90‑1 (MQ=255)                                  gaCAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAAt                                                     <  1:768442/90‑1 (MQ=255)                                      acaGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGaaa                                                 <  2:258192/90‑1 (MQ=255)                                      acaGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGaaa                                                 >  2:241345/1‑90 (MQ=255)                                       caGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGaaaa                                                <  1:763509/90‑1 (MQ=255)                                         gCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAg                                                                <  2:311162/72‑1 (MQ=255)                                         gCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAg                                                                >  1:311162/1‑72 (MQ=255)                                         gCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAaata                                              >  2:595989/1‑90 (MQ=255)                                           ttAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAaataat                                            <  2:328423/90‑1 (MQ=255)                                               aaTTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATgaga                                        <  2:797990/90‑1 (MQ=255)                                                 ttattTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGc                                      >  2:737216/1‑90 (MQ=255)                                                   attTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTa                                    >  1:249322/1‑90 (MQ=255)                                                    tttGCTTCGCTGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTAc                                   <  2:22026/90‑1 (MQ=255)                                                              tGTTAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTc                         <  2:773376/90‑1 (MQ=255)                                                                ttAATTTCACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAg                       <  2:793433/90‑1 (MQ=255)                                                                      cctcTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGaa                 <  1:73332/87‑1 (MQ=255)                                                                       cACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGaaa                >  1:378040/1‑90 (MQ=255)                                                                       cACTTTCATATTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGaaa                >  1:103738/1‑90 (MQ=255)                                                                                atTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTAtctc       <  2:219286/90‑1 (MQ=255)                                                                                atTCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTAtctc       <  2:216591/90‑1 (MQ=255)                                                                                 ttCTAGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCa      <  2:378040/90‑1 (MQ=255)                                                                                     aGTCCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGt  >  1:202995/1‑90 (MQ=255)                                                                                        |                                                                                       TTAACGTTTGTTAATGCCTGCCAGGAACATGATGACAACAGCTTAAAATTATTTGCTTCGCTGTTAATTTCACTTTCATATTCTAGTGCCTACGCAGATTTATCAGAAACAGTGTATATTAATGAAAATAATGAGAGCTACCTGAAAGCTCAGTTTGAAAAATTATCTCAACGT  >  NZ_CP009273/4250800‑4250973

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