Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,466,996	T→C	pseudogene (1605/2513 nt)	ydbA →	pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,466,996	0	T	C	100.0%	156.8 / NA	44	pseudogene (1605/2513 nt)	ydbA	pseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/24);  total (20/24)

GCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466871‑1467124                                                                                                                              |                                                                                                                                 gCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                     <  2:61091/139‑1 (MQ=255)           aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgccgtggaatcttt                                                                                                           <  1:456852/139‑13 (MQ=37)                       agagCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCt                                                                                               >  2:295931/1‑139 (MQ=255)                        gagCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTg                                                                                              >  1:467336/1‑139 (MQ=255)                           cACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGAt                                                                                           <  1:708458/139‑1 (MQ=255)                              tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGt                                                                                        >  1:345925/1‑139 (MQ=255)                                 cACCAATGGCGGTACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAACATCGGCATTGTAAATCTGGATGGTAGc                                                                                     <  2:459295/139‑1 (MQ=255)                                 cACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGc                                                                                     <  2:464255/139‑1 (MQ=255)                                   ccAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCt                                                                                   <  2:169360/139‑1 (MQ=255)                                     aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                 >  1:171129/1‑139 (MQ=255)                                        ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGTCAATATCGGCAt                                                                                                        <  2:501220/113‑1 (MQ=255)                                        ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCAGGATGGTAGCCTGACTg                                                                              >  1:391831/1‑139 (MQ=255)                                        ggcggcACAGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCAt                                                                                                        >  1:501220/1‑113 (MQ=255)                                          cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  1:392798/1‑108 (MQ=255)                                          cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  2:392798/108‑1 (MQ=255)                                          cggcaccggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGtt                                                                            >  1:365960/1‑139 (MQ=255)                                                  gcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCgg                                                                    <  1:304603/139‑1 (MQ=255)                                                       cAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTgcgc                                                               >  2:242390/1‑139 (MQ=255)                                                        aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:295931/139‑1 (MQ=255)                                                        aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:171129/139‑1 (MQ=255)                                                        aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:703162/139‑1 (MQ=255)                                                            cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgt                                                          <  2:345925/139‑1 (MQ=255)                                                               cGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGtt                                                       >  1:619336/1‑139 (MQ=255)                                                               cGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATAACGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGtt                                                       <  2:301813/139‑1 (MQ=255)                                                                ggTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              >  1:514876/1‑131 (MQ=255)                                                                ggTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGATACAGGCGGTGCGCa                                                              <  2:514876/131‑1 (MQ=255)                                                                  tAAAGACGCCAACGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTgag                                                    <  1:499541/139‑1 (MQ=255)                                                                   aaCGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTgaga                                                   <  2:365960/139‑1 (MQ=255)                                                                    acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGaa                                                  >  2:476342/1‑139 (MQ=255)                                                                         gCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTg                                             >  2:523621/1‑139 (MQ=255)                                                                         gCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTg                                             >  2:464539/1‑139 (MQ=255)                                                                          ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTgg                                            <  1:126572/139‑1 (MQ=255)                                                                                  aaTAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACg                                    >  1:597084/1‑139 (MQ=255)                                                                                    tAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGc                                  <  2:597084/139‑1 (MQ=255)                                                                                            gAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaa                          <  2:764451/139‑1 (MQ=255)                                                                                             aaaaaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaac                         >  2:114463/1‑139 (MQ=255)                                                                                               aaaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaacaa                       >  1:538756/1‑139 (MQ=255)                                                                                                aaCCACTTTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaacaac                      <  2:270610/139‑1 (MQ=255)                                                                                                           aTGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGatat           >  2:398103/1‑139 (MQ=255)                                                                                                            tGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAtt          >  2:283268/1‑139 (MQ=255)                                                                                                               aaaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttgtt       <  1:283268/139‑1 (MQ=255)                                                                                                                    aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc  <  2:538756/139‑1 (MQ=255)                                                                                                                    aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc  >  2:582969/1‑139 (MQ=255)                                                                                                                    aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc  <  1:242390/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                 GCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466871‑1467124

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