Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	255,659	A→G	C20R (TGT→CGT)	pepD ←	aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	255,659	0	A	G	97.1%	238.9 / ‑3.5	69	C20R (TGT→CGT)	pepD	aminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (35/32);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAC  >  NC_000913/255530‑255794                                                                                                                                  |                                                                                                                                        cGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTCTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGc                                                                                                                                >  2:256872/1‑139 (MQ=255)                   tCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTCTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCgg                                                                                                              <  1:256872/139‑1 (MQ=255)                         tatTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTgg                                                                                                        <  1:1182213/139‑1 (MQ=255)                           ttACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      <  2:646339/139‑1 (MQ=255)                            tacctaccTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          >  2:877218/1‑134 (MQ=255)                            tacctaccTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCtgt                                                                                                          <  1:877218/134‑1 (MQ=255)                            tacctaccTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgaga                                                                                                     >  1:988449/1‑139 (MQ=255)                                 acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                                                                                                >  2:681068/1‑139 (MQ=255)                                 acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                                                                                                >  1:197279/1‑139 (MQ=255)                                 acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                                                                                                >  1:287484/1‑139 (MQ=255)                                 acctGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAAtt                                                                                                >  1:645228/1‑139 (MQ=255)                                     gATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  1:1019850/139‑1 (MQ=255)                                      aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           >  2:1041958/1‑139 (MQ=255)                                      aTCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGAc                                                                                           >  2:1086527/1‑139 (MQ=255)                                       tCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  1:275010/137‑1 (MQ=255)                                       tCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            >  2:275010/1‑137 (MQ=255)                                        cgcgTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAg                                                                                         >  2:455143/1‑139 (MQ=255)                                         gcgTTCGACATGGAGACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCGAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGt                                                                                        >  2:740869/1‑139 (MQ=255)                                            ttCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCa                                                                                     <  2:469732/139‑1 (MQ=255)                                                 cATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGacac                                                                                <  2:287484/139‑1 (MQ=255)                                                     gAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            <  1:502475/123‑1 (MQ=255)                                                     gAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTgaga                                                                                            >  2:502475/1‑123 (MQ=255)                                                            ttCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTcc                                                                     <  2:540103/139‑1 (MQ=255)                                                                     gCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATc                                                            >  2:265575/1‑139 (MQ=255)                                                                       ccaaccaaCAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCAc                                                          <  2:584557/139‑1 (MQ=255)                                                                          accaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                       <  2:197976/139‑1 (MQ=255)                                                                          accaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCg                                                       <  1:856327/139‑1 (MQ=255)                                                                           ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    <  2:240043/125‑1 (MQ=255)                                                                           ccaacAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    >  1:240043/1‑125 (MQ=255)                                                                             aacaaTGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCaa                                                    <  1:1041958/139‑1 (MQ=255)                                                                                     aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAt                                            >  2:892568/1‑139 (MQ=255)                                                                                     aTTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTAt                                            >  2:493797/1‑139 (MQ=255)                                                                                      ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCaa                                                    <  2:415198/130‑1 (MQ=255)                                                                                      ttCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCaa                                                    >  1:415198/1‑130 (MQ=255)                                                                                       tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCaa                                                    >  1:154625/1‑129 (MQ=255)                                                                                       tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTg                                          >  2:619068/1‑139 (MQ=255)                                                                                       tCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGAAACGTTAAGTCTCCTTGTCGATCACCCGCaa                                                    <  2:154625/129‑1 (MQ=255)                                                                                        cAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  1:455143/139‑1 (MQ=255)                                                                                         aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    >  2:407482/1‑111 (MQ=255)                                                                                         aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCt                                                                    <  1:407482/111‑1 (MQ=255)                                                                                         aGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCa                                        >  2:113459/1‑139 (MQ=255)                                                                                          gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt                                                                             <  1:395335/101‑1 (MQ=255)                                                                                          gCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGtt                                                                             >  2:395335/1‑101 (MQ=255)                                                                                             aGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTc                                    <  2:362715/139‑1 (MQ=255)                                                                                              gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  1:596420/133‑1 (MQ=255)                                                                                              gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         >  2:596420/1‑133 (MQ=255)                                                                                              gTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTca                                   >  1:348199/1‑139 (MQ=255)                                                                                               ttGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                  <  2:197279/139‑1 (MQ=255)                                                                                                tGCTCTTCATGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcaca                                 <  2:988449/139‑1 (MQ=255)                                                                                                        aTGATAGGACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGc                         >  1:837424/1‑139 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      <  1:427065/55‑1 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGgag                                                                                                      >  2:427065/1‑55 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                  >  2:888881/1‑139 (MQ=255)                                                                                                               gACGGGTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTATGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAt                  >  1:119352/1‑139 (MQ=255)                                                                                                                  gggTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgtt               >  2:563659/1‑139 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCAc                                                          <  2:17182/94‑1 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCAc                                                          >  1:17182/1‑94 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                  <  1:820684/118‑1 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTcac                                  >  2:820684/1‑118 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt             >  2:787898/1‑139 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt             <  2:1094351/139‑1 (MQ=255)                                                                                                                    gTGAGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAAtgttgt             >  2:898220/1‑139 (MQ=255)                                                                                                                       aGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         >  1:915898/1‑108 (MQ=255)                                                                                                                       aGGAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGc                                         <  2:915898/108‑1 (MQ=255)                                                                                                                         gAATAGAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAaca        <  1:892568/139‑1 (MQ=255)                                                                                                                              gAACGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGa   >  2:282759/1‑139 (MQ=255)                                                                                                                               aaCGGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAc  <  1:619068/139‑1 (MQ=255)                                                                                                                               aaCGGATTTTGGCAAAAATATCCCACAGCGGCTATGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAc  <  2:119352/139‑1 (MQ=255)                                                                                                                               aaCAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAc  <  2:645228/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                                        CGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACAGATTTTGGCAAAAATATCCCACAGCGGCTGTGGAGATAATTGAGACAGTTCAGACACGTTAAGTCTCCTTGTCGATCACCCGCAAAACAGTATTGCAGGTCACAGGGTTAGCAGAAAATGTTGTCAACACAAGAC  >  NC_000913/255530‑255794

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