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%	130.8 / NA	37	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 (15/22);  total (15/22)

ACATCAATATTTCCAATAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATG  >  NZ_CP009273/386320‑386583                                                                                                                                    |                                                                                                                                     aCATCAATATTTCCAATAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGtt                                                                                                                                <  1:545043/138‑1 (MQ=255) aCATCAATATTTCCAATAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTc                                                                                                                               <  2:231103/139‑1 (MQ=255)      aaTATTTCCAATAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAAtt                                                                                                                          <  1:80660/139‑1 (MQ=255)            tCCAATAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATcc                                                                                                                    >  2:54629/1‑139 (MQ=255)              cAATAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTg                                                                                                                  <  1:81647/139‑1 (MQ=255)               aaTAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTgg                                                                                                                 <  2:134618/139‑1 (MQ=255)                aTAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGa                                                                                                                >  2:384334/1‑139 (MQ=255)                aTAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGa                                                                                                                <  2:506589/139‑1 (MQ=255)                       ccGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGa                                                                                                                <  1:469802/132‑1 (MQ=255)                       ccGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGa                                                                                                                >  2:469802/1‑132 (MQ=255)                       ccGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCa                                                                                                           <  2:291380/137‑1 (MQ=255)                                            aTCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGa                                                                                                                >  2:286404/1‑111 (MQ=255)                                            aTCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGa                                                                                                                <  1:286404/111‑1 (MQ=255)                                                        tACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCt                                                                        >  1:349847/1‑139 (MQ=255)                                                             aCAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCtgat                                                                     <  2:178491/137‑1 (MQ=255)                                                                      ccTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGa                                                          <  1:54629/139‑1 (MQ=255)                                                                      ccTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGa                                                          <  2:182534/139‑1 (MQ=255)                                                                               ccGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACt                                                                                                      >  1:40185/1‑86 (MQ=255)                                                                               ccGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACt                                                                                                      <  2:40185/86‑1 (MQ=255)                                                                                cGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGa                                                          >  1:467390/1‑129 (MQ=255)                                                                                cGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGa                                                          <  2:467390/129‑1 (MQ=255)                                                                                  aaTCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGccc                                              >  1:181651/1‑139 (MQ=255)                                                                                    tctcATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCa                                                                                                           >  1:88144/1‑76 (MQ=255)                                                                                    tctcATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCa                                                                                                           <  2:88144/76‑1 (MQ=255)                                                                                    tctcATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCgg                                            >  1:354665/1‑139 (MQ=255)                                                                                       cATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGt                                          >  1:412379/1‑138 (MQ=255)                                                                                               ggATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCa                                 <  2:180291/139‑1 (MQ=255)                                                                                               ggATAACAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCa                                 <  2:181651/139‑1 (MQ=255)                                                                                                     cAATATCTCTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTcc                           >  1:65703/1‑139 (MQ=255)                                                                                                          tctcTGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTAc                      >  2:512443/1‑139 (MQ=255)                                                                                                              tGTAAAAGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACt                     <  2:65703/136‑1 (MQ=255)                                                                                                                  aaaGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCg                                                               >  2:196771/1‑90 (MQ=255)                                                                                                                  aaaGATTCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCg                                                               <  1:196771/90‑1 (MQ=255)                                                                                                                        tCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCa                                 >  2:520505/1‑114 (MQ=255)                                                                                                                        tCCACAGTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCa                                 <  1:520505/114‑1 (MQ=255)                                                                                                                              gTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATg  <  2:305765/139‑1 (MQ=255)                                                                                                                              gTGACTCCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATg  <  2:188529/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                     ACATCAATATTTCCAATAACTCCCGCGTTGCAGGTATTGCATTATCTCAAGGTAATACTTACAACGAAACCTACACTACCGAATCTCATACTTGGGATAACAATATCTCTGTAAAAGATTCCACAGTGACTTCGGGTTCAAATTATATCCTGGATAGCAATACTTATGGCAAAACTGGTCACTTTGGCAATTCTGATGAACCGAGTGATTATGCTGGCCCGGGTGATGTTGCAATGTCCTTTACTGCTTCAGGTTCCGACTATG  >  NZ_CP009273/386320‑386583

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