Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,461,422	(G)5→4	pseudogene (1774/2526 nt)	BW25113_RS07350 →	autotransporter domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,461,418	0	G	.	100.0%	56.1 / NA	14	pseudogene (1770/2526 nt)	BW25113_RS07350	autotransporter domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (6/8);  total (6/8)

AATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGAC  >  NZ_CP009273/1461284‑1461540                                                                                                                                       |                                                                                                                           aaTCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAgggg                                                                                                                         <  1:514607/138‑1 (MQ=255)       gATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACTGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACa                                                                                                                   <  1:351507/137‑1 (MQ=255)                                 tGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTa                                                                                                                >  1:370959/1‑114 (MQ=255)                                 tGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTa                                                                                                                <  2:370959/114‑1 (MQ=255)                                 tGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGt                                                                                         >  1:272506/1‑137 (MQ=255)                                               agaACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACa                                                                         <  1:245100/139‑1 (MQ=255)                                                    aTTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAggggcaa                                                                                                                      <  1:448740/90‑4 (MQ=255)                                                    aTTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAggggcaa                                                                                                                      >  2:448740/1‑87 (MQ=255)                                                    aTTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAgggg                                                                                                                         >  1:421831/1‑87 (MQ=255)                                                    aTTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAgggg                                                                                                                         <  2:421831/87‑1 (MQ=255)                                                    aTTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGtt                                                                    >  1:195259/1‑139 (MQ=255)                                                         tGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCa                                                               <  1:339618/139‑1 (MQ=255)                                                                                 aGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATc                                                             <  1:467867/117‑1 (MQ=255)                                                                                 aGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATc                                                             >  2:467867/1‑117 (MQ=255)                                                                                   gagaTATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAAACAACAGTGGGAAGGTTTCg                                     >  1:368476/1‑139 (MQ=255)                                                                                            ttgttTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTggcagg                             <  2:272506/138‑1 (MQ=255)                                                                                                 tCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAgggg                                                                                                                         >  1:546402/1‑42 (MQ=255)                                                                                                 tCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAgggg                                                                                                                         <  2:546402/42‑1 (MQ=255)                                                                                                                aTTGGCGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGAttt        >  1:12201/1‑139 (MQ=255)                                                                                                                     cGTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGa   <  2:195259/139‑1 (MQ=255)                                                                                                                      gTGCTCATCAACGGTGA‑GGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTTAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGGAGGCAGTATGCAGGTTGGCGATTTCTCGAc  <  1:561880/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                           AATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTCCGATACTGGATCGATTGGCGTGCTCATCAACGGTGAGGGGGCAACAGTATCCAATACGGGTGATGTTAACGTTAGCAATGAAGCGACAGGGTTCAGCATCACAACCAACAGTGGGAAGGTTTCGCTGGCAGGCAGTATGCAGGTTGGCGATTTCTCGAC  >  NZ_CP009273/1461284‑1461540

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