Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	931,091	T→C	A320A (GCA→GCG)	trxB ←	thioredoxin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	931,091	0	T	C	97.3%	119.7 / ‑4.0	37	A320A (GCA→GCG)	trxB	thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (21/15);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

ATGATCAATATGCCGCTCACAAAGCCCAGCAGACGAGAAATATTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTT‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGAC  >  NC_000913/930847‑931335                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      aTGATCAATATGCCGCTCACAAAGCCCAGCAGACGAGAAATATTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCa                                                                                                                                                                                                                                                <  2:86800/251‑1 (MQ=255)   gATCAATATGCCGCTCACAAAGCCCAGCAGACGAGAAATATTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGc                                                                                                                                                                                                                                              >  1:192140/1‑251 (MQ=255)     tCAATATGCCGCTCACAAAGCCCAGCAGACGAGAAATATTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCt                                                                                                                                                                                                                                             >  1:131220/1‑250 (MQ=255)                  cacaAAGCCCAGCCGGCGAGAAATCTTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGt                                                                                                                                                                                                                                >  2:224589/1‑250 (MQ=255)                        gCCCAGCAGACGAGAAATATTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTGTTGTTACTGATTTGTAAAATTAGTAC‑GCGTCAGCTAATCCAGCGAGGTAGCGTTc                                                                                                                                                                                                                         <  2:371631/251‑1 (MQ=255)                        gCCCAGCAGACGAGAAATATTCAGCCAACCTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTc                                                                                                                                                                                                                         <  1:293813/251‑1 (MQ=255)                                        atatTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGcc                                                                                                                                                                                                         >  1:15320/1‑251 (MQ=255)                                                ccAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGcc                                                                                                                                                                                                 >  2:128140/1‑251 (MQ=255)                                                ccAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGcc                                                                                                                                                                                                 >  1:1616/1‑251 (MQ=255)                                                     gTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTAc                                                                                                                                                                                            >  1:139706/1‑251 (MQ=255)                                                                          gTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGc                                                                                                                                                                       <  1:247128/251‑1 (MQ=255)                                                                                               ctTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATTCATCAc                                                                                                                                                    <  2:325880/249‑1 (MQ=255)                                                                                                     gACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTg                                                                                                                                            <  2:331226/251‑1 (MQ=255)                                                                                                       cGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCg                                                                                                                                          <  2:160348/251‑1 (MQ=255)                                                                                                         agaTTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGc                                                                                                                                        >  2:326865/1‑251 (MQ=255)                                                                                                         agaTTGATTCATTGCACTCGTAGCAGGTGAGTTATCATAATTATTTGCATAGCAATGGTACAACTGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTTCTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTATCGTTCCGCATCAAGTGCTGCCATTCAGCCTGTCCCTTCCGAAGTAATTGACTTGCGATAAATGTGATCCAGCACGTCTCCTGCGGc                                                                                                                                        <  2:348240/251‑1 (MQ=255)                                                                                                              tATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCaa                                                                                                                                      >  1:228728/1‑248 (MQ=255)                                                                                                              tATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCaa                                                                                                                                      >  1:82190/1‑248 (MQ=255)                                                                                                                               gCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTc                                                                                                                  >  2:106269/1‑251 (MQ=255)                                                                                                                               gCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTc                                                                                                                  >  2:276857/1‑251 (MQ=255)                                                                                                                                 aGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGGCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAAGGTGATCCAGCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTg                                                                                                                <  1:187430/251‑1 (MQ=255)                                                                                                                                                   attTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGAATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAAt                                                                                              >  2:345913/1‑251 (MQ=255)                                                                                                                                                       gCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATAccc                                                                                          <  2:81168/251‑1 (MQ=255)                                                                                                                                                               aTGTTCCAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCGGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCGGGTCTGGGTGGCATTACCATGAATACCCGACTGTAc                                                                                  <  2:115388/251‑1 (MQ=255)                                                                                                                                                                   tACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTg                                                                              >  1:275597/1‑251 (MQ=255)                                                                                                                                                                                       ggCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTttccagttc                                                          >  1:22498/1‑251 (MQ=255)                                                                                                                                                                                               aTAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGcc                                                  >  1:72074/1‑251 (MQ=255)                                                                                                                                                                                                    cGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTc                                              <  2:341233/250‑1 (MQ=255)                                                                                                                                                                                                      cATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCCTTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGa                                            >  2:154790/1‑250 (MQ=255)                                                                                                                                                                                                       aTGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGaaa                                          <  1:128140/251‑1 (MQ=255)                                                                                                                                                                                                              cGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAg                                   <  1:276857/251‑1 (MQ=255)                                                                                                                                                                                                               gcattcttTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCAGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGATGCCCTTCGAAAATCGc                                     >  2:336450/4‑248 (MQ=255)                                                                                                                                                                                                                   tcttTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCTTCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTAtt                               <  1:139662/250‑1 (MQ=255)                                                                                                                                                                                                                   tcttTACTTTTGTTACTGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTATTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAACTCGCAGTATTc                              >  2:14839/1‑251 (MQ=255)                                                                                                                                                                                                                                   tGATTTGTAAAATTATTTC‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGAt                >  1:244161/1‑249 (MQ=255)                                                                                                                                                                                                                                            aaaatattTTCGCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCaacaaac      <  2:197593/247‑1 (MQ=255)                                                                                                                                                                                                                                               ttattTC‑GCTTCAGCTAAACCATCGAGGTAGCGTTCCGAATCAAATGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGCTGTAGCCGTTTTCCAGTTCCAGCTGCACTTCGAAAATCGCAGTATTCGGGCTGTGACAGATAGCAACAAACAGAc  >  1:239563/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      ATGATCAATATGCCGCTCACAAAGCCCAGCAGACGAGAAATATTCAGCCAACGTTGGGAGATGACGCTTTGCTGTTTTAACCAGCGGGTTAACTCTTTTTGACGAGATTTATTCATTGCACGCTTAGCAGGTGAGTTATCAGAATTATTTGCAGAGCAATGTTACAACGGGGAAAAAATAAAGGCGACCCATAGTCGCATGGTGTCGCCTTCTTTACTTTTGTTACTGATTTGTAAAATTATTTT‑GCGTCAGCTAAACCATCGAGGTAGCGTTCCGCATCAAGTGCTGCCATGCAGCCTGTACCGGCCGAAGTAATGGCCTGGCGATAAATGTGATCCATCACGTCGCCTGCGGCAAAGACGCCAGGAATGCTGGTCTGGGTGGCATTACCATGAATACCCGACTGTACTTTGATGTAGCCGTTTTCCAGTTCCAGCTGCCCTTCGAAAATCGCAGTATTCGGGCTGTGACCGATAGCAACAAACAGAC  >  NC_000913/930847‑931335

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