Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	386,451	T→C	pseudogene (745/1527 nt)	BW25113_RS01915 →	outer membrane autotransporter barrel domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	386,451	0	T	C	100.0%	128.5 / NA	34	pseudogene (745/1527 nt)	BW25113_RS01915	outer membrane autotransporter barrel domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/16);  total (18/16)

AAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCC  >  NZ_CP009273/386368‑386538                                                                                    |                                                                                        aaGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTc                                                                                  >  2:482297/1‑91 (MQ=255) aaGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTc                                                                                  >  1:448610/1‑91 (MQ=255)      aaTACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAAtt                                                                             >  1:629186/1‑91 (MQ=255)       aTACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATta                                                                            <  2:263280/91‑1 (MQ=255)         aCTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATGACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATtata                                                                          <  1:48763/91‑1 (MQ=255)              cAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATcc                                                                       <  2:540057/89‑1 (MQ=255)              cAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATcc                                                                       >  1:540057/1‑89 (MQ=255)              cAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTg                                                                     >  2:416753/1‑91 (MQ=255)                      ccTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCaa                                                             <  1:390174/91‑1 (MQ=255)                       cTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAAt                                                            <  1:234518/91‑1 (MQ=255)                              aCCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATg                                                     >  2:29808/1‑91 (MQ=255)                              aCCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATg                                                     >  1:276818/1‑91 (MQ=255)                                             tGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACt                                      >  2:424941/1‑91 (MQ=255)                                                 aTAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTgg                                  >  1:719942/1‑91 (MQ=255)                                                     cAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAAt                              <  2:719942/91‑1 (MQ=255)                                                         atCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACtt                                                        >  1:37387/1‑61 (MQ=255)                                                         atCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACtt                                                        <  2:37387/61‑1 (MQ=255)                                                         atCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACtt                                                        <  2:186949/61‑1 (MQ=255)                                                         atCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACtt                                                        >  1:186949/1‑61 (MQ=255)                                                         atCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGCATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCtg                          >  1:814368/1‑91 (MQ=255)                                                             ctGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCtgatga                      >  1:864891/1‑91 (MQ=255)                                                               gTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAAc                    >  1:229838/1‑91 (MQ=255)                                                                  aaaGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGa                 <  2:276818/91‑1 (MQ=255)                                                                  aaaGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGa                 <  2:864891/91‑1 (MQ=255)                                                                   aaGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAg                <  1:29808/91‑1 (MQ=255)                                                                       ttCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGAt            <  1:168469/91‑1 (MQ=255)                                                                       ttCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGAt            <  2:438238/91‑1 (MQ=255)                                                                       ttCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGAt            >  2:846939/1‑91 (MQ=255)                                                                         ccACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTa          >  1:63243/1‑91 (MQ=255)                                                                            caGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGc       <  2:707807/91‑1 (MQ=255)                                                                                gACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCtgatga                      >  1:621899/1‑72 (MQ=255)                                                                                gACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCtgatga                      <  2:621899/72‑1 (MQ=255)                                                                                gACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGc   >  2:96648/1‑91 (MQ=255)                                                                                 aCTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGcc  <  2:159427/91‑1 (MQ=255)                                                                                    |                                                                                        AAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCC  >  NZ_CP009273/386368‑386538

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