Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,500,814	Δ1 bp	intergenic (‑23/‑446)	yjgX ← / → yjgZ	pseudogene fragment/uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,500,811	0	T	.	100.0%	314.7 / NA	67	intergenic (‑20/‑449)	yjgX/yjgZ	pseudogene fragment/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (38/29);  total (38/29)

TCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAGGGG  >  NC_000913/4500713‑4500905                                                                                                   |                                                                                               tCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtt                                                                                               >  1:600612/1‑100 (MQ=255)  ccGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGttt                                                                                              >  2:142990/1‑100 (MQ=255)    tctaGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttct                                                                                            <  2:316343/97‑3 (MQ=255)    gCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttc                                                                                             <  2:229011/99‑2 (MQ=255)     ccAGAAAATCACCACGAATACCAGAATCGCGGCACTGATCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCTTGTAGTGtttctg                                                                                           <  2:223048/100‑4 (MQ=255)      cAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATAATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttctg                                                                                           >  2:879747/1‑96 (MQ=255)       aGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGtttctaga                                                                                         >  1:316343/1‑95 (MQ=255)            aTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTc                                                                                   >  1:1376652/1‑100 (MQ=255)            aTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTc                                                                                   <  1:872417/100‑1 (MQ=255)            aTCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTc                                                                                   <  1:872429/100‑1 (MQ=255)              caccacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTg                                                                                 >  2:102820/1‑100 (MQ=255)              caccacGAAGACCAGAATCGGGGCACTGTCCAGAAAATGATCCTCAGGAATCACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTg                                                                                 <  2:25663/100‑1 (MQ=255)               accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGa                                                                                >  2:110691/1‑100 (MQ=255)               accacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGa                                                                                <  1:843444/100‑1 (MQ=255)                ccacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGaa                                                                               >  2:1218358/1‑100 (MQ=255)                 cacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAAc                                                                              >  2:1361764/1‑100 (MQ=255)                 cacGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAAc                                                                              >  1:824925/1‑100 (MQ=255)                  acGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAAcc                                                                             >  2:171920/1‑100 (MQ=255)                  acGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAAcc                                                                             <  1:387400/100‑1 (MQ=255)                    gAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCgg                                                                           <  1:1381227/100‑1 (MQ=255)                      aTACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGtt                                                                         <  2:316297/100‑1 (MQ=255)                         ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTa                                                                      >  2:49175/1‑100 (MQ=255)                         ccAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TGTCTGTAGGTTCTGAACCGGTTCTa                                                                      >  2:1324919/1‑100 (MQ=255)                              aTCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCggg                                                                 >  1:1004571/1‑100 (MQ=255)                                    gCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTtgtg                                                           <  1:142990/100‑1 (MQ=255)                                     cACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGa                                                          <  2:1322590/100‑1 (MQ=255)                                          aCCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTg                                                     <  1:1125665/100‑1 (MQ=255)                                           ccAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGa                                                    >  1:824249/1‑100 (MQ=255)                                            cAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGAt                                                   >  2:97437/1‑100 (MQ=255)                                            cAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGAt                                                   >  1:21663/1‑100 (MQ=255)                                             ngAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                  >  1:143654/2‑100 (MQ=255)                                             aGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATg                                                  >  1:1263972/1‑100 (MQ=255)                                              gAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                 <  1:1281029/100‑1 (MQ=255)                                              gAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                 >  2:200316/1‑100 (MQ=255)                                              gAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                 <  2:1362145/100‑1 (MQ=255)                                              gAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATgg                                                 <  1:277690/100‑1 (MQ=255)                                               aaaaTGATAAACAGTAAACACTGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGc                                                <  1:1490490/100‑1 (MQ=255)                                                aaaTGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCGAGCGGGGTTGTGAGGCTGATGGCa                                               <  2:1180381/100‑1 (MQ=255)                                                      tAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGcc                                         >  2:297111/1‑100 (MQ=255)                                                      tAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGcc                                         >  1:1332159/1‑100 (MQ=255)                                                       aaaCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCg                                        >  2:639679/1‑100 (MQ=255)                                                         aCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCggg                                      <  1:897986/100‑1 (MQ=255)                                                         aCAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCggg                                      >  1:1286901/1‑100 (MQ=255)                                                          cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                     >  2:869168/1‑100 (MQ=255)                                                          cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                     <  1:1361771/100‑1 (MQ=255)                                                          cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                     >  2:869175/1‑100 (MQ=255)                                                          cAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGc                                     <  1:102820/100‑1 (MQ=255)                                                            ggAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGctc                                   >  1:1209407/1‑100 (MQ=255)                                                             gAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGctct                                  >  2:665732/1‑100 (MQ=255)                                                                 cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCTGGGCTCTTCCt                              >  2:1062764/1‑100 (MQ=255)                                                                 cacGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCt                              >  2:1312746/1‑100 (MQ=255)                                                                   cGGGGCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGc                            <  1:1421169/100‑1 (MQ=255)                                                                    ggggCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCa                           <  1:49175/100‑1 (MQ=255)                                                                      ggCGGTTCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATg                         <  2:824244/100‑1 (MQ=255)                                                                           ttCACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTaca                    >  2:658583/1‑100 (MQ=255)                                                                             cACTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACt                  >  1:60819/1‑100 (MQ=255)                                                                               cTTTGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTgg                >  1:452203/1‑100 (MQ=255)                                                                                tttGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGt               >  1:1499583/1‑100 (MQ=255)                                                                                tttGCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGt               >  1:1499049/1‑100 (MQ=255)                                                                                   gCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACt            >  2:892464/1‑100 (MQ=255)                                                                                   gCTTTCTCCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACt            <  1:1077314/100‑1 (MQ=255)                                                                                       tctcCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCaa        >  2:556650/1‑100 (MQ=255)                                                                                        ctcCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAag       >  1:445836/1‑100 (MQ=255)                                                                                        ctcCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAag       >  2:1410497/1‑100 (MQ=255)                                                                                        ctcCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAag       <  2:500980/100‑1 (MQ=255)                                                                                         tcCTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAaga      >  2:811519/1‑100 (MQ=255)                                                                                          ccTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAagan     <  2:1261155/100‑2 (MQ=255)                                                                                          ccTGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAagag     <  1:280134/100‑1 (MQ=255)                                                                                            tGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAggg   >  1:1229843/1‑100 (MQ=255)                                                                                            tGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAggg   <  1:1410222/100‑1 (MQ=255)                                                                                            tGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAggg   <  1:1059574/100‑1 (MQ=255)                                                                                            tGTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATAGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAggg   <  1:1136758/100‑1 (MQ=255)                                                                                             gTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAgggg  >  1:404010/1‑100 (MQ=255)                                                                                             gTAGTG‑TTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAgggg  <  1:1312753/100‑1 (MQ=255)                                                                                                   |                                                                                               TCCGCCAGAAAATCACCACGAATACCAGAATCGCGGCACTGACCAGAAAATGATAAACAGGAAACACGGGGCGGTTCACTTTGCTTTCTCCTGTAGTGTTTTCTGTAGGTTCTGAACCGGTTCTAGCGGGGTTGTGAGGCTGATGGCAATCGCCGGGCTCTTCCTGCATGTTACACTGGTACTGCAAGAGGGG  >  NC_000913/4500713‑4500905

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