Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,519,014	G→A	pseudogene (294/504 nt)	insB1 →	IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,519,014	0	G	A	100.0%	211.8 / NA	77	pseudogene (294/504 nt)	insB1	IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (41/36);  total (41/36)

CGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGA  >  NC_000913/4518883‑4519144                                                                                                                                    |                                                                                                                                   cGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACgg                                                                                                                             >  2:72179/1‑139 (MQ=255)          tggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    <  1:72179/139‑1 (MQ=255)           ggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                                   <  1:170092/139‑1 (MQ=255)                     cgcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    <  2:814379/128‑1 (MQ=255)                     cgcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    >  1:814379/1‑128 (MQ=255)                          aTGACAGGCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc                                                                                                   >  2:321236/1‑139 (MQ=11)                          aTGACAGGCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc                                                                                                   >  2:997825/1‑139 (MQ=11)                          aTGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                                    >  2:1072341/1‑139 (MQ=255)                          aTGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                                    >  2:1359428/1‑139 (MQ=255)                           tGACAGGCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGcc                                                                                                  <  1:321236/139‑1 (MQ=11)                                  ctccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAggg                                                                                           >  1:463649/4‑139 (MQ=11)                                  aTACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                                            >  2:1248194/1‑139 (MQ=255)                                   tccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                                          <  1:648374/137‑1 (MQ=11)                                     cGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaa                                                                                         >  2:1064558/1‑139 (MQ=255)                                      ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCAGAAGGGaaa                                                                                        >  2:874925/1‑139 (MQ=255)                                        aggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    <  2:1016946/109‑1 (MQ=255)                                        aggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    >  1:1016946/1‑109 (MQ=255)                                        aggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGc                                                                                      >  2:486640/1‑139 (MQ=255)                                           aCGGTTGTGGCTCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCATGAAGGGAAAGCTGc                                                                                   >  1:120974/1‑139 (MQ=255)                                           aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  >  2:1336282/1‑108 (MQ=255)                                           aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  1:1336282/108‑1 (MQ=255)                                            cGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                                  <  1:1064558/139‑1 (MQ=255)                                                tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGa                                                                                                                            >  1:662359/1‑93 (MQ=255)                                                tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGa                                                                                                                            <  2:662359/93‑1 (MQ=255)                                                     cgcACGTATTGGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGTCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAACCAGc                                                                        <  1:300740/139‑1 (MQ=14)                                                      gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                        <  2:1043421/139‑1 (MQ=255)                                                      gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGATATCAGc                                                                        <  1:154986/139‑1 (MQ=255)                                                       cACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                      >  2:157428/1‑139 (MQ=14)                                                       cACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                      <  2:783460/139‑1 (MQ=14)                                                                  tGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                                      <  1:836747/97‑1 (MQ=255)                                                                  tGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                                      >  2:836747/1‑97 (MQ=255)                                                                  tGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       >  2:1145454/1‑79 (MQ=18)                                                                  tGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       <  1:1145454/79‑1 (MQ=18)                                                                    aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       <  2:435049/77‑1 (MQ=18)                                                                    aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       >  1:435049/1‑77 (MQ=17)                                                                    aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAg                                                         >  2:601588/1‑139 (MQ=11)                                                                      cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  >  1:832792/1‑80 (MQ=18)                                                                      cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                                  <  2:832792/80‑1 (MQ=17)                                                                      cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                                            <  1:338979/102‑1 (MQ=18)                                                                      cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                                            >  2:338979/1‑102 (MQ=17)                                                                             tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                <  1:242851/139‑1 (MQ=11)                                                                             tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                <  2:703487/139‑1 (MQ=11)                                                                              ggcgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAAt                                                                                                        >  1:345272/8‑82 (MQ=18)                                                                              ggcgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        >  1:755564/8‑114 (MQ=17)                                                                              ggcgacgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        <  2:755564/107‑1 (MQ=17)                                                                              ggCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAAt                                                                                                        <  2:345272/82‑1 (MQ=18)                                                                                   cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtat                                          >  2:368050/3‑135 (MQ=14)                                                                                   cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa                                           <  1:312275/139‑1 (MQ=255)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                        >  2:445355/1‑108 (MQ=255)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                        <  1:445355/108‑1 (MQ=255)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        >  1:1029216/1‑107 (MQ=18)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                                        <  2:1029216/107‑1 (MQ=18)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                                 <  1:571397/130‑1 (MQ=14)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                                 >  2:571397/1‑130 (MQ=14)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtataa                                        >  2:999538/1‑133 (MQ=14)                                                                                     cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATaa                                        >  1:478130/1‑139 (MQ=11)                                                                                         aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       <  2:60787/57‑1 (MQ=255)                                                                                         aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                                       >  1:60787/1‑57 (MQ=255)                                                                                           cGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGAGCATTGAGCGACATaatctga                                   >  2:803348/1‑139 (MQ=255)                                                                                             tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                                      <  1:203627/70‑1 (MQ=255)                                                                                             tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                                      >  2:203627/1‑70 (MQ=255)                                                                                             tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat                                 >  2:1300534/1‑139 (MQ=255)                                                                                                   gAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTgaga                           >  1:145835/1‑139 (MQ=255)                                                                                                         gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAt                     <  1:1072341/139‑1 (MQ=255)                                                                                                           ctgtggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                >  2:557332/5‑110 (MQ=11)                                                                                                           ctgtggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                                <  1:557332/106‑1 (MQ=11)                                                                                                               ggCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    >  1:752538/1‑38 (MQ=38)                                                                                                               ggCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                                    <  2:752538/38‑1 (MQ=38)                                                                                                               ggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCt                                                                                                 >  2:441435/1‑57 (MQ=255)                                                                                                               ggCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCt                                                                                                 <  1:441435/57‑1 (MQ=255)                                                                                                                   tatGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCAGAGACAACATCTGGCAAGGc           <  1:803348/137‑1 (MQ=255)                                                                                                                     tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTg                                                                                    >  1:222594/1‑64 (MQ=255)                                                                                                                     tGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTg                                                                                    <  2:222594/64‑1 (MQ=255)                                                                                                                          tCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACg    <  1:486640/139‑1 (MQ=255)                                                                                                                            gTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGa  >  1:955706/1‑139 (MQ=255)                                                                                                                             tGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                      >  2:774644/1‑70 (MQ=255)                                                                                                                             tGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                                      <  1:774644/70‑1 (MQ=255)                                                                                                                                    |                                                                                                                                   CGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGA  >  NC_000913/4518883‑4519144

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