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

TGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508316‑3508583                                                                                                                                        |                                                                                                                                     tGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGTGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATCacgcacgctcg                                                                                                                                   <  1:967662/139‑4 (MQ=255)    tCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccg                                                                                                                                   >  2:1623738/1‑135 (MQ=255)    tCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccg                                                                                                                                   <  1:1623738/136‑2 (MQ=255)     ccGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtc                                                                                                                              >  1:832464/1‑139 (MQ=255)           aaaaCGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                        <  1:1009741/139‑1 (MQ=255)           aaaaCGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                        <  2:411812/139‑1 (MQ=255)                              aCCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  2:940703/1‑139 (MQ=255)                                cGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   <  1:646516/139‑1 (MQ=255)                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 <  2:987996/139‑1 (MQ=255)                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 <  2:312996/139‑1 (MQ=255)                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGaa                                                                                                 >  1:650862/1‑139 (MQ=255)                                  gTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTCGCAAATGaa                                                                                                 >  2:1068813/1‑139 (MQ=255)                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATGTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  1:1142429/1‑139 (MQ=255)                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGatc                                                                                                >  2:952115/1‑137 (MQ=255)                                   tACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:1142429/139‑1 (MQ=255)                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:1276421/134‑1 (MQ=255)                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  1:1276421/1‑134 (MQ=255)                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           >  1:1422017/1‑139 (MQ=255)                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           >  2:361385/1‑139 (MQ=255)                                        gCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATa                                                                                           >  1:807558/1‑139 (MQ=255)                                           accacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTg                                                                                        <  1:361385/139‑1 (MQ=255)                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                        >  1:301718/1‑105 (MQ=255)                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                        <  2:301718/105‑1 (MQ=255)                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCa                                                                                                                 >  1:721468/1‑112 (MQ=255)                                             cacGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCa                                                                                                                 <  2:721468/112‑1 (MQ=255)                                                gTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCaa                                                                                   <  1:1068813/139‑1 (MQ=255)                                                gTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCaa                                                                                   <  2:832464/139‑1 (MQ=255)                                                    aTCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcg                                                                               >  1:255919/1‑139 (MQ=255)                                                    aTCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcg                                                                               >  2:584349/1‑139 (MQ=255)                                                          gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                        <  2:1330871/92‑1 (MQ=255)                                                          gTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAgg                                                                                                                        >  1:1330871/1‑92 (MQ=255)                                                           tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcgg                                                                              <  1:992026/133‑1 (MQ=255)                                                           tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcgg                                                                              >  2:992026/1‑133 (MQ=255)                                                           tctcGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTcc                                                                        >  2:295814/1‑139 (MQ=255)                                                                cgcATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACgtt                                                                   >  2:658077/1‑139 (MQ=255)                                                                          cAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCAGATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                         <  1:153079/139‑1 (MQ=255)                                                                          cAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGATCGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAg                                                         <  1:952115/139‑1 (MQ=255)                                                                            aGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGcc                                                       <  1:940703/139‑1 (MQ=255)                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  1:302154/1‑92 (MQ=255)                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  2:302154/92‑1 (MQ=255)                                                                             gCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCg                                                      >  2:131725/1‑139 (MQ=255)                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccgg                                                                                                                                  >  2:982120/1‑60 (MQ=255)                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGcccgg                                                                                                                                  <  1:982120/62‑3 (MQ=255)                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    >  1:949121/1‑139 (MQ=255)                                                                               aTAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    <  1:1146415/139‑1 (MQ=255)                                                                                 aaaTTGCGCATTGCGCCCCATTTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  2:949121/139‑1 (MQ=255)                                                                                  aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagag                                                 >  1:95573/1‑139 (MQ=255)                                                                                  aaTTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgagag                                                 >  2:1230032/1‑139 (MQ=255)                                                                                    ttgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               <  1:295814/139‑1 (MQ=255)                                                                                    ttgcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               <  1:1230032/139‑1 (MQ=255)                                                                                      gcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTc                                                                                                                  >  2:1553630/1‑70 (MQ=255)                                                                                      gcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTc                                                                                                                  <  1:1553630/70‑1 (MQ=255)                                                                                      gcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  1:1016616/134‑1 (MQ=255)                                                                                      gcgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  2:1016616/1‑134 (MQ=255)                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtca                                                                                                                             <  2:1052880/58‑1 (MQ=255)                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGtca                                                                                                                             >  1:1052880/1‑58 (MQ=255)                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTGATCTCCTGGCAAATg                                                                                                   >  1:140676/1‑84 (MQ=255)                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATg                                                                                                   <  2:140676/84‑1 (MQ=255)                                                                                       cgcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCaa                                            >  1:592149/1‑139 (MQ=255)                                                                                        gcattgcgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAat                                           >  1:549421/1‑139 (MQ=255)                                                                                           ttgcgcCCCATGTAACGGTTGGTCATCTCAACCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               >  2:1417499/1‑132 (MQ=255)                                                                                           ttgcgcCCCATGTAACGGTTGGTCATCTCAACCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGAt                                               <  1:1417499/132‑1 (MQ=255)                                                                                              cgcCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAg                                     <  1:584349/139‑1 (MQ=255)                                                                                                 cccATGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATg                                  <  2:95573/139‑1 (MQ=255)                                                                                                    aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  2:1563531/1‑120 (MQ=255)                                                                                                    aTGTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  1:1563531/120‑1 (MQ=255)                                                                                                      gTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAat                                           >  2:301427/1‑125 (MQ=255)                                                                                                      gTAACGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCAGCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAat                                           <  1:301427/125‑1 (MQ=255)                                                                                                        aaCGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCgg                                                     <  2:1629554/113‑1 (MQ=255)                                                                                                        aaCGGTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCgg                                                     >  1:1629554/1‑113 (MQ=255)                                                                                                           ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  1:470850/113‑1 (MQ=255)                                                                                                           ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  2:470850/1‑113 (MQ=255)                                                                                                           ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  2:1224553/139‑1 (MQ=255)                                                                                                           ggTTGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCg                        <  1:131725/139‑1 (MQ=255)                                                                                                              tGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAgg                     <  2:1617586/139‑1 (MQ=255)                                                                                                              tGGTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAgg                     >  1:1617586/1‑139 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGtt                                                                                                                   <  1:357922/44‑1 (MQ=38)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGtt                                                                                                                   >  2:357922/1‑44 (MQ=38)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      >  1:1564441/1‑73 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGAt                                                                                      <  2:1564441/73‑1 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcggc                                                                             >  2:1597479/1‑82 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAggcggc                                                                             <  1:1597479/82‑1 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    >  2:1515212/1‑107 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGga                                                    <  1:1515212/107‑1 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                    >  1:742166/1‑139 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                    >  1:1513034/1‑139 (MQ=255)                                                                                                               ggTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAggg                    >  2:524330/1‑139 (MQ=255)                                                                                                                gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  >  2:442601/1‑108 (MQ=255)                                                                                                                gTCATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGgaga                                                  <  1:442601/108‑1 (MQ=255)                                                                                                                  cATCTCAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGgtgt                 <  2:1422017/139‑1 (MQ=255)                                                                                                                     ctcAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                <  2:730086/57‑1 (MQ=255)                                                                                                                     ctcAATCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAAc                                                                                                >  1:730086/1‑57 (MQ=255)                                                                                                                         aTCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGt          >  1:1502534/1‑139 (MQ=255)                                                                                                                         aTCACATTACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGt          <  2:807558/139‑1 (MQ=255)                                                                                                                               ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     <  1:194520/42‑1 (MQ=255)                                                                                                                               ttACGCACG‑CCCGGGTCATCAGGTCGTTCATCTCCTGGCaaa                                                                                                     >  2:194520/1‑42 (MQ=255)                                                                                                                                 acgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:549421/139‑1 (MQ=255)                                                                                                                                 acgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:1502534/139‑1 (MQ=255)                                                                                                                                 acgcacg‑cCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCaa  <  2:255919/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     TGTTCCGGGAAAAACGCGTCCTGGTAATAACCGGTACAGGCCACCACGTTGATCCCCGTCTCGCGCATTACATCAAGCATAAATTGCGCATTGCGCCCCATGTAACGGTTGGTCATCTCAATCACATTACGCACGCCCCGGGTCATCAGGTCGTTCATCTCCTGGCAAATGAACGCATACTGATCAAGGCGGCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAA  >  NC_000913/3508316‑3508583

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