Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,931,779	T→C	F155S (TTT→TCT)	kup →	K(+):H(+) symporter Kup

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,931,779	0	T	C	96.8%	205.3 / ‑4.3	63	F155S (TTT→TCT)	kup	K(+):H(+) symporter Kup
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (22/39);  total (23/40)
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.83e-01

GTGATGCGGGCAGATAACGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATTTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCT  >  NC_000913/3931529‑3932014                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             gTGATGCGGGCAGATAACGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGCTACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTattatt                                                                                                                                                                                                                                              >  1:15917/1‑251 (MQ=255)       cGGGCAGATAACGCCGGTGAGGGGGGGATCCTTACGTTGATGTCGCTTGCCGTGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGAt                                                                                                                                                                                                                                        <  1:396637/251‑1 (MQ=255)            aGATAACGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑Caa                                                                                                                                                                                                                                   >  1:101671/1‑250 (MQ=255)              aTAACGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑Caaaa                                                                                                                                                                                                                                 >  1:639558/1‑250 (MQ=255)                   gCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGCTACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATgg                                                                                                                                                                                                                           >  2:118845/1‑251 (MQ=255)                      gtgaggggggggTTCCGGACGTTGATGTCGCTTGCCTGGGGTAATCCGTCGGCGCGAACCAACTCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCGCACTTTATTATCTATGATT‑CAAAAACATGGCAc                                                                                                                                                                                                                        <  2:15917/246‑1 (MQ=255)                       ggggaggggggATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACg                                                                                                                                                                                                                       <  1:668602/246‑1 (MQ=255)                         gAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACgg                                                                                                                                                                                                                      <  1:447692/250‑1 (MQ=255)                             gggggATCCTGACGTTGATGTCGCTTGCCTGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGGTTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGGGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGCTAGTTCCCCTCTCAATTATCGTTCTCACCTTCTTATCTATGATT‑CCAAAGCATGGCACAGCTATg                                                                                                                                                                                                                 >  2:159177/1‑251 (MQ=255)                                         cGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTg                                                                                                                                                                                                     <  1:66295/251‑1 (MQ=255)                                            tGATGTCGCTTGCCGGGCGTAATACATCGGCGCGAACCACATAAATGCTGGTGCTTCTGGGGCTAATCGGCGGACGCTTTTTCTATGGTGCAGTCGACATAAAACCCGCTATTTCGGAGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCCGCTGGATACCTGGATAGTTCACCTCTCAATTATCGTTCTACCTTTATTATCTTTGTTT‑CAGAAACATGCAACGGCTATGGTCggacacatgttt                                                                                                                                                                                                  >  2:643547/1‑241 (MQ=255)                                               ggtcgtgttgcGGGCGTAATACGTAGGCGCGACCCACATCAATGCTGGTTATTATGGGGCTAATCGGGGGAAGCGTTTTCTAGGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAATGCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCGCGCAATGATCGTTCTCACTTTATTATCTATGATG‑CAAAAACATGGCACGGCTATGGTCGGTAATCTGTTTgcg                                                                                                                                                                                               <  2:630293/241‑1 (MQ=255)                                                  cGCTTGCCGTGCTTAATACGTCGGCGCGAACCCTATCAATGCGGGTGATTATTGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGAAGACATACCACCCGCTCTTTCGGTGATGTCAGCCCTTGAAGGTCTGGAAAGCGTTGCCCCGCAGCTGGATTCCTTGATAGTTCCCCTCTCAATGCTCGTCGTCACTTTATTATCTATGATT‑CAAAATCATGGCTCGGCTATGGTCGGTAAGCGGTTGGCGcc                                                                                                                                                                                             <  2:334667/250‑1 (MQ=255)                                                      tGCCGGGCGTAATCCGTCTGCGCGCACTACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGAt                                                                                                                                                                                          <  1:525797/249‑1 (MQ=255)                                                        ccGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATg                                                                                                                                                                                      >  1:233018/1‑251 (MQ=255)                                                               tAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCt                                                                                                                                                                                                                    >  1:36501/1‑214 (MQ=255)                                                               tAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCt                                                                                                                                                                                                                    <  2:36501/214‑1 (MQ=255)                                                                     gTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGtttt                                                                                                                                                                         <  1:483429/251‑1 (MQ=255)                                                                                  acaTCAATGCTGGTGATTATGAGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CACAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAgg                                                                                                                                                            >  1:574847/1‑251 (MQ=255)                                                                                          gCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGAGAAGACGTCATAACACCCGCTATTTCGGAGATGTCAGCCATTGAAGGACTGGAAATAGTCGCCCAGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCAAGGCTATGGTGGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGGTGTTGATTCTGGCAGGACTGGGGt                                                                                                                                                    >  1:267916/1‑251 (MQ=255)                                                                                           catgtgtttATGGGGTTAATCGGCGGCAGCTTTTTCTATGGTGATGTCGTCATAACACCCTGTATTGCGGAGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATTCCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGTTGTTTGCGCCGATCATGCTGACCTGTTTTTGGATTCTGGCAGGTCTTGGGtt                                                                                                                                                   <  2:480740/244‑1 (MQ=255)                                                                                             ggTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTa                                                                                                                                                  <  1:439722/250‑1 (MQ=255)                                                                                                gATTATGGGGCTAATCGGCGGCAGCTTGTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGGAGCTGGATGCCTGGATAGTTCCCCCCTCAATTGTCGTTCTCACTTTATTATTTATGATG‑CAAAAACATGGCTCGGCTGTGGTCGGTAAGCTGGTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGt                                                                                                                                               <  1:285198/250‑1 (MQ=255)                                                                                                   tATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGc                                                                                                                                            <  2:639558/250‑1 (MQ=255)                                                                                                         gCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCGCACTTTAGTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTg                                                                                                                                     <  1:643547/251‑1 (MQ=255)                                                                                                          cTAATCGGCGGCAGCTTTTTCTGAGGTGAAGTCGTCCTAACACCCTCTATTTCGGTGATGTCAGCCCTTGAAGGTCTGGAATTCGACGCCCCCCAGCTGGATACCAGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAACCATGGCACGGCTATGGTATGTAAGCTGTTTGCTCCGATCATGCTGACCTGGTTTTTGGTTCTGGCAGGACTGGGTTTCCGTacctttcttgc                                                                                                                                    >  2:540050/1‑241 (MQ=255)                                                                                                                      aGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGt                                                                                                                        <  1:342301/251‑1 (MQ=255)                                                                                                                              cTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGttttt                                                                                                                                                                        >  1:28340/1‑195 (MQ=255)                                                                                                                              cTATGGTGAAGTCGTCATAACACCCGCTATTCCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCAGGCGGACCTGGttttt                                                                                                                                                                        <  2:28340/195‑1 (MQ=255)                                                                                                                                aTGGTTAAGCCGTCATACCCCCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATg                                                                                                                <  2:151483/249‑1 (MQ=255)                                                                                                                                aTGGTGAAGTCGTCATAACACCCGCTATTTTGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTAAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATgc                                                                                                               <  1:539906/250‑1 (MQ=255)                                                                                                                                   gggAAGTCGGCATAATGCCCGCTATTTCGGGGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCTCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACGTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGTTTTCGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTg                                                                                                           <  2:268827/249‑1 (MQ=255)                                                                                                                                          cgtcATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCa                                                                                                     <  1:270807/250‑1 (MQ=255)                                                                                                                                                aaCACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTggg                                                                                              <  1:557346/251‑1 (MQ=255)                                                                                                                                                          ttttCGGTGATGCCAGCGATTGAAGGTCTGGAAATCTTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCGCGCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGGAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCAtt                                                                                     <  2:285591/249‑1 (MQ=255)                                                                                                                                                          aTTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCTCTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCAtt                                                                                     <  1:256748/250‑1 (MQ=255)                                                                                                                                                              cGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATttctt                                                                                 <  1:354528/250‑1 (MQ=255)                                                                                                                                                                  gATGGCAGCCATTGAAGGGCTGGAAATCGTCGCCCCGCAGCTTGATACCCGGATAGTTCCCCTCTCAATTATCGTTCTCACTTGATTATCTATGATT‑CAAAAACATGGCACGGCGATGGTCGTTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATGCTGGCATGACTGGGGTTACGTAGCATGATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGTGCTGTGCATTTCTTCCTTg                                                                            <  2:243874/251‑1 (MQ=255)                                                                                                                                                                   aTGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACt                                                                                                                                                         >  2:104707/1‑173 (MQ=255)                                                                                                                                                                   aTGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACt                                                                                                                                                         <  1:104707/173‑1 (MQ=255)                                                                                                                                                                     gTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAAt                                                                         <  1:458384/251‑1 (MQ=255)                                                                                                                                                                        aGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTTGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGAGT‑CAAAAACATGGCACGGCTATGGTCGGTAATCTGGTGGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACa                                                                      <  2:120304/251‑1 (MQ=255)                                                                                                                                                                        aGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTTGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACa                                                                      <  1:100163/251‑1 (MQ=255)                                                                                                                                                                        aGCCATTGAAGGACTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGAGGGCGGTGACTTTCTTCCTTGAAttc                                                                       >  2:17931/1‑248 (MQ=255)                                                                                                                                                                            aTTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAAc                                                                  >  1:571111/1‑251 (MQ=255)                                                                                                                                                                             ttGAAGGTCTGGAAATCGTCGCCCCGCACCTCGATCCCTGGATAGTTCCCCTCTCAATTATCGATCTCACTTAATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTACGCTGTTTGCGCCGATCAAGCTGACCTGGTTTTTGATTCTGGCAGGCCTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGTTGAGTGCGATGACTCCAATGAGGTCGGGGAATTTCTTCTTTGAATACaaac                                                                   >  2:330417/1‑248 (MQ=255)                                                                                                                                                                                   gTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCAt                                                                                      >  2:123845/1‑224 (MQ=255)                                                                                                                                                                                   gTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCAt                                                                                      <  1:123845/224‑1 (MQ=255)                                                                                                                                                                                   gTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCt                                                           <  1:644637/251‑1 (MQ=255)                                                                                                                                                                                          aaTCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTg                                                    <  1:422165/251‑1 (MQ=255)                                                                                                                                                                                                   cccGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAgggg                                           <  1:207092/251‑1 (MQ=255)                                                                                                                                                                                                    ccGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGATCTCACTTTATTTTTAATGCTTACAAAAA‑TCGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTAATTCTGGCAGGACTGGGGTTACGTAGAATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGAGGTTGCAGCATTTCTTCCTTGAATACAAAACGTTTTCTTTTATTGCATTAGGGGc                                          >  1:395983/1‑251 (MQ=255)                                                                                                                                                                                                       cAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATgc                                                                                                               >  2:672375/1‑179 (MQ=255)                                                                                                                                                                                                       cAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATgc                                                                                                               <  1:672375/179‑1 (MQ=255)                                                                                                                                                                                                        aGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAgtg                                      <  1:195941/251‑1 (MQ=255)                                                                                                                                                                                                         gCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGCCACGGCTCTGGTCGGTAAGCTTTTTGCGCCGATCATGCTGACCTGGTTTTTGTTTCTGGCAGGACTGGGGTTACGTAGCATTATTTCTAACCCGGAAGTGCTGCATGCGCTGTATCCAATGGTGGCGTTACATTTCTTCCTTGAATTCAAAACGGTTTCTTTTATTGCATTCGGGGCAgtgg                                     >  2:549696/1‑251 (MQ=255)                                                                                                                                                                                                              attccTGGTTAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAAGCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGtgttg                                <  2:671705/248‑1 (MQ=255)                                                                                                                                                                                                              aTACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTGCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGtgttg                                <  1:281348/251‑1 (MQ=255)                                                                                                                                                                                                                tccTGGATAGTTCCCCTCTCAATTAACGTTCTCATTGTAGTATCTATGATT‑CAAAAACATGGCATGGCTATGGTCGGTAAGCTGTTTGCGCCGATCAAGCTGACCTTGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCCGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTGGAATACAAAACGTTTTCTTTTATTGCATTAGGGGCAGTGGtgttgt                               <  1:274565/249‑1 (MQ=255)                                                                                                                                                                                                                    ggATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGAtt                          <  2:208104/251‑1 (MQ=255)                                                                                                                                                                                                                          ttCCCCTCTCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACgggg                    >  1:430685/1‑251 (MQ=255)                                                                                                                                                                                                                               ctctCAATTATCGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCg                >  1:164555/1‑250 (MQ=255)                                                                                                                                                                                                                                          cGTTCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTCTTGCATTAGGGGCAGTGGTGTTGTCGATAACGGTGGTCTAGGCCCTGTATg    >  2:12516/1‑251 (MQ=255)                                                                                                                                                                                                                                            ttCTCACTTTATTATCTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCt  >  2:502454/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             GTGATGCGGGCAGATAACGCCGGTGAGGGGGGGATCCTGACGTTGATGTCGCTTGCCGGGCGTAATACGTCGGCGCGAACCACATCAATGCTGGTGATTATGGGGCTAATCGGCGGCAGCTTTTTCTATGGTGAAGTCGTCATAACACCCGCTATTTCGGTGATGTCAGCCATTGAAGGTCTGGAAATCGTCGCCCCGCAGCTGGATACCTGGATAGTTCCCCTCTCAATTATCGTTCTCACTTTATTATTTATGATT‑CAAAAACATGGCACGGCTATGGTCGGTAAGCTGTTTGCGCCGATCATGCTGACCTGGTTTTTGATTCTGGCAGGACTGGGGTTACGTAGCATTATTGCTAACCCGGAAGTGCTGCATGCGCTGAATCCAATGTGGGCGGTGCATTTCTTCCTTGAATACAAAACGGTTTCTTTTATTGCATTAGGGGCAGTGGTGTTGTCGATTACGGGGGTCGAGGCGCTGTATGCT  >  NC_000913/3931529‑3932014

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