Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,702,757	C→T	F175F (TTC→TTT)	add →	adenosine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,702,757	0	C	T	94.6%	111.3 / ‑6.8	38	F175F (TTC→TTT)	add	adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (17/18);  minor base G (0/2);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCA  >  NC_000913/1702517‑1703000                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    cGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCAGGTGATGAACTTGGTTTTCCGGGCAGt                                                                                                                                                                                                                                            >  1:252440/1‑250 (MQ=255)      tGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGGGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGACTGTTc                                                                                                                                                                                                                                      >  2:103407/1‑251 (MQ=255)        cGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCt                                                                                                                                                                                                                                    <  1:145685/251‑1 (MQ=255)              tCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTtctc                                                                                                                                                                                                                               <  1:130277/250‑1 (MQ=255)                            ggcaatggcaCATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACcgcgcgc                                                                                                                                                                                                                >  1:124114/1‑251 (MQ=255)                                atggcaCATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGa                                                                                                                                                                                                            <  1:103407/251‑1 (MQ=255)                                           gCTGCCTGTAGCGGGTGTTGCCGACGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCCGCCAGAGCTGGAGGCCTTTTTAGCCCACCGGGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGc                                                                                                                                                                                                 >  2:164671/1‑251 (MQ=255)                                            cTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGc                                                                                                                                                                                                 >  1:7049/1‑250 (MQ=255)                                              gCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCat                                                                                                                                                                                               <  1:77016/250‑1 (MQ=255)                                               ccTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCatat                                                                                                                                                                                             >  1:214023/1‑251 (MQ=255)                                                  gTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATAtt                                                                                                                                                                                            <  1:164671/249‑1 (MQ=255)                                                   tcgagggTGTTGTCGAAGCGGGGATCGCTGGCGTACGTGAAGGTGGCAGCACCTTTGGGGTGCAGGCGAAGCTTAACGGCATTATGAGCCGGACCTTCGGCGAAGCCGACTGTCAGGAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCT‑CACTTAGTTTAGCAGGTGATTAGTGTGGTTTTCCGGCGAGTCTGTTCCATTCTCACTTCAAAAGCACGAGTGATGCGGGATGGCagtttagc                                                                                                                                                                                         >  2:26265/5‑244 (MQ=255)                                                     gCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGGGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGCTTACCGTCA‑TTGGTTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGCTGCGGGCTGGCATCTTACCGt                                                                                                                                                                                       >  2:285174/1‑251 (MQ=255)                                                         gtgtTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTGCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCAt                                                                                                                                                                                   <  2:145624/251‑1 (MQ=255)                                                          gttggtCGAAGCGGTGATCGATGGGGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTGCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCTCCGTGACCAGATTACCG‑CCCTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCCCTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATg                                                                                                                                                                                  <  2:247424/247‑1 (MQ=255)                                                                            cGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACGGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTgccggg                                                                                                                                                                 <  1:285174/250‑1 (MQ=255)                                                                                    tACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGc                                                                                                                                                        >  1:317962/1‑251 (MQ=255)                                                                                                  cGCATCTTTTGTGGGCAGGCGAAGTTTAGCGGCATTATGAGCCGGACCTTTGGCGAAGCCGCCTGTCAGCAAGAGCTTGAGGCCTGTTTATCCCACCGTGACCAGATTTCCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCGCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCTGAAAGCATCTGGCAGGCGAt                                                                                                                                          <  1:332567/251‑1 (MQ=255)                                                                                                      ccTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGt                                                                                                                                      <  1:277719/251‑1 (MQ=255)                                                                                                         ttGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATGCGTGa                                                                                                                                    <  1:287489/250‑1 (MQ=255)                                                                                                                  aGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAAGCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGACCTGGGTGCg                                                                                                                          >  1:117433/1‑251 (MQ=255)                                                                                                                  aGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCg                                                                                                                          >  1:83866/1‑251 (MQ=255)                                                                                                                          cTTATCGACATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTAt                                                                                                                  >  1:139172/1‑251 (MQ=255)                                                                                                                              ccgtcATTAGGAGCCGGACCTTCTGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGa                                                                                                              <  1:219738/247‑1 (MQ=255)                                                                                                                                   aTTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCGTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATgg                                                                                                         <  1:139751/251‑1 (MQ=255)                                                                                                                                              gACCTTCGGCGAAGCCGCCTGTCAGCAATAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCa                                                                                              <  2:214023/251‑1 (MQ=255)                                                                                                                                               cccTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCAt                                                                                             <  1:26265/250‑1 (MQ=255)                                                                                                                                                                       cAAGAGCTGGAGGCCTTTTTAGCCCAACGTGACCAGATTACCG‑CACTTGATTTAGCCGGAGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAttt                                                                     >  1:129920/1‑251 (MQ=255)                                                                                                                                                                        aagatttGGAGGCTTTTTTAGCCCACCGTTACCAGATTACCG‑CTCTTGATTTAGCCGGGGATGAACTTGGTTTTCCGGGAAGTCGGTTCCTTTCGCACTTCAACCGCGCGCGTGAGGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGAtttt                                                                    <  2:215313/245‑1 (MQ=255)                                                                                                                                                                          gagCTGGAGGCCTTTTTAGCCCCCCGGGACCAGATTACCG‑CACGTGATTTAGCCGGAGATGAACTTGGTTTGCCGTGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGAAGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGTCGATTCGTGAACTGGGTGCGGAGCGTATTTGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTtct                                                                  <  1:211139/251‑1 (MQ=255)                                                                                                                                                                            ggtgtaGGCCTTTTTAGACCACCGTGACCAGATTAGCG‑CATGTTAGTTAGTCGGTTATGAATTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGAGCGCGTGATGCGGGCTGGCAGATTGCCTTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCg                                                                <  2:185337/246‑1 (MQ=255)                                                                                                                                                                                            gCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGCACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGGATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCGGATGGATTTTCTCGCCGAGCAACAAATTgg                                                >  2:263437/1‑251 (MQ=255)                                                                                                                                                                                                aacGTGCCCAGCTTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCATTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGACAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGAATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGTTGGATTTTCTCGCCGAGCAAAAAATTGGAAtt                                            >  2:221244/3‑251 (MQ=255)                                                                                                                                                                                                  cGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTg                                           >  1:159237/1‑250 (MQ=255)                                                                                                                                                                                                   gTGACCAGAGTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGaa                                         <  1:244288/251‑1 (MQ=255)                                                                                                                                                                                                        cAGATTTCCG‑AACCTGAGTTAGCCGGTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCGGCCGGGACGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTAGTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTg                                    <  2:139172/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggTGATGAACTTGGTTTTCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGAcc              >  1:185012/1‑249 (MQ=255)                                                                                                                                                                                                                                            tGGTTTGCCGGGACGTCGGTTCCTTTCTCACGTCAACCGCGCGCGGGATGCGTGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGTCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTTGCa  <  2:83866/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    CGAGCTGCGTTTTTCACCAGGCTACATGGCAATGGCACATCAGCTGCCTGTAGCGGGTGTTGTCGAAGCGGTGATCGATGGCGTACGTGAAGGTTGCCGCACCTTTGGTGTGCAGGCGAAGCTTATCGGCATTATGAGCCGGACCTTCGGCGAAGCCGCCTGTCAGCAAGAGCTGGAGGCCTTTTTAGCCCACCGTGACCAGATTACCG‑CACTTGATTTAGCCGGTGATGAACTTGGTTTCCCGGGAAGTCTGTTCCTTTCTCACTTCAACCGCGCGCGTGATGCGGGCTGGCATATTACCGTCCATGCAGGCGAAGCTGCCGGGCCGGAAAGCATCTGGCAGGCGATTCGTGAACTGGGTGCGGAGCGTATTGGACATGGCGTAAAAGCCATTGAAGATCGGGCGCTGATGGATTTTCTCGCCGAGCAACAAATTGGTATTGAATCCTGTCTGACCTCCAATATTCAGACCAGCACCGTAGCA  >  NC_000913/1702517‑1703000

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