Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,508,454	Δ1 bp	coding (137/879 nt)	php ←	phosphotriesterase homology protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,451	0	C	.	100.0%	331.7 / NA	79	coding (140/879 nt)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (41/38);  total (41/38)

GTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAA  >  NC_000913/3508317‑3508585                                                                                                                                       |                                                                                                                                       gTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccgg                                                                                                                                    >  1:577618/1‑137 (MQ=255)        gAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtcatca                                                                                                                            >  2:218933/1‑139 (MQ=255)           aaaCGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGt                                                                                                                         <  2:322532/139‑1 (MQ=255)                           taaCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCa                                                                                                         <  1:211477/139‑1 (MQ=255)                               cGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                     >  1:1332402/1‑139 (MQ=255)                                  tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTc                                                                                                                    <  1:1081854/121‑1 (MQ=255)                                  tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTc                                                                                                                    >  2:1081854/1‑121 (MQ=255)                                  tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                  <  1:218933/139‑1 (MQ=255)                                      ggCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCAt                                                                                              >  2:622123/1‑139 (MQ=255)                                       gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       <  1:1287747/129‑1 (MQ=255)                                       gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       <  2:1006460/129‑1 (MQ=255)                                       gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       >  1:1006460/1‑129 (MQ=255)                                       gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       >  2:1287747/1‑129 (MQ=255)                                       gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                             >  2:14425/1‑139 (MQ=255)                                                 tggtCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGt                                                                                                                         <  1:758633/98‑1 (MQ=255)                                                 tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGt                                                                                                                         >  2:758633/1‑101 (MQ=255)                                                 tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                   >  1:30341/1‑139 (MQ=255)                                                   aTCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcg                                                                                 <  2:30341/139‑1 (MQ=255)                                                          tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                          >  1:928739/1‑91 (MQ=255)                                                          tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                          <  2:928739/91‑1 (MQ=255)                                                          tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                  >  1:707064/1‑115 (MQ=255)                                                          tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                  <  2:707064/115‑1 (MQ=255)                                                             cgcgCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                       <  1:622123/139‑1 (MQ=255)                                                              gcgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                      <  2:577618/139‑1 (MQ=255)                                                               cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                             <  2:1203442/115‑1 (MQ=255)                                                               cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                             >  1:1203442/1‑115 (MQ=255)                                                               cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgtt                                                                     >  2:20655/1‑139 (MQ=255)                                                               cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgtt                                                                     >  2:355867/1‑139 (MQ=255)                                                                   ttACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       <  2:1187327/101‑1 (MQ=255)                                                                   ttACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       >  1:1187327/1‑101 (MQ=255)                                                                            gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                        <  2:1332402/139‑1 (MQ=255)                                                                            gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                        >  2:1070082/1‑139 (MQ=255)                                                                                aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       <  2:993765/88‑1 (MQ=255)                                                                                aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       >  1:993765/1‑88 (MQ=255)                                                                                aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                    <  1:14425/139‑1 (MQ=255)                                                                                 aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagag                                                   <  1:20655/139‑1 (MQ=255)                                                                                   ttgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                                 <  1:1070082/139‑1 (MQ=255)                                                                                      cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccg                                                                                                                                     >  2:386240/1‑52 (MQ=255)                                                                                      cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccg                                                                                                                                     <  1:386240/53‑2 (MQ=255)                                                                                      cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                                 >  1:1282163/1‑136 (MQ=255)                                                                                      cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                                 <  2:1282163/136‑1 (MQ=255)                                                                                       gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                        >  2:397786/1‑128 (MQ=255)                                                                                       gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                        <  1:397786/128‑1 (MQ=255)                                                                                        cattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       >  2:1340571/1‑80 (MQ=255)                                                                                        cattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                       <  1:1340571/80‑1 (MQ=255)                                                                                           tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                             <  1:1369301/87‑1 (MQ=255)                                                                                           tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                             >  2:1369301/1‑87 (MQ=255)                                                                                           tgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGc                                         >  1:620573/1‑139 (MQ=255)                                                                                                    tGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCa                                >  1:369482/1‑139 (MQ=255)                                                                                                     gTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                      <  2:1351117/100‑1 (MQ=255)                                                                                                     gTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                      >  1:1351117/1‑100 (MQ=255)                                                                                                         cGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTc                                                                                                                        >  2:1336254/1‑46 (MQ=39)                                                                                                         cGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTc                                                                                                                        <  1:1336254/46‑1 (MQ=39)                                                                                                          ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                          <  2:369482/139‑1 (MQ=255)                                                                                                            ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAg                        >  1:843397/1‑139 (MQ=255)                                                                                                             tGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                           >  2:1341539/1‑55 (MQ=255)                                                                                                             tGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                           <  1:1341539/55‑1 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCTTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                      >  1:84164/1‑139 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAAt                                                                                                      <  1:368314/59‑1 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAAt                                                                                                      >  2:368314/1‑59 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgtt                                                                     >  2:947575/1‑92 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgtt                                                                     <  1:947575/92‑1 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                    >  1:884061/1‑109 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                    <  2:884061/109‑1 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCa                                               >  2:1181241/1‑114 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCa                                               <  1:1181241/114‑1 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                          >  2:378101/1‑135 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                          <  1:378101/135‑1 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                      >  1:172616/1‑139 (MQ=255)                                                                                                              ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                      >  2:517060/1‑139 (MQ=255)                                                                                                               gTCATCTCAATCACATTACGCACG‑CCCGGGtcatc                                                                                                                             >  2:592092/1‑35 (MQ=37)                                                                                                               gTCATCTCAATCACATTACGCACG‑CCCGGGtcatc                                                                                                                             <  1:592092/35‑1 (MQ=37)                                                                                                                 cATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAGCGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCaa                                              >  2:82656/1‑112 (MQ=255)                                                                                                                 cATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAGCGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCaa                                              <  1:82656/112‑1 (MQ=255)                                                                                                                       aaTCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCg             <  2:620573/139‑1 (MQ=255)                                                                                                                         tCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTc           <  1:1074947/139‑1 (MQ=255)                                                                                                                         tCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTc           >  2:1074947/1‑139 (MQ=255)                                                                                                                              ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATc      >  2:353301/1‑139 (MQ=255)                                                                                                                               tACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                   <  2:166971/61‑1 (MQ=255)                                                                                                                               tACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                                   >  1:166971/1‑61 (MQ=255)                                                                                                                               tACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                                 <  1:1358453/95‑1 (MQ=255)                                                                                                                               tACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                                 >  2:1358453/1‑95 (MQ=255)                                                                                                                                acgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa    <  1:355867/139‑1 (MQ=255)                                                                                                                                  gcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaaaa  <  1:839997/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       GTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAA  >  NC_000913/3508317‑3508585

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