Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,366	A→G	pseudogene (721/5422 nt)	yneO ←	AIDA‑I family autotransporter YneO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,366	0	A	G	100.0%	196.5 / NA	54	pseudogene (721/5422 nt)	yneO	AIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/31);  total (23/31)

TGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTG  >  NC_000913/1597242‑1597502                                                                                                                             |                                                                                                                                         tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTAttttt                                                                                                                           >  2:181838/1‑140 (MQ=255)  gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCGACCCCCCCTCCCAAATAGTCGCTTCCCCAGTTCCTTGTTTACAATTTTGAAATTCATTTGGGGAATCATTAtttttt                                                                                                                          >  2:681110/1‑140 (MQ=255)  gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTAtttttt                                                                                                                          >  2:481264/1‑140 (MQ=255)  gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTAtttttt                                                                                                                          >  1:792685/1‑140 (MQ=255)  gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTAtttttt                                                                                                                          >  2:90235/1‑140 (MQ=255)        tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGc                                                                                                                    >  2:540980/1‑140 (MQ=255)        tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGc                                                                                                                    <  1:278368/140‑1 (MQ=255)        tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGc                                                                                                                    >  2:788894/1‑140 (MQ=255)        tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGc                                                                                                                    <  1:532350/140‑1 (MQ=255)         cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGcc                                                                                                                   <  2:995187/140‑1 (MQ=255)                      ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTttatt                                                                                                                       <  2:194415/123‑1 (MQ=255)                      ccGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTttatt                                                                                                                       >  1:194415/1‑123 (MQ=255)                        ggTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCtttt                                                                                                    >  2:1023902/1‑140 (MQ=255)                         gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATATTCGCCTCCCCAGTTCCTTGATTCCCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTc                                                                                                   >  2:653992/1‑140 (MQ=255)                         gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTc                                                                                                   >  1:298273/1‑140 (MQ=255)                             cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg                                                                                               <  1:729279/140‑1 (MQ=255)                             cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg                                                                                               <  1:923856/140‑1 (MQ=255)                             cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACg                                                                                               <  1:992149/140‑1 (MQ=255)                              aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGa                                                                                              >  1:170184/1‑140 (MQ=255)                                  ccACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                                                                                          <  1:481264/140‑1 (MQ=255)                                    cccGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCt                                                                                        <  2:170184/139‑1 (MQ=255)                                    aCCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCt                                                                                        >  1:694061/1‑140 (MQ=255)                                      cGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                                                                                      <  1:788894/140‑1 (MQ=255)                                      cGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                                                                                      <  1:90235/140‑1 (MQ=255)                                       gATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGc                                                                                     >  2:894104/1‑140 (MQ=255)                                            tCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTa                                                                                <  2:587257/140‑1 (MQ=255)                                              gTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACt                                                                              <  1:1023902/140‑1 (MQ=255)                                                aTTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                                                                                          <  1:857650/126‑1 (MQ=255)                                                aTTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                                                                                          >  2:857650/1‑126 (MQ=255)                                                aTTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGc                                                                            <  2:36255/140‑1 (MQ=255)                                                                  cccccccTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTaaaa                                                          >  1:19472/1‑140 (MQ=255)                                                                      cccTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGc                                                      >  2:166988/1‑140 (MQ=255)                                                                       ccTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGcg                                                     <  1:937347/140‑1 (MQ=255)                                                                        cTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGcga                                                    <  2:19472/140‑1 (MQ=255)                                                                             gAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGAtt                                               >  2:311189/1‑140 (MQ=255)                                                                                   tCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTc                                         <  1:311189/140‑1 (MQ=255)                                                                                    cGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCg                                        <  1:166988/140‑1 (MQ=255)                                                                                        tCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTg                                              >  2:609427/1‑130 (MQ=255)                                                                                        tCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTg                                              <  1:609427/130‑1 (MQ=255)                                                                                          cccAGTTCCTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGt                                  >  1:355967/1‑140 (MQ=255)                                                                                                 ccTTGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCtt                           <  1:77517/140‑1 (MQ=255)                                                                                                    tGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATc                        >  2:917473/1‑140 (MQ=255)                                                                                                    tGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATc                        >  1:364385/1‑140 (MQ=255)                                                                                                    tGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATc                        >  2:898127/1‑140 (MQ=255)                                                                                                    tGATTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATc                        >  2:734526/1‑140 (MQ=255)                                                                                                      aTTACCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCCCCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCg                      <  2:694061/140‑1 (MQ=255)                                                                                                         aCCAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTaa                   >  2:697547/1‑140 (MQ=255)                                                                                                           cAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAata                 <  2:792685/140‑1 (MQ=255)                                                                                                           cAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAata                 <  1:540980/140‑1 (MQ=255)                                                                                                           cAATTTGAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAata                 <  2:364385/140‑1 (MQ=255)                                                                                                                 gAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAAt           <  1:384444/140‑1 (MQ=255)                                                                                                                 gAAATTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAAt           <  1:734526/140‑1 (MQ=255)                                                                                                                  aaaTTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATg          <  2:355967/140‑1 (MQ=255)                                                                                                                  aaaTTCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATg          <  2:1000280/140‑1 (MQ=255)                                                                                                                      tCAATTGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAAcc      <  1:697547/140‑1 (MQ=255)                                                                                                                          ttGGGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCtgtg  <  1:894104/140‑1 (MQ=255)                                                                                                                             |                                                                                                                                         TGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTG  >  NC_000913/1597242‑1597502

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