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%	289.5 / NA	89	pseudogene (312/921 nt)	yjiP	pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (52/37);  total (52/37)

ATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGAAG  >  NC_000913/4569185‑4569444                                                                                                                             |                                                                                                                                        atatTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCa                                                                                                                           <  1:436769/139‑1 (MQ=255)         tGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAAc                                                                                                                   >  1:998009/1‑139 (MQ=255)         tGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAAc                                                                                                                   >  1:1226827/1‑139 (MQ=255)            tCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGcc                                                                                                                >  2:160041/1‑139 (MQ=255)              ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                                  <  2:1261168/119‑1 (MQ=255)              ggTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                                  >  1:1261168/1‑119 (MQ=255)                tAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGcc                                                                                                                >  1:914603/1‑135 (MQ=255)                tAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGcc                                                                                                                <  2:914603/135‑1 (MQ=255)                 aaaGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           <  1:611961/139‑1 (MQ=255)                          cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCata                                                                                                                                  >  1:945449/1‑107 (MQ=255)                          cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATagag                                                                                                                             >  1:953021/1‑112 (MQ=255)                          cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATagag                                                                                                                             <  2:953021/112‑1 (MQ=255)                          cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTg                                                                                                  >  1:28377/1‑139 (MQ=255)                          cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTg                                                                                                  >  2:1038850/1‑139 (MQ=255)                          cgAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCACTATTCCATGGCGGTGATGTAACGCCata                                                                                                                                  <  2:945449/107‑1 (MQ=255)                                 gATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           <  2:28377/139‑1 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               >  1:959611/1‑105 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               <  2:959611/105‑1 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               >  1:1250461/1‑105 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCa                                                                                                               <  2:1250461/105‑1 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           >  2:40032/1‑125 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           <  1:40032/125‑1 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  1:1142630/1‑139 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  2:117585/1‑139 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  1:449738/1‑139 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  2:82596/1‑139 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  2:163304/1‑139 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  2:189627/1‑139 (MQ=255)                                               tgtgGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTAt                                                                             >  2:755675/1‑139 (MQ=255)                                                gtggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATc                                                                            >  1:826629/1‑139 (MQ=255)                                                 tggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                      <  2:1006895/128‑1 (MQ=255)                                                 tggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATg                                                                                      >  1:1006895/1‑128 (MQ=255)                                                  ggtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCAc                                                                          <  2:998009/139‑1 (MQ=255)                                                   gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                       >  2:199291/1‑77 (MQ=255)                                                   gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTa                                                                                                        >  2:474200/1‑108 (MQ=255)                                                   gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTa                                                                                                        <  1:474200/108‑1 (MQ=255)                                                   gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACg                                                                         >  1:1080355/1‑139 (MQ=255)                                                   gtgATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACg                                                                         >  1:1208825/1‑139 (MQ=255)                                                     gATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           >  2:96217/1‑103 (MQ=255)                                                     gATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           <  1:96217/103‑1 (MQ=255)                                                     gATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGt                                                                       >  1:990675/1‑139 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAacg                                                                                                                                       >  1:1129558/1‑74 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           <  1:150779/102‑1 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCg                                                                                                           >  2:150779/1‑102 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           <  1:1109496/118‑1 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           >  2:1109496/1‑118 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATTCCGa                                                                                        >  2:872166/1‑121 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGa                                                                                        <  1:872166/121‑1 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGc                                                                                     <  1:247847/124‑1 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGc                                                                                     >  2:247847/1‑124 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCa                                                                           >  2:913801/1‑134 (MQ=255)                                                      aTTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCa                                                                           <  1:913801/134‑1 (MQ=255)                                                       ttGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCGCGATGCTATTTTATCACGGTAg                                                                     >  2:569363/1‑139 (MQ=255)                                                       ttGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAg                                                                     >  1:248050/1‑139 (MQ=255)                                                          aaCATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                                  <  2:248050/139‑1 (MQ=255)                                                           aCATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGt                                                                 >  2:402070/1‑139 (MQ=255)                                                             aTCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                                  >  2:90325/1‑136 (MQ=255)                                                             aTCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCg                                                                  <  1:90325/136‑1 (MQ=255)                                                             aTCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAg                                                               >  2:1265204/1‑139 (MQ=255)                                                                 gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                           <  2:1226827/139‑1 (MQ=255)                                                                 gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                           <  2:449738/139‑1 (MQ=255)                                                                 gagCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCt                                                           <  1:189627/139‑1 (MQ=255)                                                                       cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           >  2:224680/1‑101 (MQ=255)                                                                       cgAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATcc                                                                                           <  1:224680/101‑1 (MQ=255)                                                                           gACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGt                                                 <  1:1265204/139‑1 (MQ=255)                                                                                      ggCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTgctgg                                      >  2:323758/1‑139 (MQ=255)                                                                                      ggCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTgctgg                                      <  1:82596/139‑1 (MQ=255)                                                                                          tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTc                                                <  2:615490/125‑1 (MQ=255)                                                                                          tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTc                                                >  1:615490/1‑125 (MQ=255)                                                                                          tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGctggctgg                                  >  2:145610/1‑139 (MQ=255)                                                                                          tttCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGctggctgg                                  >  2:805322/1‑139 (MQ=255)                                                                                            tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTAcc                                                                                                      >  2:376742/1‑69 (MQ=255)                                                                                            tCGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTAcc                                                                                                      <  1:376742/69‑1 (MQ=255)                                                                                             cGCCTGATGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATCAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCAGTGCTGGCTGGACg                               >  1:793202/1‑139 (MQ=255)                                                                                                   aTGCGCTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTg                         <  1:402070/139‑1 (MQ=255)                                                                                                        cTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGtt                                                                                                   <  1:611984/60‑1 (MQ=255)                                                                                                        cTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGtt                                                                                                   >  2:611984/1‑60 (MQ=255)                                                                                                        cTATTCCATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGccgac                    <  1:1038850/139‑1 (MQ=255)                                                                                                              cATGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGACTGGACGAATTTGCCGACCCGACt              >  2:525906/1‑139 (MQ=255)                                                                                                               aTGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTa             <  1:755675/139‑1 (MQ=255)                                                                                                                tGGCGGTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTAc            >  2:73956/1‑139 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGAt                                                                                       <  2:559954/59‑1 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGAt                                                                                       >  1:559954/1‑59 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGt                                                                       <  2:447361/75‑1 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGt                                                                       >  1:447361/1‑75 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc       >  1:653117/1‑139 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc       >  1:1283464/1‑139 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc       >  2:665550/1‑139 (MQ=255)                                                                                                                     gTGATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCAc       >  2:377376/1‑139 (MQ=255)                                                                                                                       gATGTAACGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACgg     <  2:793202/139‑1 (MQ=255)                                                                                                                          gTAACGCCATATAGAGCATGATAAACGCCAGCCGCTAACGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGAAg  <  2:826629/139‑1 (MQ=255)                                                                                                                             |                                                                                                                                        ATATTTTATGGTCGGTAAAGACCCGCGAAGGCGATGGCTATATCTATGTGGTGATTGAACATCAGAGCCGCGAGGACATTCATATGGCCTTTCGCCTGATGCGCTATTCCATGGCGGTGATGTAGCGCCATATAGAGCATGATAAACGCCAGCCGCTACCGTTGGTCATCCCGATGCTATTTTATCACGGTAGCCGTAGTCCTTACCCCTGGTCCCTGTGCTGGCTGGACGAATTTGCCGACCCGACTACCGCACGGAAG  >  NC_000913/4569185‑4569444

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: ‑