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%	140.6 / NA	40	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/20);  total (20/20)

CGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466861‑1467124                                                                                                                                        |                                                                                                                                 cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCAcc                                                                                                                               >  1:691274/1‑139 (MQ=255)                cGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAat                                                                                                                >  1:248651/1‑139 (MQ=255)                cGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAat                                                                                                                >  2:610890/1‑139 (MQ=255)                     aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  2:386572/1‑139 (MQ=255)                            aGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGACTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTa                                                                                                    <  1:364427/139‑1 (MQ=255)                             ggCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTaa                                                                                                   <  1:184047/139‑1 (MQ=255)                                gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTGACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   >  1:461262/1‑120 (MQ=255)                                gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   <  2:461262/120‑1 (MQ=255)                                gagaGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  2:773897/128‑1 (MQ=255)                                gagaGCACTATCACCAATAGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTATTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  1:773897/1‑128 (MQ=255)                                        tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTGCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGt                                                                                        >  1:196040/1‑139 (MQ=255)                                               aaTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                 <  2:751881/139‑1 (MQ=255)                                                        accggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAg                                                                        <  2:280650/139‑1 (MQ=255)                                                        accggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAg                                                                        >  2:327491/1‑139 (MQ=255)                                                          cggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGc                                                                      <  1:386572/139‑1 (MQ=255)                                                             cacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGt                                                                   <  1:686228/139‑1 (MQ=255)                                                                  aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:691274/139‑1 (MQ=255)                                                                  aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:648032/139‑1 (MQ=255)                                                                      cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              >  1:168960/1‑135 (MQ=255)                                                                      cAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:168960/135‑1 (MQ=255)                                                                       aaCGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGgtg                                                         <  1:395272/139‑1 (MQ=255)                                                                          ggTAACGACGCCACCGCGAATAACAGTGGAAAAACCTCTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTg                                                      >  2:712825/1‑139 (MQ=255)                                                                          ggTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  1:561901/1‑86 (MQ=255)                                                                          ggTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  2:561901/86‑1 (MQ=255)                                                                            tAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTgag                                                    >  1:50811/1‑139 (MQ=255)                                                                              acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGaa                                                  >  2:164948/1‑139 (MQ=255)                                                                                 acgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACAt                                               >  1:630343/1‑139 (MQ=255)                                                                                  cgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACAtt                                              >  1:355928/1‑139 (MQ=255)                                                                                      aCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATg                                                                                          >  1:745945/1‑91 (MQ=255)                                                                                      aCCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATg                                                                                          <  2:745945/91‑1 (MQ=255)                                                                                             aTAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACgg                                   >  1:726320/1‑139 (MQ=255)                                                                                                       aaaaaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaac                         >  2:180942/1‑139 (MQ=255)                                                                                                         aaaCCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaacaa                       >  2:732131/1‑139 (MQ=255)                                                                                                            ccACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAAcaacaa                    <  2:494302/139‑1 (MQ=255)                                                                                                                tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                 <  1:763191/74‑1 (MQ=255)                                                                                                                tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGa                                                                                 >  2:763191/1‑74 (MQ=255)                                                                                                                       ggAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttg         <  1:180942/139‑1 (MQ=255)                                                                                                                         aaaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttgtt       <  2:686689/139‑1 (MQ=255)                                                                                                                              aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc  <  2:50811/139‑1 (MQ=255)                                                                                                                              aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc  <  1:732131/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 CGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466861‑1467124

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