Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,769,064	A→G	pseudogene (74/185 nt)	BW25113_RS24190 →	DUF905 domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,769,064	0	A	G	100.0%	99.0 / NA	28	pseudogene (74/185 nt)	BW25113_RS24190	DUF905 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/16);  total (12/16)

GCAGGGCCACTGGTATCTGCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGATGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCAT  >  NZ_CP009273/2768982‑2769150                                                                                   |                                                                                       gCAGGGCCACTGGTATCTGCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGaa                                                                                 <  1:530424/90‑1 (MQ=255)    gggCCACTGGTATCTGCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGt                                                                              >  1:94913/1‑90 (MQ=255)      gCCACTGGTATCTTCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCg                                                                            >  1:26965/1‑90 (MQ=255)                 cTGCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCg                                                                 <  1:726794/90‑1 (MQ=255)                 cTGCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCg                                                                 <  1:764484/90‑1 (MQ=255)                 cTGCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCg                                                                 <  1:798108/90‑1 (MQ=255)                             ccGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGc                                                                  >  2:715981/1‑77 (MQ=255)                             ccGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGc                                                                  <  1:715981/77‑1 (MQ=255)                                aCATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGa                                                  >  2:872475/1‑90 (MQ=255)                                  aTTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAAc                                                <  2:859025/90‑1 (MQ=255)                                       cATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGAc                                           <  1:391318/90‑1 (MQ=255)                                       cATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGAc                                           <  1:13307/90‑1 (MQ=255)                                          tGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGcc                                        >  2:838760/1‑90 (MQ=255)                                                  acCAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAgg                                >  1:817429/1‑90 (MQ=255)                                                    cAGGGCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGtt                              >  1:77304/1‑90 (MQ=255)                                                      gggCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGtttt                            >  2:1023976/1‑90 (MQ=255)                                                        gCATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTaa                          <  1:607320/90‑1 (MQ=255)                                                         cATTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAAc                         <  1:280950/90‑1 (MQ=255)                                                          aTTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCgg                                             >  1:62580/1‑69 (MQ=255)                                                          aTTCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCgg                                             <  2:62580/69‑1 (MQ=255)                                                           ttCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCg                       <  2:525508/90‑1 (MQ=255)                                                           ttCACGTCCGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCg                       <  2:745888/90‑1 (MQ=255)                                                                  ccGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCg                       <  2:878783/83‑1 (MQ=255)                                                                  ccGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCg                       >  1:878783/1‑83 (MQ=255)                                                                   cGCCTGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATcc               <  2:77304/90‑1 (MQ=255)                                                                       tGATAGTTCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAg           <  2:817429/90‑1 (MQ=255)                                                                              tCGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGc    >  2:590122/1‑90 (MQ=255)                                                                               cGTGGTGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCa   >  2:109368/1‑90 (MQ=255)                                                                                gtggtgCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCAt  <  2:26965/90‑1 (MQ=255)                                                                                   |                                                                                       GCAGGGCCACTGGTATCTGCCGGTATCACCGACATTAACATTGAAGACGACCAGGGCATTCACGTCCGCCTGATAGTTCGTGATGCCGAAGGTCGGATGGTTGGCGGGCATGGAACTTTGAACCGGACGCCGGTGAAGGTTTTAACCGCTATATCCACAGGTCAGGCAT  >  NZ_CP009273/2768982‑2769150

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