Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,003,502	G→T	P226Q (CCG→CAG)	rarD ←	putative transporter RarD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,003,502	0	G	T	97.8%	125.3 / NA	47	P226Q (CCG→CAG)	rarD	putative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (17/28);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAG  >  NC_000913/4003258‑4003731                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           cGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAACAACAGctgtacgg                                                                                                                                                                                                                                      >  1:738074/1‑251 (MQ=255)    cATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAg                                                                                                                                                                                                                                   <  2:188325/251‑1 (MQ=255)             tCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGccgtacgggagggaaacaac                                                                                                                                                                                                                          >  2:739402/1‑239 (MQ=255)              cAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATAc                                                                                                                                                                                                                          <  1:8409/250‑1 (MQ=255)                      atttatttTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGgcggcg                                                                                                                                                                                                                  <  1:777961/250‑1 (MQ=255)                         tattttttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGTGCCAATGTACTGGAAAAAGCCTAGCGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCTATc                                                                                                                                                                                                               <  1:475495/250‑1 (MQ=255)                         tattttttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATc                                                                                                                                                                                                               <  1:334501/250‑1 (MQ=255)                         tattttttatttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATc                                                                                                                                                                                                               <  1:475482/250‑1 (MQ=255)                                   ttGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTAAGGTAGGGACAAAAACGGCGGGGGTCAGCAGTAAAtt                                                                                                                                                                                                    >  1:615066/1‑251 (MQ=255)                                         gTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATca                                                                                                                                                                                                              <  1:59131/235‑1 (MQ=255)                                         gTTCTACGCTGCGTATAAATCGCAACCATCACACAAATTGCCAGCGCAACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACCACTGTACGGGAGTGACAATAACGGCGGCGATCAGCAGTAAATTCaccgg                                                                                                                                                                                              >  1:286174/1‑247 (MQ=255)                                              aCGCTGAGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAACCGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCATAACGTTGAGAGACGCAAGCGCGTGGCAGCGGAGGTAAAAAACAAAAACCAGACGGTAGAGAAAATACGGGCGGACATTAGCAGTAAAttatgcgcactcg                                                                                                                                                                                         >  1:93300/1‑240 (MQ=255)                                                   gCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGACACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGt                                                                                                                                                                                     <  1:179593/249‑1 (MQ=255)                                                      tataAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACACAAGCTGTACGGTAGTGGCAATACCGGCGGCGATAAGCAGTAAATTAAGCGACATC‑GGGtttttc                                                                                                                                                                                >  1:105874/1‑249 (MQ=255)                                                             cGCATCCAGCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTCCACCATAAAACGTCGCAGCCAGCAGGAACATCAGCGACGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTGGCCCCatat                                                                                                                                                                         <  1:429514/251‑1 (MQ=255)                                                                atccatcACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAACAACAGCTGTACGGGAGGGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACAAC‑GGGTTTTGCCCCATATgtc                                                                                                                                                                      >  1:533482/1‑249 (MQ=255)                                                                                gCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTTCCCCATATTGCTGGT‑TGAGCTGTCggg                                                                                                                                                     >  1:671828/1‑250 (MQ=255)                                                                                  cAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAAAAGCTGTACGGGAGTGAAAATACCGGCGGCGATAAACAGTAAAATCAGCGACACC‑GGGTTTTTCCCCCTATGGCTGGT‑TGAGCTTTCGGCa                                                                                                                                                    >  1:493858/1‑250 (MQ=255)                                                                                   aGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGGTTTGCCCCATATGGCAGGG‑TGAGCTGTCGGCaa                                                                                                                                                   >  1:816135/1‑250 (MQ=255)                                                                                                            tCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGGAAAACACAAAAGCTGTAAGGGAGTGACAATAACGGCGGCGGTCAGCCGGAAAATCAGCGGCATC‑GGGGTTTTCCCCCTATGGCTGGG‑TTAGCCTTCCGCAATA‑ACAAACAG‑GGAAAATtccgccccc                                                                                                                       >  2:101772/1‑248 (MQ=255)                                                                                                               ccATCTTATCGTCACCCTGTTTTGCACCATAACACGTCACAGCCAGCAGGACCACCAGCGTCGGGCCAATTTACTGTAAAATGCCTAACGTCGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTATTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGTCATC‑GGGTTTTGCCCCATATGGCTGGT‑TGATCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACggg                                                                                                                     <  1:787981/250‑1 (MQ=255)                                                                                                                        cGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACCACAGCTGTACGGGAGGGACAAAAACGGCGGCGGTCAGCAGTAAAATCAGCGACATC‑GGGTTTTGCCCCCTAAGGCTGGT‑TTAGCTTACGGCAATA‑GCAAAAAG‑GTAAATTTCCGGCACCGGcagcaggca                                                                                                            >  1:691790/1‑247 (MQ=255)                                                                                                                             cccGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGGAGTGAAAATAACGGCGGCGATCAGCAGTAAATTTAGCGACATC‑GGGGTTTGCCCCATATTGCTGGT‑TGAGATGGCGGCAATA‑ACAAAAAA‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGt                                                                                                     >  1:789707/1‑251 (MQ=255)                                                                                                                                ggTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTc                                                                                                  <  1:101772/251‑1 (MQ=255)                                                                                                                                ggTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAAAAAAACAGCTGGAAGGTAGTGACAATACCGGCGGGGGTCAGCCGGAAAATAAGCGCCATC‑GGGTTTTGCCCCCTATGGGTGGGGTGAGCTGT‑GGCACTA‑AAAAACAGAGCAAATTgccgcccggggcggaggcaaccaggtgt                                                                                                    >  2:751104/1‑229 (MQ=255)                                                                                                                                      tCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTa                                                                                             <  1:387136/250‑1 (MQ=255)                                                                                                                                       cACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTACCGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAAc                                                                                           <  1:82633/251‑1 (MQ=255)                                                                                                                                                    cacaGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGGTTTGCCCCATATGGCTGGT‑TTAGCTGTCGGCAATA‑ACAAAAAG‑GGAAATTGCCGCCACGGGCAGCCGCCAC‑ATGGGTTCGATTAACATGCCGGttttgg                                                                              >  1:428354/1‑248 (MQ=255)                                                                                                                                                      cagccagcAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTGGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGTCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCt                                                                            <  2:738074/251‑1 (MQ=255)                                                                                                                                                       agccagcAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCTCTTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGTTCAGCAGTAAATTCATCTACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCtt                                                                           <  1:577201/251‑1 (MQ=255)                                                                                                                                                       agccagcAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAACTGTACGGGAGTGACAAAACCGGCGGCGATCAGCAGTAAATTCAGCGACACCCGGGTTTT‑CCCCATATTGCTGGT‑TGAGCTGGCGGCAAAA‑ACAAACAG‑GGAAATTGCCCCCCCGGGCAGCACCCACCATG‑TTTCGATTAACATACCGGTTTGCGCtt                                                                           >  1:119020/1‑251 (MQ=255)                                                                                                                                                            gcagGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACg                                                                      <  1:547978/251‑1 (MQ=255)                                                                                                                                                                    atgagtggggggcaagtttatcggAAAAAGTCTAACGTTGAGAGACGCTAGCTCGTGTGAGCGGCGGTAAAACACAACAGCTGGACGGTAGTGACAATACCGGCGGCGATCATCAGTAAATGCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GGAATCAG‑GTAAATTGCCGCCAGGGGCAGCAGCCAC‑ATGGTTGCGAGTAACATGCCGGTTTGCGCTGCAACGGCAAtctt                                                              <  2:743910/229‑1 (MQ=255)                                                                                                                                                                          gTCGGGCCAATGGACAGGAAAAAGCCGAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAAGTCAGCGACATC‑GGGTTGTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCg                                                        <  2:691790/251‑1 (MQ=255)                                                                                                                                                                               gCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCa                                                   <  1:585296/251‑1 (MQ=255)                                                                                                                                                                                 cAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGctgta                                                                                                                                                                                                                                         <  2:709429/72‑1 (MQ=255)                                                                                                                                                                                 cAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGAAAAACACAACAGctgta                                                                                                                                                                                                                                         >  1:709429/1‑72 (MQ=255)                                                                                                                                                                                       ccTGGAAAAAGCCTAACGGTGAGAGACGCAAGCGCGGGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAg                                           <  1:443750/250‑1 (MQ=255)                                                                                                                                                                                              aaaGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCaa                                    <  1:269491/251‑1 (MQ=255)                                                                                                                                                                                                aGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGTTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCaaaa                                  <  2:615066/251‑1 (MQ=255)                                                                                                                                                                                                     aacgtggagggccGCAAGGGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTaa                              <  2:789707/238‑1 (MQ=255)                                                                                                                                                                                                     aaCGGTGAGAGACCAAAGCGCGTGGCAGCGGCGGTAAAACTCAACAGCTGTACGTTAGTGAAAATACCGGCGGCGATCAGCAGTAAATCCAGCTACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATGTCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGTCAATCTTCTTGCGTACCAGACCTTAGAAGGCAAAACTaaa                             <  2:760346/251‑1 (MQ=255)                                                                                                                                                                                                        gTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTCCGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCAACGGGCAGCAGCCAC‑ATGGTTTCTATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCGTGCGTACCAGACCGTAGAAGGCAAAACTAAATGc                          <  2:611257/251‑1 (MQ=255)                                                                                                                                                                                                          tGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGcc                         <  1:279655/250‑1 (MQ=255)                                                                                                                                                                                                                   cAAGCGCGTGGCAGCGGCGGTAAAACACAACAGGTGTACGGAAGTGACAATACCCGCGGCGATCAACAGTAAATTCAGCGACATA‑GGGATTTGCCCCATAAGGCCGGG‑TGAGCTGTAGGCAATAACCAAAAAG‑G‑AAATTGCCGCCACGGGCAACAACCAA‑AAGGGTTCGATTAACAAACCGGTTTGCGCTTCAAAGGCCATCTTCTTGAGTAACAAGCCCTGGAAGGCAAAAAAAAATGGCAGTCCCaacg               >  1:506402/1‑248 (MQ=255)                                                                                                                                                                                                                    aaGCGCGTGGCAGCGGCGGTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGGTCAGCAGTAAATTCAGCGACATC‑GGGGTTTGCCCCATATGGCTGGT‑TGAACTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCAAGGGCAGCAGCCAC‑ATGGGTTCGATTAACATGCCGGGTTGCGCTTCAACGGCAAACTTCTTGCGTACCCGACCGTAGAAGGCAAAAATAAATTCCAGTCCCAgcgc              >  1:513502/1‑251 (MQ=255)                                                                                                                                                                                                                              cAGCGGCGGTAAAACCCAAAAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGGCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGt    <  2:428354/251‑1 (MQ=255)                                                                                                                                                                                                                              cAGCGGCGGTAAAACACAACAGCTGTACGGGAGTGACAATACCGGCGGCGATCAGCAGGAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAACTGTCGGCAATA‑GCAAACAG‑GGAAATTGCCGCCCCCGGCAGCAGCCAC‑ATGGGTTTGATTAAAATGCCGGGTTGCGCTTCAACGGCAATCTTCTTGCGGACCAGACCGGAGAAGGCAAAACTCAATGCCAGGCCCACCGCGaaaaaaggg    >  2:287508/1‑243 (MQ=255)                                                                                                                                                                                                                                gcggcgGTAAAACACAACAGCTGGACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGAAATC‑GGGTTTTGCCCCATAGGGCTGTT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAg  <  2:306150/251‑1 (MQ=255)                                                                                                                                                                                                                                     ggTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAAt                                                                                                                                 >  2:222368/1‑120 (MQ=255)                                                                                                                                                                                                                                     ggTAAAACACAACAGCTGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAAt                                                                                                                                 <  1:222368/120‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                           CGGCATATTTTTTCAGAATATATTTATTTTTTATTTGGACGTTCTACGCTGCGTATAAATCGCATCCATCACAAAAATTGCCAGCGCCACCCAAATAAAGGCGAAAGTCACCATCTTATCGGCACCCGGTTTTTCACCATAAAACGTCACAGCCAGCAGGAACATCAGCGTCGGGCCAATGTACTGGAAAAAGCCTAACGTTGAGAGACGCAAGCGCGTGGCAGCGGCGGTAAAACACAACAGCGGTACGGTAGTGACAATACCGGCGGCGATCAGCAGTAAATTCAGCGACATC‑GGGTTTTGCCCCATATGGCTGGT‑TGAGCTGTCGGCAATA‑GCAAACAG‑GTAAATTGCCGCCACGGGCAGCAGCCAC‑ATGGTTTCGATTAACATGCCGGTTTGCGCTTCAACGGCAATCTTCTTGCGTACCAGACCGTAGAAGGCAAAACTAAATGCCAGTCCCAGCGCGATAATAGGTAG  >  NC_000913/4003258‑4003731

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