Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,626,825	C→T	V186I (GTT→ATT)	dcp ←	dipeptidyl carboxypeptidase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,626,825	0	C	T	98.7%	247.3 / NA	76	V186I (GTT→ATT)	dcp	dipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/43);  minor base G (0/1);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCAAGGACAAAACGTT  >  NC_000913/1626694‑1626959                                                                                                                                    |                                                                                                                                        cgcCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccg                                                                                                                                  <  2:571472/139‑1 (MQ=255)  gcCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccgc                                                                                                                                 >  1:1111429/1‑139 (MQ=255)      gttgttGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGa                                                                                                                             <  2:375466/139‑1 (MQ=255)             ggtggtATTCAGCAGCGGAGTCAGCCATTTGTTACCCAGACCTTTCTCGCGAGCCGCCTCTGCCGACAGCGAAATCTCTTGCTCACTCATGCCTGCGAGCTGCGCGATATTGTTCACAAGCGGACGGACGGatttattt                                                                                                                      <  1:334032/139‑1 (MQ=255)                       agcagcGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTa                                                                                                            >  2:863740/1‑139 (MQ=255)                          agcGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATc                                                                                                         <  2:197027/139‑1 (MQ=255)                           gcGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                                        <  1:555081/139‑1 (MQ=255)                            cGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGc                                                                                                       <  1:1526406/139‑1 (MQ=255)                            cGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGc                                                                                                       >  1:294361/1‑139 (MQ=255)                                        tttGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActgg                                                                                           <  2:177812/139‑1 (MQ=255)                                        tttGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTCCCAGCTGCGCGATATCGTTCACAATCGGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActgg                                                                                           <  1:1763996/139‑1 (MQ=255)                                                   gACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑t                                                                               >  2:413831/1‑139 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGt                                                                                                             >  1:86320/1‑109 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGt                                                                                                             <  2:86320/109‑1 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATc                                                                                                         <  1:1788332/113‑1 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATc                                                                                                         >  2:1788332/1‑113 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                                        >  1:1716530/1‑114 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                                        <  2:1716530/114‑1 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActg                                                                                            <  2:377987/126‑1 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActg                                                                                            >  1:377987/1‑126 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑Tc                                                                              >  1:131367/1‑139 (MQ=255)                                                    aCCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑Tc                                                                              >  2:1657901/1‑139 (MQ=255)                                                      cTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGc                                                                            <  2:294361/139‑1 (MQ=255)                                                      cTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGc                                                                            >  2:680577/1‑139 (MQ=255)                                                      cTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGc                                                                            <  2:1490373/139‑1 (MQ=255)                                                       tttCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCa                                                                           >  2:1321841/1‑139 (MQ=255)                                                        ttCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTaa                                                                                                >  1:597905/1‑118 (MQ=255)                                                        ttCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTaa                                                                                                <  2:597905/118‑1 (MQ=255)                                                        ttCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          <  1:824726/139‑1 (MQ=255)                                                        ttCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          <  1:1657901/139‑1 (MQ=255)                                                         tctcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccgcc                                                                                                                                <  2:336012/85‑1 (MQ=255)                                                         tctcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAccgcc                                                                                                                                >  1:336012/1‑85 (MQ=255)                                                         tctcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGc                                                                         >  2:836008/1‑139 (MQ=255)                                                         tctcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGc                                                                         >  1:612048/1‑139 (MQ=255)                                                           tcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCtt                                                                       <  2:612048/139‑1 (MQ=255)                                                           tcGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCtt                                                                       <  2:131367/139‑1 (MQ=255)                                                              cgAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAc                                                                                                                                    >  2:267521/1‑76 (MQ=255)                                                              cgAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGAc                                                                                                                                    <  1:267521/76‑1 (MQ=255)                                                                aGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                                        >  2:1594645/1‑102 (MQ=255)                                                                aGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCg                                                                                                        <  1:1594645/102‑1 (MQ=255)                                                                  ccgccTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGAATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTAtt                                                                >  2:1157487/1‑139 (MQ=255)                                                                   cgccTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTc                                                               >  1:116988/1‑139 (MQ=255)                                                                    gccTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCa                                                              <  1:101400/139‑1 (MQ=255)                                                                    gccTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCa                                                              <  1:680577/139‑1 (MQ=255)                                                                          gccgccAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACtt                                                        <  2:660110/139‑1 (MQ=255)                                                                               cAGCGCAATCTCTTGCGCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAAt                                                   <  1:456992/139‑1 (MQ=255)                                                                                aGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTTCGCAGCTTCTGTATTCAGTACTTTTAAt                                                   >  1:1856303/1‑139 (MQ=255)                                                                                aGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCt                                                                     >  2:304268/1‑120 (MQ=255)                                                                                aGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCt                                                                     <  1:304268/120‑1 (MQ=255)                                                                                   gcAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAAt                                                   >  1:681325/1‑135 (MQ=255)                                                                                   gcAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAAt                                                   <  2:681325/135‑1 (MQ=255)                                                                                   gcAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTg                                               <  1:863740/139‑1 (MQ=255)                                                                                        ctctTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTa                                          <  2:1328586/139‑1 (MQ=255)                                                                                         tctTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTAt                                         <  1:1321841/139‑1 (MQ=255)                                                                                                 actcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTg                                               >  1:103389/1‑125 (MQ=255)                                                                                                 actcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTg                                               <  2:103389/125‑1 (MQ=255)                                                                                                 actcaTTCCTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTgc                                 <  2:473398/139‑1 (MQ=255)                                                                                                        ccTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGt                                                                                      <  2:691949/80‑1 (MQ=255)                                                                                                        ccTGCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGt                                                                                      >  1:691949/1‑80 (MQ=255)                                                                                                           gCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAAct                                                                                             <  1:1457401/70‑1 (MQ=255)                                                                                                           gCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAAct                                                                                             >  2:1457401/1‑70 (MQ=255)                                                                                                           gCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          >  1:483850/1‑88 (MQ=255)                                                                                                           gCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          <  2:483850/88‑1 (MQ=255)                                                                                                           gCCAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTgg                       <  2:116988/139‑1 (MQ=255)                                                                                                             cAGCTGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCt                     <  1:836008/139‑1 (MQ=255)                                                                                                                 tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActgg                                                                                           <  1:1120729/66‑1 (MQ=255)                                                                                                                 tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAActgg                                                                                           >  2:1120729/1‑66 (MQ=255)                                                                                                                 tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGc                                                                            >  1:1490654/1‑80 (MQ=255)                                                                                                                 tGCGCGATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGc                                                                            <  2:1490654/80‑1 (MQ=255)                                                                                                                   cgcgATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          >  2:1338771/1‑80 (MQ=255)                                                                                                                   cgcgATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          >  1:1597804/1‑80 (MQ=255)                                                                                                                   cgcgATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          <  2:1597804/80‑1 (MQ=255)                                                                                                                   cgcgATATCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAg                                                                          <  1:1338771/80‑1 (MQ=255)                                                                                                                          tCGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGGTTCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCAAGGACa         <  2:1856303/139‑1 (MQ=255)                                                                                                                           cGTTCACAATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCAAGGACaaa       <  1:413831/139‑1 (MQ=255)                                                                                                                                acaATCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCAAGGACAAAACGtt  >  2:606547/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                        CGCCGGTTGTTGGGTGGTATTCAGCAGCGGAATCAGCCATTTGTTATCCAGACCTTTCTCGCGAGCCGCCTCTGCCGCCAGCGCAATCTCTTGCTCACTCATTCCTGCCAGCTGCGCGATATCGTTCACAACCAGACCGCCGGATTTATTTGCTGCCAGTAATCGCTGGTTAAACTGGCTGGTCAGGG‑TCGCAGCTTCTGTATTCAGTACTTTTAATTTTGCTTTATCAGCTTGCGCAAGTTTGGCTCCGGCAAGGACAAAACGTT  >  NC_000913/1626694‑1626959

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: ‑