Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,761,230	G→A	S7F (TCC→TTC)	selA ←	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,761,230	0	G	A	97.4%	128.2 / NA	38	S7F (TCC→TTC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/14);  minor base C (0/1);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGC  >  NC_000913/3760984‑3761475                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      gTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAAcg                                                                                                                                                                                                                                                   >  1:218229/1‑251 (MQ=255) gTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAAcg                                                                                                                                                                                                                                                   >  1:55249/1‑251 (MQ=255)                         cGCAGCGCGCTCTGCTCTTAGTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGGCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCt                                                                                                                                                                                                                           <  2:249455/251‑1 (MQ=255)                             gcgcgcTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCt                                                                                                                                                                                                                       <  1:99840/251‑1 (MQ=255)                                      gctcttctgTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAg                                                                                                                                                                                                                <  1:256222/246‑1 (MQ=255)                                      gcgcTTCTGTCGGCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGTGCTTCGTCGAGCATCTGACGCAACAATTCCCCCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGAGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTTTCATAGGAAACCTGGCTGATCaa                                                                                                                                                                                                                 <  1:160521/248‑1 (MQ=255)                                                  tCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACgg                                                                                                                                                                                                  >  1:59183/1‑251 (MQ=255)                                                        gggCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAgg                                                                                                                                                                                            <  2:4218/251‑1 (MQ=255)                                                                          cccAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATc                                                                                                                                                                                                                   >  2:74712/1‑210 (MQ=255)                                                                          cccAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATc                                                                                                                                                                                                                   <  1:74712/210‑1 (MQ=255)                                                                           ccAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACgg                                                                                                                                                                         >  2:101215/1‑251 (MQ=255)                                                                                      accacGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTg                                                                                                                                                              <  1:245007/251‑1 (MQ=255)                                                                                            caggcagCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATgg                                                                                                                                                        >  1:27022/1‑251 (MQ=255)                                                                                                 agCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCag                                                                                                                                                   <  1:24738/251‑1 (MQ=255)                                                                                                          cTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACat                                                                                                                                           >  1:260728/1‑250 (MQ=255)                                                                                                                    tCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATc                                                                                                                                 >  1:257539/1‑250 (MQ=255)                                                                                                                           tcGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCttt                                                                                                                         >  1:425/1‑251 (MQ=255)                                                                                                                            cgcgCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCttt                                                                                                                         >  1:246239/1‑250 (MQ=255)                                                                                                                                              tGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCa                                                                                                       >  1:80412/1‑250 (MQ=255)                                                                                                                                               gACGCAACAATTCCACCACGCGGGTGGGACCATAAGTAGCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAAGAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGGTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATGTCGCGCATCAAGTCTTTGGCGGTTCGGTGCgggc                                                                                                        <  2:249544/248‑4 (MQ=255)                                                                                                                                                 cGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCTCGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATAGCGCGCATCAAGCCTTTGGCGGTTCGTTGCGAGCAAAAc                                                                                                   <  2:11604/251‑1 (MQ=255)                                                                                                                                                      caaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAATAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATATCCGCGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGGGTTTGGCGGTTTTGTTCGAGCAAAATTTTcg                                                                                              <  2:240155/251‑1 (MQ=255)                                                                                                                                                       aaTTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTcgc                                                                                             >  2:14277/1‑251 (MQ=255)                                                                                                                                                           ccaccaCGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGAGTGTAGCGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCt                                                                                         <  1:257075/251‑1 (MQ=255)                                                                                                                                                                                        cGCAAAGAAAGGAAGGAGCTATCGCGCACTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGGGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAgg                                                            >  2:255782/1‑250 (MQ=255)                                                                                                                                                                                               aaagGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGGTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCg                                                     >  2:72452/1‑251 (MQ=255)                                                                                                                                                                                                aaggaaggAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCgg                                                    >  1:39435/1‑251 (MQ=255)                                                                                                                                                                                                aaggaaggAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCgg                                                    >  1:13013/1‑251 (MQ=255)                                                                                                                                                                                                   gaaggaGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCCAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCg                                                 >  2:209227/1‑251 (MQ=255)                                                                                                                                                                                                              cgcgCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAAc                                       <  1:230029/250‑1 (MQ=255)                                                                                                                                                                                                                   aaTAAGCGATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTCCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCCCACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTTGGTGCGAGCAAACATTTCGCGGCTAACCAGGTTTTAGCCCAGTTTGACTCAGTATGGGCGGGCGAcaacccaacccaccc                                 >  2:145557/1‑246 (MQ=255)                                                                                                                                                                                                                          gATCAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAAcgcg                          >  2:77303/1‑251 (MQ=255)                                                                                                                                                                                                                             cAATAGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACgcggc                        <  1:244797/250‑1 (MQ=255)                                                                                                                                                                                                                                 aGCCGGAAGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGaa                    >  1:220078/1‑250 (MQ=255)                                                                                                                                                                                                                                    cGGAAGTTGACTATAGAGGCAACTCGGTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCGCTCACACGGAGAAGGGCGTTTAACATAGCCACGGATTGTAACGTGAGATGGTTCAGGAGGACATATCGCGGATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTGCGCGGCTAAACAGGTTTTCGGCGATTTTGATTAAGTGTGGGCTGTCGACACACCAACCTGGCGCAACGCGGCCGAAATTg                <  2:155819/251‑1 (MQ=255)                                                                                                                                                                                                                                        aGTTGACTATAGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTAAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGGGGGGGAAAAACTTTCGCGGATAAACAGGTATTCGAACAGTTTGACTAAGAGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTgaga             >  2:82633/1‑250 (MQ=255)                                                                                                                                                                                                                                               tataGAGGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGAc     >  1:135197/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gagGAAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGACGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTCAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCCCAACGCTGCGGAAATTGAGATATCCGACGGc  >  2:209964/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      GTTCCCGTCAGGTTGATCACCGGACGCAGCGCGCTCTGCGCTTCTTTCGTCAACCGGGCATCGACTTCTTGCGCCCAGTTTTCACACCACGCAGGCAGCGTCTGGCTGCCACGAATCACTTCTCGCGCTTCGTCGAGCATCTGACGCAACAATTCCACCACGCGGGTGTGACCATAAGTATCACGCAAAGAAAGGAAGGAGCTATCGCGCAATAAGCGATCAATAGCCGGAAGTTGACTATAGAGGGAACGCGTTTCGGTTGTCATAGGAAACCTGGCTGATCAAGGCCCTCTCACACGGAGAAGGGCGTTTAACATAACCACGGATTGTAACGTGAGATGGGTCAGGAGGACATATCGCGCATCAAGCCTTTGGCGGTTCGGTGCGAGCAAAACTTTCGCGGCTAAACAGGTTTTCGACCAGTTTGACTAAGTGTGGGCGGTCGACACACCAACCTGGCGCAACGCGGCGGAAATTGAGATATCCGACGGC  >  NC_000913/3760984‑3761475

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