Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	780,624	C→T	noncoding (71/76 nt)	lysT →	tRNA‑Lys

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	780,624	0	C	T	97.9%	146.2 / NA	48	noncoding (71/76 nt)	lysT	tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/25);  minor base A (1/0);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGT  >  NC_000913/780378‑780862                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               aaCGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACtcacc                                                                                                                                                                                                                                            >  1:408712/1‑251 (MQ=255)       aTGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGCCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGAGGGTCGATAGCTCAGTTGGTAGAGCAGTTGCCTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCCCGACTCACCAATCGc                                                                                                                                                                                                                                      >  2:309253/1‑251 (MQ=255)         gTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCt                                                                                                                                                                                                                                     <  1:474254/250‑1 (MQ=255)         gTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCt                                                                                                                                                                                                                                     >  1:41016/1‑250 (MQ=255)            aTGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAgg                                                                                                                                                                                                                                 <  1:21412/251‑1 (MQ=255)             tGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGAAGAGCAGTTGACTTTTAATCAATTGGTCACAGGTTCGAATCCTGCACGATccacatgttttgacttc                                                                                                                                                                                                                                >  2:68700/1‑238 (MQ=255)                    cacacGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCg                                                                                                                                                                                                                         <  1:252599/251‑1 (MQ=255)                       ccGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTTCATCGAAATTTTGGTTGCTCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGgta                                                                                                                                                                                                                      <  2:166055/250‑1 (MQ=255)                       acGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATCTATGCCGCCGTTGCCACGGGATATCAACCAAACCGAAAGCAACGAAACAGTGGGTCGTTAGCTAAGTTGGTAGCGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccaccagttgtaaagtcggaagcgga                                                                                                                                                                                                                       >  2:302385/1‑230 (MQ=38)                         gACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGgtagt                                                                                                                                                                                                                    >  1:38588/1‑251 (MQ=255)                          aCCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGAGGGTCGTTAGCTCAGTTGGTAGAGCCGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACACACCAATTGCAAAggcgacagcggatga                                                                                                                                                                                                                    >  1:311918/1‑237 (MQ=38)                              gAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAAc                                                                                                                                                                                                               <  1:242500/251‑1 (MQ=255)                               aaGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCATGCACGACTCACCAATCGCTAAGGTGGAAGCGGAGGTAAAACg                                                                                                                                                                                                              >  1:335524/1‑251 (MQ=255)                                   cGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGa                                                                                                                                                                                                           >  2:406073/1‑250 (MQ=255)                                     cATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAgg                                                                                                                                                                                                        >  1:175891/1‑251 (MQ=255)                                      attactTCTGGTCGTCTCGTGCGATTCATAAGGAGTTGAGGGCTCCCCATCGAAATTGTTGTTGGGTACAAAGT‑TGAAACCAGTAATATATGCCGCCGTTGCCACGGTATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTTGTAGAGCAGTTGACTTTTAATCAATTGGTCGAAGGTTCGAATCCTGCACGACTCAACAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGa                                                                                                                                                                                                       <  2:111874/251‑1 (MQ=25)                                         attTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTACGGTGGAGGCGGTAGTAAAACGTGAAGGATa                                                                                                                                                                                                     >  1:103876/1‑250 (MQ=255)                                               ggTCGTGTCGTGCGAATCATAAGCAGTTGATTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAGCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGt                                                                                                                                                                                                 <  2:234388/248‑1 (MQ=255)                                               ggTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGc                                                                                                                                                                                              >  1:483340/1‑251 (MQ=255)                                                        gtgCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAAc                                                                                                                                                                                     >  1:480580/1‑251 (MQ=255)                                                          gCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACg                                                                                                                                                                                    <  2:175891/250‑1 (MQ=255)                                                               tCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTCATATATGCCGCCGTTGCCACGGAATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGATCGCACGTTCGAATCCTGCACGACTCACCAATCGTTAAAGCGGAATCGGTAGTAAAACGTGACGGATAACGTTACATCAGCAACGGcacg                                                                                                                                                                               >  1:316013/1‑247 (MQ=37)                                                                  tAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCAAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAgg                                                                                                                                                                           >  1:506606/1‑251 (MQ=255)                                                                          ggAGTGATCTACATCGAAATTTTTGTTGCGCTAAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGAc                                                                                                                                                                    <  2:81948/249‑1 (MQ=255)                                                                                 tCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCg                                                                                                                                                             >  1:127967/1‑250 (MQ=255)                                                                                           aaTTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGc                                                                                                                                                  <  1:357946/251‑1 (MQ=255)                                                                                                     gcgcTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCGCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACccacc                                                                                                                                        <  1:504000/251‑1 (MQ=255)                                                                                                               tCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACAt                                                                                                                                >  1:364286/1‑250 (MQ=255)                                                                                                                      tCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCAAGACCCACCACTAACATagttagtt                                                                                                                        >  1:283391/1‑251 (MQ=255)                                                                                                                                      ccgTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCgtag                                                                                                        <  2:86018/251‑1 (MQ=255)                                                                                                                                        gTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCgtagt                                                                                                       >  1:441089/1‑250 (MQ=255)                                                                                                                                         gtGCCAAGGGAGAGCAAAAAAACCGAAAGCAAACAAAAAGTGGGTAGTTAGCCAAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGGCGCCGTTACGAATCCTGCACGACTCACCAAGCGCTAAGGTGGAAGCTGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGTCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCgtagt                                                                                                       <  2:204238/248‑1 (MQ=21)                                                                                                                                          tGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATc                                                                                                    >  1:175945/1‑251 (MQ=255)                                                                                                                                           gCCACGGGAAAACAAACAAACCGAAAGGAACGAAAAAGTTGGGCGTTAGCTCAGTTGGTAGAGCAGTTTGCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCTGTAAGGGGGAAGGGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACCGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATc                                                                                                    <  1:309253/250‑1 (MQ=21)                                                                                                                                             cACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCggg                                                                                                 >  1:319278/1‑251 (MQ=255)                                                                                                                                                      acaaaaaaaCCGAAAGCAACGAAAAATTGGGTCTGTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAGCCTGCACGACTCCCCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCGCCACTACCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCt                                                                                        <  2:309750/250‑1 (MQ=25)                                                                                                                                                                                  gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg                                                            >  1:160586/1‑251 (MQ=255)                                                                                                                                                                                  gggTCGGTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAGCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAagg                                                            <  2:441089/251‑1 (MQ=255)                                                                                                                                                                                    gTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGTGCGCAGGTTCGAATCCTGCACGACTCACCAATCGCGAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGTGTAGTATCGGGTTAGTAGCGCAGCTGGGAGAGCACCTCCCTTCCAaggag                                                          <  2:203812/251‑1 (MQ=25)                                                                                                                                                                                          aGCTCAGTTGGTATAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGTCTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGAGTAGCTCAGCTGGGAGAGCACCGCCCTTACAAGGAggggggc                                                    <  2:319278/251‑3 (MQ=255)                                                                                                                                                                                           gCTCAGTGGGTAGAGCAGTTGACTTTTAATCAATTGGGCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCGGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcg                                                   <  1:97149/251‑1 (MQ=255)                                                                                                                                                                                               aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcggcgg                                               >  1:428037/1‑251 (MQ=255)                                                                                                                                                                                               aGTTGGTAGAGCAGTTGACTTTTAAGCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTcggcgg                                               <  1:350683/251‑1 (MQ=255)                                                                                                                                                                                                  tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGAt                                                                                             >  1:369410/1‑202 (MQ=255)                                                                                                                                                                                                  tGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTACGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGAt                                                                                             <  2:369410/202‑1 (MQ=255)                                                                                                                                                                                                   ggTAGAGCAGTTGACTTTGAATCAATTGGTCGCAGGTTCGAATCCTTCACGACTCACCAATCGCTAAGGGGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCg                                           <  2:483340/251‑1 (MQ=255)                                                                                                                                                                                                           aGTTGACTTTTAATCAATGGGTCGCAGGTTCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGtc                                   <  2:175945/251‑1 (MQ=255)                                                                                                                                                                                                                                        tCGAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAg       <  1:448078/250‑1 (MQ=255)                                                                                                                                                                                                                                          gAATCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCCGGTCGTTAGCTCAGTTg    >  1:348739/1‑251 (MQ=255)                                                                                                                                                                                                                                           aaTCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTgg   >  1:233908/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTCCTGCACGACTCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGt  <  1:406073/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               AACGCGATGTAATGCATAACACACGACCAGAAGTCGCATTATTTCTGGTCGTGTCGTGCGAATCATAAGCAGTTGAGTGATCTACATCGAAATTTTTGTTGCGCTCAA‑GTCTGAAATCAGTAATATATGCCGCCGTTGCCACGGGATATCAAACAAACCGAAAGCAACGAAAAAGTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATCGCTAAGGTGGAAGCGGTAGTAAAACGTGAAGGATAACGTTGCATGAGCAACGGCCCGAAGGGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGT  >  NC_000913/780378‑780862

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