Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,281	A→G	pseudogene (806/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,281	0	A	G	100.0%	226.7 / NA	59	pseudogene (806/5422 nt)	yneO	AIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/28);  total (31/28)

CTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCT  >  NC_000913/1597144‑1597416                                                                                                                                          |                                                                                                                                        cTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTa                                                                                                                                       <  1:1189923/140‑1 (MQ=255)    tGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATc                                                                                                                                    <  2:1141164/140‑1 (MQ=255)           aaTCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTAtttt                                                                                                                             <  2:1102507/140‑1 (MQ=255)                gTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGcc                                                                                                                        >  2:306232/1‑140 (MQ=255)                gTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGcc                                                                                                                        <  2:698437/140‑1 (MQ=255)                gTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGcc                                                                                                                        >  1:202216/1‑140 (MQ=255)                 tCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCg                                                                                                                       <  1:38525/140‑1 (MQ=255)                 tCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCg                                                                                                                       >  2:211363/1‑140 (MQ=255)                 tCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCg                                                                                                                       >  1:723962/1‑140 (MQ=255)                        ttGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTaaccaa                                                                                                                <  2:239063/140‑1 (MQ=255)                          ggAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTaaccaacc                                                                                                              <  1:951207/140‑1 (MQ=255)                           gAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAAccc                                                                                                             >  1:159794/1‑140 (MQ=255)                           gAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAAccc                                                                                                             >  1:1279918/1‑140 (MQ=255)                                      tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGaa                                                                                                  >  1:339478/1‑140 (MQ=255)                                       gtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAAt                                                                                                 >  1:1265866/1‑140 (MQ=255)                                         gCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAg                                                                                               >  2:757463/1‑140 (MQ=255)                                          ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                                                              >  1:751358/1‑140 (MQ=255)                                          ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                                                              <  2:339478/140‑1 (MQ=255)                                             ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                                                                           <  2:1265866/140‑1 (MQ=255)                                             ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGc                                                                                           <  1:620009/140‑1 (MQ=255)                                              tACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGcc                                                                                          >  1:1084654/1‑140 (MQ=255)                                                 cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcc                                                                                       <  2:202216/140‑1 (MQ=255)                                                   aaaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcccc                                                                                     <  2:776935/140‑1 (MQ=255)                                                      tatCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGt                                                                                  <  1:528633/140‑1 (MQ=255)                                                        tCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTc                                                                                >  2:585148/1‑140 (MQ=255)                                                             tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                                                                           >  2:571452/1‑140 (MQ=255)                                                             tataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCACGAATAGTCGCCTCCCCAGTTCCTTGa                                                                           >  2:49087/1‑140 (MQ=255)                                                              ataGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCACGAATAGTCGCCTCCCCAGTTCCTTGAt                                                                          <  1:49087/140‑1 (MQ=255)                                                                 gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                                       <  1:306232/140‑1 (MQ=255)                                                                  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                                      >  2:606066/1‑140 (MQ=255)                                                                    tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                                    >  2:900477/1‑140 (MQ=255)                                                                      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaact                                                                  <  1:585148/140‑3 (MQ=255)                                                                      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                                  >  1:359439/1‑140 (MQ=255)                                                                      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAtt                                                                  >  2:757691/1‑140 (MQ=255)                                                                       gCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAttt                                                                 >  1:772507/1‑140 (MQ=255)                                                                        cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                                <  2:723962/140‑1 (MQ=255)                                                                           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                                             <  2:772507/140‑1 (MQ=255)                                                                           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                                             <  2:1084654/140‑1 (MQ=255)                                                                             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                                           >  1:1228377/1‑140 (MQ=255)                                                                             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                                           <  1:1259240/140‑1 (MQ=255)                                                                               tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                                       <  2:298371/126‑1 (MQ=255)                                                                               tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                                       >  1:298371/1‑126 (MQ=255)                                                                                         gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAAt                                               <  1:720701/140‑1 (MQ=255)                                                                                            ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                                            <  1:216/140‑1 (MQ=255)                                                                                                  cTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                      >  2:621167/1‑140 (MQ=255)                                                                                                   tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                                     >  1:422501/1‑140 (MQ=255)                                                                                                    gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                    >  2:1274590/1‑140 (MQ=255)                                                                                                    gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                                    >  2:66255/1‑140 (MQ=255)                                                                                                    gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAA‑CCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                                   >  2:513215/1‑140 (MQ=255)                                                                                                          tCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc                              >  1:908424/1‑140 (MQ=255)                                                                                                          tCAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc                              >  2:1187749/1‑140 (MQ=255)                                                                                                           cAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                             <  2:1279918/140‑1 (MQ=255)                                                                                                           cAGGGTTCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc                             <  2:873828/140‑1 (MQ=255)                                                                                                                ttCCGATACCGGTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg                        <  1:606066/140‑1 (MQ=255)                                                                                                                          ggTGGCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCtttt              <  1:1187749/140‑1 (MQ=255)                                                                                                                              gCATTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAAc          >  1:51750/1‑140 (MQ=255)                                                                                                                                aTTGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGa        >  1:810501/1‑140 (MQ=255)                                                                                                                                    ccACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc    <  2:1232092/140‑1 (MQ=255)                                                                                                                                      aCCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCt  >  2:695578/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        CTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCT  >  NC_000913/1597144‑1597416

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