Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,569,309	G→A	pseudogene (312/921 nt)	yjiP →	pseudogene, transposase_31 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,569,309	0	G	A	100.0%	206.2 / NA	65	pseudogene (312/921 nt)	yjiP	pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/29);  total (36/29)

CTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGG  >  NC_000913/4569175‑4569441                                                                                                                                       |                                                                                                                                     cTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGcc                                                                                                                                  <  2:1242310/139‑1 (MQ=255)                 tATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATaa                                                                                                                  >  1:1477026/1‑139 (MQ=255)                        ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                               >  1:1338756/1‑119 (MQ=255)                        ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                               <  1:352185/119‑1 (MQ=255)                        ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                               >  2:352185/1‑119 (MQ=255)                        ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                               <  2:1338756/119‑1 (MQ=255)                        ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAg                                                                                                           >  1:492103/1‑139 (MQ=255)                           aaaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                        <  2:492103/139‑1 (MQ=255)                                      aaGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGt                                                                                             <  2:575409/139‑1 (MQ=255)                                               gCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGAt                                                                                    <  2:1477026/139‑1 (MQ=255)                                                        aTGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTa                                                                           >  1:1179044/1‑139 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                             <  1:1410128/88‑1 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCAtata                                                                                                                             >  2:1410128/1‑88 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGCGCATGATAAACGCCAGc                                                                                                          >  2:266584/1‑107 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGc                                                                                                          <  1:266584/107‑1 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                       >  2:1307182/1‑110 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGc                                                                                                       <  1:1307182/110‑1 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCAt                                                                                          >  1:1094510/1‑123 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCAt                                                                                          <  2:1094510/123‑1 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTAtattat                                                                          >  2:5917/1‑139 (MQ=255)                                                         tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                          >  1:978458/1‑139 (MQ=255)                                                          gtggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGc                                                                                  <  1:1333650/130‑1 (MQ=255)                                                          gtggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGc                                                                                  >  2:1333650/1‑130 (MQ=255)                                                             gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACg                                                                      >  1:183825/1‑139 (MQ=255)                                                             gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACg                                                                      >  1:1178566/1‑139 (MQ=255)                                                              tgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCGTGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACgg                                                                     >  2:57401/1‑139 (MQ=255)                                                              tgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACgg                                                                     >  1:81142/1‑139 (MQ=255)                                                                aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                            <  1:1320857/98‑1 (MQ=255)                                                                aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                            >  2:1320857/1‑98 (MQ=255)                                                                aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                   <  1:108286/123‑1 (MQ=255)                                                                aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                   >  2:108286/1‑123 (MQ=255)                                                                aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTa                                                                   >  1:1071568/1‑139 (MQ=255)                                                                      cATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTa                                                             <  2:1178566/139‑1 (MQ=255)                                                                           gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                        <  2:183825/139‑1 (MQ=255)                                                                           gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                        <  2:63801/139‑1 (MQ=255)                                                                           gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                        <  1:260834/139‑1 (MQ=255)                                                                                        aTTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                        >  2:1259741/1‑78 (MQ=255)                                                                                        aTTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                        <  1:1259741/78‑1 (MQ=255)                                                                                         ttCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCg                                                                                                  <  1:697229/83‑1 (MQ=255)                                                                                         ttCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCg                                                                                                  >  2:697229/1‑83 (MQ=255)                                                                                              atGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTgct                                     >  2:731593/1‑139 (MQ=255)                                                                                                  ccTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATa                                                                                                                   <  1:637737/57‑1 (MQ=255)                                                                                                  ccTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATa                                                                                                                   >  2:637737/1‑57 (MQ=255)                                                                                                    tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAacgc                                                                                                                                   >  1:41231/1‑39 (MQ=38)                                                                                                    tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAacgc                                                                                                                                   <  2:41231/39‑1 (MQ=255)                                                                                                      tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGt                                                                                             <  2:925883/75‑1 (MQ=255)                                                                                                      tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGt                                                                                             >  1:925883/1‑75 (MQ=255)                                                                                                                        cATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGa                              >  2:795346/1‑120 (MQ=255)                                                                                                                        cATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGa                              <  1:795346/120‑1 (MQ=255)                                                                                                                          tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTAc         >  1:649130/1‑139 (MQ=255)                                                                                                                          tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTAc         >  2:1037465/1‑139 (MQ=255)                                                                                                                           ggCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGcc                    <  1:963732/127‑1 (MQ=255)                                                                                                                           ggCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGcc                    >  2:963732/1‑127 (MQ=255)                                                                                                                            gCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGcc                    >  1:1433002/1‑126 (MQ=255)                                                                                                                            gCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGcc                    <  2:1433002/126‑1 (MQ=255)                                                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                   >  2:474537/1‑62 (MQ=255)                                                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                   <  1:474537/62‑1 (MQ=255)                                                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGc      >  2:431401/1‑139 (MQ=255)                                                                                                                             cGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGc      >  1:1207894/1‑139 (MQ=255)                                                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGcc                                                                <  2:270237/79‑1 (MQ=255)                                                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGcc                                                                >  1:270237/1‑79 (MQ=255)                                                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGcc                    <  2:930215/123‑1 (MQ=255)                                                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGcc                    >  1:930215/1‑123 (MQ=255)                                                                                                                               gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc    >  1:506847/1‑139 (MQ=255)                                                                                                                                 gATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACgg  <  1:57401/139‑1 (MQ=255)                                                                                                                                  aTGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                        >  2:321465/1‑52 (MQ=255)                                                                                                                                  aTGTAACGACATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                        <  1:321465/52‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     CTACACTCCGATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGG  >  NC_000913/4569175‑4569441

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