Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,068,285	Δ1 bp	intergenic (+31/+226)	ymdF → / ← rutG	KGG family protein/pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,068,285	0	T	.	100.0%	322.2 / NA	79	intergenic (+31/+226)	ymdF/rutG	KGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (36/43);  total (36/43)

GGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACT  >  NC_000913/1068156‑1068412                                                                                                                                  |                                                                                                                                ggACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGt                                                                                                                      >  1:716383/1‑139 (MQ=255)            cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCAt                                                                                                           >  1:693065/1‑139 (MQ=255)               gggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAc                                                                                                        >  2:58208/1‑139 (MQ=255)                 gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAg                                                                                                                           <  2:35065/118‑1 (MQ=255)                 gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAg                                                                                                                           >  1:35065/1‑118 (MQ=255)                 gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAccg                                                                                                      >  1:1440424/1‑139 (MQ=255)                  tAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAccgc                                                                                                     >  2:494411/1‑139 (MQ=255)                   aaTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAccgcc                                                                                                    >  2:1715204/1‑139 (MQ=255)                          aGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             >  1:1017296/1‑139 (MQ=255)                            aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  2:464373/139‑1 (MQ=255)                            aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  2:386290/139‑1 (MQ=255)                            aaTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  2:1017296/139‑1 (MQ=255)                              tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggt                                                                                         >  2:606412/1‑139 (MQ=255)                              tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggt                                                                                         >  1:127021/1‑139 (MQ=255)                                 cccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                                      <  1:1073225/139‑1 (MQ=255)                                 cccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                                      <  2:1440424/139‑1 (MQ=255)                                      aGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                                                 <  2:716383/139‑1 (MQ=255)                                              tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         <  2:693065/139‑1 (MQ=255)                                                       ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAg                                                                >  2:5623/1‑139 (MQ=255)                                                       ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAg                                                                >  2:1761398/1‑139 (MQ=255)                                                         tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCt                                                              >  2:833764/1‑139 (MQ=255)                                                                  tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         <  2:865329/119‑1 (MQ=255)                                                                  tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         >  1:865329/1‑119 (MQ=255)                                                                  tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATg                                                     >  2:1641122/1‑139 (MQ=255)                                                                  tggtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATg                                                     >  1:1280341/1‑139 (MQ=255)                                                                    gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   <  1:494411/139‑1 (MQ=255)                                                                    gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   <  1:1641122/139‑1 (MQ=255)                                                                    gtAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   <  2:1280341/139‑1 (MQ=255)                                                                     tAAGAGCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcgc                                                  <  1:5623/139‑1 (MQ=255)                                                                       agagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                                <  1:1715204/139‑1 (MQ=255)                                                                       agagCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                                <  1:58208/139‑1 (MQ=255)                                                                          gCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTg                                                                                    >  2:1707822/1‑100 (MQ=255)                                                                          gCAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTg                                                                                    <  1:1707822/100‑1 (MQ=255)                                                                           cAGTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCt                                            >  1:1047695/1‑139 (MQ=255)                                                                             gTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTa                                          <  2:148504/139‑1 (MQ=255)                                                                             gTCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTa                                          <  2:1047695/139‑1 (MQ=255)                                                                              tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  1:360694/89‑1 (MQ=255)                                                                              tCACGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           >  2:360694/1‑89 (MQ=255)                                                                                aCGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             <  2:1872671/85‑1 (MQ=255)                                                                                aCGGCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             >  1:1872671/1‑85 (MQ=255)                                                                                  ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                     <  1:1761398/139‑1 (MQ=255)                                                                                  ggCAAAAGCGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAgg                                     <  1:1363034/139‑1 (MQ=255)                                                                                          cGACAACTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCag                             <  2:127021/139‑1 (MQ=255)                                                                                              aaCTAACCGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaaga                         <  1:902788/139‑1 (MQ=255)                                                                                                   accGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAccgc                                                                                                     >  2:904455/2‑58 (MQ=255)                                                                                                   accGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCAccgc                                                                                                     <  1:904455/57‑1 (MQ=255)                                                                                                   accGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                                <  2:558690/110‑1 (MQ=255)                                                                                                   accGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                                >  1:558690/2‑111 (MQ=255)                                                                                                   accGGGCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                    <  1:1432898/138‑1 (MQ=255)                                                                                                        gCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAAt                                                                     <  1:1526146/85‑1 (MQ=255)                                                                                                        gCTAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAAt                                                                     >  2:1526146/1‑85 (MQ=255)                                                                                                          tAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTa                                                                                               >  2:117292/1‑57 (MQ=255)                                                                                                          tAATCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTa                                                                                               <  1:117292/57‑1 (MQ=255)                                                                                                           aaTCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             <  1:309203/58‑1 (MQ=255)                                                                                                           aaTCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             >  2:824661/1‑58 (MQ=255)                                                                                                           aaTCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             >  2:309203/1‑58 (MQ=255)                                                                                                           aaTCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             <  1:824661/58‑1 (MQ=255)                                                                                                            aTCAATGACGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         <  2:1085583/77‑1 (MQ=255)                                                                                                            aTCAATGACGAAT‑CATTTTTG‑CTGTAGCTAGTCAAAAAGCCATCAGCGCCGGTTACCCGGTGGTTGATACTGATGAc                                                                         >  1:1085583/1‑77 (MQ=255)                                                                                                                   aCGAAT‑CATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCa    <  1:1656150/139‑1 (MQ=255)                                                                                                                      gaatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACt  <  1:833764/135‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGt                                                                                                 >  2:1053274/4‑43 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGt                                                                                                 <  1:1053274/40‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             <  2:1323744/44‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                             >  1:1323744/4‑47 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           >  2:713951/4‑49 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                           <  1:713951/46‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   <  2:1129816/86‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   <  1:753306/86‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   >  1:1129816/4‑89 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                   >  2:753306/4‑89 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                                >  1:1692814/4‑92 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGa                                                <  2:1692814/89‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCa                              >  1:222609/4‑110 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCa                              <  2:222609/107‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                    <  2:688469/117‑1 (MQ=255)                                                                                                                       aatcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCa                    >  1:688469/4‑120 (MQ=255)                                                                                                                        atcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                         <  2:1053399/96‑1 (MQ=255)                                                                                                                        atcATTTTTG‑CTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTAt                                         >  1:1053399/3‑98 (MQ=255)                                                                                                                                  |                                                                                                                                GGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACT  >  NC_000913/1068156‑1068412

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