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%	209.7 / NA	58	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 (28/30);  total (28/30)

CGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466861‑1467124                                                                                                                                        |                                                                                                                                 cGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCAcc                                                                                                                               >  1:821465/1‑139 (MQ=255)      aCCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAAt                                                                                                                          >  1:1023454/1‑139 (MQ=255)          ggCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                      <  1:532916/139‑1 (MQ=255)           gCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                     <  2:1023454/139‑1 (MQ=255)           gCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCgg                                                                                                                     <  2:1116696/139‑1 (MQ=255)                     aTAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  2:1382415/139‑1 (MQ=255)                                   agCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTgg                                                                                             >  1:1042767/1‑139 (MQ=255)                                   agCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTgg                                                                                             >  2:1266697/1‑139 (MQ=255)                                        tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                      >  2:236412/1‑109 (MQ=255)                                        tATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcgcg                                                                                                                      <  1:236412/109‑1 (MQ=255)                                                aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATgg                                                                                         >  1:555354/1‑130 (MQ=255)                                                aTGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATgg                                                                                         <  2:555354/130‑1 (MQ=255)                                                 tGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAatat                                                                                                              <  2:1132259/108‑1 (MQ=255)                                                 tGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAatat                                                                                                              >  1:1132259/1‑108 (MQ=255)                                                  ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  2:44252/110‑1 (MQ=255)                                                  ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  2:1140877/1‑110 (MQ=255)                                                  ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  1:1140877/110‑1 (MQ=255)                                                  ggcggcACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  1:44252/1‑110 (MQ=255)                                                        accggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           <  1:976902/104‑1 (MQ=255)                                                        accggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCgg                                                                                                           >  2:976902/1‑104 (MQ=255)                                                          cggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCAGGATGGTAGCCTGACTGTTACAGGc                                                                      <  2:821465/139‑1 (MQ=255)                                                           ggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTAc                                                                          <  2:214011/134‑1 (MQ=255)                                                           ggcacTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTAc                                                                          >  1:214011/1‑134 (MQ=255)                                                                 cAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAA‑CACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:1241460/139‑1 (MQ=255)                                                                  aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:1266697/139‑1 (MQ=255)                                                                  aGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  1:102806/139‑1 (MQ=255)                                                                    aTCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATg                                                            >  2:991860/1‑139 (MQ=255)                                                                     tCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATgg                                                           <  1:991860/139‑1 (MQ=255)                                                                       aaCGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              >  1:1210757/1‑134 (MQ=255)                                                                       aaCGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCa                                                              <  2:1210757/134‑1 (MQ=255)                                                                         cGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGtt                                                       >  2:170811/1‑139 (MQ=255)                                                                            tAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTgag                                                    >  2:366807/1‑139 (MQ=255)                                                                              acgacgCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGaa                                                  >  1:522234/1‑139 (MQ=255)                                                                                   gCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                    >  1:824699/1‑68 (MQ=255)                                                                                   gCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCggg                                                                                                                    <  2:824699/68‑1 (MQ=255)                                                                                    ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                        >  2:350962/1‑63 (MQ=255)                                                                                    ccACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATcg                                                                                                                        <  1:350962/63‑1 (MQ=255)                                                                                          cgAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACaa                                      >  1:567841/1‑139 (MQ=255)                                                                                            aaTAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACg                                    >  2:69691/1‑139 (MQ=255)                                                                                               aaCAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCa                                 <  1:170811/139‑1 (MQ=255)                                                                                               aaCAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCa                                 <  2:567841/139‑1 (MQ=255)                                                                                                    tGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGtt                            >  1:279427/1‑139 (MQ=255)                                                                                                                tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTa                                                                           >  1:334991/1‑80 (MQ=255)                                                                                                                tGTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTa                                                                           <  2:334991/80‑1 (MQ=255)                                                                                                                 gTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   <  1:724584/39‑1 (MQ=255)                                                                                                                 gTTGATGGAAAAGATTCCACGGGCACCAAAATCGCGGGc                                                                                                                   >  2:724584/1‑39 (MQ=255)                                                                                                                    gATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAatat                                                                                                              <  1:612976/41‑1 (MQ=255)                                                                                                                    gATGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAatat                                                                                                              >  2:612976/1‑41 (MQ=255)                                                                                                                     aTGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGc                                                                      <  1:293670/80‑1 (MQ=255)                                                                                                                     aTGGAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGc                                                                      >  2:293670/1‑80 (MQ=255)                                                                                                                       ggAAAAGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGACTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGAAGATAttg         >  1:895821/1‑139 (MQ=255)                                                                                                                         aaaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATAttgtt       <  1:366807/139‑1 (MQ=255)                                                                                                                          aaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaa                          >  1:102658/1‑119 (MQ=255)                                                                                                                          aaaGATTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTaa                          <  2:102658/119‑1 (MQ=255)                                                                                                                              aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGTACGGTTAACAACAAAGGAGATATTGTTGTTTc  <  2:895821/139‑1 (MQ=255)                                                                                                                              aTTCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTc  <  2:522234/139‑1 (MQ=255)                                                                                                                               ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCa                                 >  2:39467/1‑107 (MQ=255)                                                                                                                               ttCCACGGGCACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCa                                 <  1:39467/107‑1 (MQ=255)                                                                                                                                        |                                                                                                                                 CGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCACGGGTACCAAAATCGCGGGCAATATCGGCATTGTAAATCTGGATGGTAGCCTGACTGTTACAGGCGGTGCGCATGGTGTTGAGAACATTGGTGACAACGGCACGGTTAACAACAAAGGAGATATTGTTGTTTC  >  NC_000913/1466861‑1467124

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