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%	355.6 / NA	87	coding (140/879 nt)	php	phosphotriesterase homology protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (43/44);  total (43/44)

AAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATC  >  NC_000913/3508326‑3508581                                                                                                                              |                                                                                                                                   aaaaCGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                      <  1:1399296/139‑1 (MQ=255) aaaaCGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                      <  2:551477/139‑1 (MQ=255)       cgTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTc                                                                                                                <  1:1479550/139‑1 (MQ=255)                aataaCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                       <  2:693702/139‑1 (MQ=255)                aataaCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                       <  2:257038/139‑1 (MQ=255)                    aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   >  2:79570/1‑139 (MQ=255)                    aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   <  1:79570/139‑1 (MQ=255)                      cGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGtt                                                                                                                 >  1:40319/1‑123 (MQ=255)                      cGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGtt                                                                                                                 <  2:40319/123‑1 (MQ=255)                        gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGt                                                                                                                     <  1:643106/117‑1 (MQ=255)                        gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGt                                                                                                                     >  2:643106/1‑117 (MQ=255)                         tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              >  1:729796/1‑139 (MQ=255)                              gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtcatca                                                                                                                        >  1:1055293/1‑108 (MQ=255)                              gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtcatca                                                                                                                        <  2:1055293/108‑1 (MQ=255)                              gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                         >  1:1381278/1‑139 (MQ=255)                                caccacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACt                                                                                       >  1:698706/1‑139 (MQ=255)                                   cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                    >  1:1114901/1‑139 (MQ=255)                                   cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                    >  1:981801/1‑139 (MQ=255)                                        tGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                               >  1:1255930/1‑139 (MQ=255)                                         gATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggc                                                                              >  1:604282/1‑139 (MQ=255)                                             cccGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCa                                                                          >  2:130153/1‑139 (MQ=255)                                                gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTc                                                                       <  1:775039/139‑1 (MQ=255)                                                           tACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                            <  2:1114901/139‑1 (MQ=255)                                                           tACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGtttt                                                            <  2:1255930/139‑1 (MQ=255)                                                              aTCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTaa                                                         <  2:981801/139‑1 (MQ=255)                                                                 aaGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTaaa                                                        <  1:1467002/137‑1 (MQ=255)                                                                 aaGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTaaa                                                        >  2:1467002/1‑137 (MQ=255)                                                                   gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   >  1:754512/1‑92 (MQ=255)                                                                   gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                   <  2:754512/92‑1 (MQ=255)                                                                     aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTcc                                                                                                          <  2:102935/83‑1 (MQ=255)                                                                     aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTcc                                                                                                          >  1:102935/1‑83 (MQ=255)                                                                     aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCa                                                                                           <  2:496694/98‑1 (MQ=255)                                                                     aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCa                                                                                           >  1:496694/1‑98 (MQ=255)                                                                       aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAt                                                                                                              >  1:49328/1‑77 (MQ=255)                                                                       aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAt                                                                                                              <  2:49328/77‑1 (MQ=255)                                                                       aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  1:770858/90‑1 (MQ=255)                                                                       aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 >  2:770858/1‑90 (MQ=255)                                                                       aaaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                >  2:520513/1‑139 (MQ=255)                                                                        aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagag                                               <  2:604282/139‑1 (MQ=255)                                                                         attgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                  >  2:1129176/1‑56 (MQ=255)                                                                         attgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGccc                                                                                                                                  <  1:1129176/56‑1 (MQ=255)                                                                          ttgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                             <  2:355214/139‑1 (MQ=255)                                                                             cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                   <  1:30978/114‑1 (MQ=255)                                                                             cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACg                                                                   >  2:30978/1‑114 (MQ=255)                                                                             cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                >  1:349099/1‑133 (MQ=255)                                                                             cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                <  2:349099/133‑1 (MQ=255)                                                                              gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                       >  2:1265113/1‑125 (MQ=255)                                                                              gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                       <  1:1265113/125‑1 (MQ=255)                                                                               cattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAata                                        <  2:729796/139‑1 (MQ=255)                                                                                    cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                               >  2:559889/1‑79 (MQ=255)                                                                                    cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                               <  1:559889/79‑1 (MQ=255)                                                                                    agacagaTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCTGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAg                                   <  1:1486298/133‑1 (MQ=255)                                                                                         cATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCt                              >  1:1705128/1‑139 (MQ=255)                                                                                          aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCaa                                                                                 <  2:1290420/87‑1 (MQ=255)                                                                                          aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCaa                                                                                 >  1:1290420/1‑87 (MQ=255)                                                                                            gTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCAt                           <  2:698706/139‑1 (MQ=255)                                                                                             tAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATg                          <  1:130153/139‑1 (MQ=255)                                                                                              aaCGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                               <  2:1058322/69‑1 (MQ=255)                                                                                              aaCGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                               >  1:1058322/1‑69 (MQ=255)                                                                                                cGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                               >  2:12886/1‑83 (MQ=255)                                                                                                cGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAgg                                                                               <  1:12886/83‑1 (MQ=255)                                                                                                 ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              >  1:209024/1‑67 (MQ=255)                                                                                                 ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                              <  2:209024/67‑1 (MQ=255)                                                                                                   ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCa                                                                          <  2:1207485/85‑1 (MQ=255)                                                                                                   ttGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCa                                                                          >  1:1207485/1‑85 (MQ=255)                                                                                                     ggacATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTAGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                  >  2:1004952/4‑91 (MQ=255)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                      >  2:623146/1‑39 (MQ=37)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                      <  1:623146/39‑1 (MQ=37)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                       >  2:1686506/1‑54 (MQ=39)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTgg                                                                                                       <  1:1686506/54‑1 (MQ=39)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 >  2:1047578/1‑60 (MQ=255)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  1:1047578/60‑1 (MQ=255)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCCGCCCACGTTGTTTTTAAAGCCGAGGAGATCAATATGCAGATGCTCATGGGctgtct                  >  2:1486298/1‑134 (MQ=255)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                  <  1:1004952/91‑1 (MQ=255)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAata                                        <  2:394519/117‑1 (MQ=255)                                                                                                     ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAata                                        >  1:394519/1‑117 (MQ=255)                                                                                                      gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                             <  2:1549482/111‑1 (MQ=255)                                                                                                      gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                             >  1:1549482/1‑111 (MQ=255)                                                                                                      gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGgt                 >  2:1493114/1‑139 (MQ=255)                                                                                                               aTCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 >  2:649404/1‑50 (MQ=255)                                                                                                               aTCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                 <  1:649404/50‑1 (MQ=255)                                                                                                                tCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTc       >  2:962330/1‑139 (MQ=255)                                                                                                                     ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                  >  2:1374006/1‑75 (MQ=255)                                                                                                                     ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgt                                                                  <  1:1374006/75‑1 (MQ=255)                                                                                                                     ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCa                            >  2:680642/1‑113 (MQ=255)                                                                                                                     ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCa                            <  1:680642/113‑1 (MQ=255)                                                                                                                     ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                      >  1:788741/1‑119 (MQ=255)                                                                                                                     ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                      <  2:788741/119‑1 (MQ=255)                                                                                                                     ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATc  >  1:645456/1‑139 (MQ=255)                                                                                                                              |                                                                                                                                   AAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATC  >  NC_000913/3508326‑3508581

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