Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,132,077	G→A	A111V (GCG→GTG)	yghQ ←	putative transport protein YghQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,132,077	0	G	A	97.7%	145.6 / NA	44	A111V (GCG→GTG)	yghQ	putative transport protein YghQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/19);  minor base C (1/0);  total (25/19)
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

CGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCGCCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCT  >  NC_000913/3131839‑3132315                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               cGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGcgccgc                                                                                                                                                                                                                                    >  2:332209/1‑251 (MQ=255)   caacaaccgagacAAAATACCAGGCTACGACGCTCCCCGCTGCGAGCAGAACAGGTTTCGACGACTGCTGTACCGCCATTACATCAACGAGATCTACCGCACGCCGACTGCCGGTCGACGTGTAGAAAGCCATTGCAGGAATGAGCGAGCAATAGAGCGCAGCCAGCCAAAAACT‑ATGGTCATCCAGACCTAATGGATGAGAAGGGGATGACAGTAAGGCATTGCCACCGACGATCACCACcgcacagctg                                                                                                                                                                                                                                  >  2:239703/13‑244 (MQ=255)             gTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCa                                                                                                                                                                                                                        <  1:287806/251‑1 (MQ=255)                 aaaTACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCg                                                                                                                                                                                                                     >  1:164687/1‑250 (MQ=255)                               gacgCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAAACGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATgaga                                                                                                                                                                                                      <  2:127861/251‑1 (MQ=255)                                acgCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAgac                                                                                                                                                                                                     <  1:332209/251‑1 (MQ=255)                                 cgCTCACCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGCATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGAAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAgacg                                                                                                                                                                                                    >  2:144225/1‑251 (MQ=255)                                      cccGCTGCGCGCAGAAAAGGTTCCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGAACTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACAtt                                                                                                                                                                                               <  1:145443/251‑1 (MQ=255)                                       ccGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGACGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTg                                                                                                                                                                                              >  1:36922/1‑251 (MQ=255)                                           tGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGCCAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCgg                                                                                                                                                                                           >  2:340715/1‑250 (MQ=255)                                             cgcgcAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGGCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAAGGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAAt                                                                                                                                                                                        <  2:90518/251‑1 (MQ=255)                                                 cAGAAAAGGTTGCGTCGCCTGCTGTACATCAATTAAAACGAAGCGATCTACCGCACGCAGAAGGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGGCATCCAGACCTAATTAATGGGAAAGGAATGGCAGGAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCTGAATtgct                                                                                                                                                                                    <  2:315644/251‑1 (MQ=255)                                                  ggAAAAGGTTCCGGCGCCTGCTGTACAGCAAATAAACCGAAGCGATCTAACGTACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCTAAAAACT‑TTGGTCACCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGACACCGAAAATCACCGCCGCGCCGCTGATGATATCCAGCGAGAAGGAAAATGATACGACATTGctggttttgtt                                                                                                                                                                                   <  2:191501/250‑11 (MQ=255)                                                    aaaaGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCg                                                                                                                                                                                 >  1:119194/1‑251 (MQ=255)                                                            tcgtcgCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGattatt                                                                                                                                                                          <  2:254320/250‑1 (MQ=255)                                                                 gCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGACTGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGt                                                                                                                                                                    >  2:108942/1‑251 (MQ=255)                                                                 gCCTGCTGTACAGCAATTAAAGCGAAGCGATCTACCGCACGCAGAATGCCGGGCGGCGTGGAGGAAGCCATGGAAGGAATGAGCGTGCAATAGAGCGGTGCCAGCCAAGAACTATTTG‑AAGCCAGTCCTAATGATTGGGAAAGGAATGGCAGTAAGGCAGTGCCCACGACAAGCACGTTCGCGACTCTGACGAGATGCAGCGAGAAGGAAAATGATCCGAGATGGGGGAATTGGTGCGGTTTATTGTTGGt                                                                                                                                                                    <  2:417083/251‑1 (MQ=255)                                                                           cAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGCCGTGGAGCAAGCCATTGAAGCAATGAGCGTGCAATAGAGCGCTGCCACCCACAAACC‑TTGGCCATCCAGACCTAATGAACGGGAAAGGAATGGCAGCAAGGCAATGCCACCGCCAATCACCACCGCGCCGCTGACGATATCCAGAGAGAAGGAAACTGAGACGCCATTGCGGCATTGCTGCGGATTATTGTTGGTTAATGCTggg                                                                                                                                                          >  1:69463/1‑250 (MQ=255)                                                                            aGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGCGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGACAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATCTCCAGCGAGAAGGAAAATGAGACGACATTGAGGAATTGCTGCGGATTATTGTTGGTTAATGCTGgtg                                                                                                                                                         >  1:303900/1‑251 (MQ=255)                                                                              cAATTAAAGCGAAGCGATCTGCCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑GTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGtt                                                                                                                                                       <  2:336188/251‑1 (MQ=255)                                                                                                       acgcaGAAGGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTAGTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATg                                                                                                                              <  2:39020/251‑1 (MQ=255)                                                                                                                      cGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTtaataa                                                                                                               <  1:358321/251‑1 (MQ=255)                                                                                                                              aggaAGCCATTGAAGTAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGGCCTAATGAAGGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTg                                                                                                       <  1:403475/251‑1 (MQ=255)                                                                                                                                 aaGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCAACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGGCGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTTGTGTTCCGTTCTGAACCTCCAGTTGCCATTTCTTAAACTTAATAAAATCGCTGATc                                                                                                    <  2:288750/251‑1 (MQ=255)                                                                                                                                    ccATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGCATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGAc                                                                                                 >  2:271533/1‑251 (MQ=255)                                                                                                                                     cATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACt                                                                                                <  2:181761/251‑1 (MQ=255)                                                                                                                                          aaGGAATGATCGCGCAATCGAGCGCTGCCAGCCAAAAACT‑ATGGTCAACCCGACCTAATGAATGGGAAAGGCATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGAAATCCAGCGAGAAGGAAACTGAGACGACATTGCGGCATTGCTGCGGATTATTGTTGGTTAATGCAGGTGCTCAGTACTGAACCACCAGCTGCCATGTCTGAAACTTAATAACATCGCTGATCGaattagcc                                                                                           >  2:178694/1‑244 (MQ=255)                                                                                                                                            ggAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGCGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGt                                                                                          >  1:170842/1‑250 (MQ=255)                                                                                                                                                tGAGCGTGCAATAGAGCGCTGACAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCCCCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGCCTTGGCGTACGAt                                                                                     >  1:392409/1‑251 (MQ=255)                                                                                                                                                     gTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGGTGCCATGTCTGAAACGTAATAAAATCGCTGATCGACTTGGCGCACGATTGcac                                                                                <  2:74387/251‑1 (MQ=255)                                                                                                                                                       gCAATAAAGCGCTGCCAGCCCAAAACT‑TTGGTCATACAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTCTTGTTGGTTAATGCTGGTGTTCCGTACTGCAACACCAGTTGCCATGTCTGAAACTTACTAAAATCGCTGATCGACTTGGCGTACGATTGcaaa                                                                               >  2:100140/1‑248 (MQ=255)                                                                                                                                                        cAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAAGAAAATCGCTGATCGACTTGGCGTACGATTGCACAAt                                                                             <  2:104411/251‑1 (MQ=255)                                                                                                                                                                cgcTGCCAGCCAAAACCT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGCCAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGACCCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGCTCGACTTGGCGTACGATTGCACAATCACCAGt                                                                      >  2:63148/1‑250 (MQ=255)                                                                                                                                                                 gcTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTAc                                                                    <  2:119194/251‑1 (MQ=255)                                                                                                                                                                                      gTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGCTTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATggcggggggc                                                 >  1:388470/1‑246 (MQ=255)                                                                                                                                                                                                     tgaatgGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGCCGACCTTGCGGAATTGCTGCGGATTCTTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATACAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGGCATCCCTTTACCGGcc                                  >  1:299832/1‑249 (MQ=255)                                                                                                                                                                                                     tgaatgGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGcac                                 >  1:341716/1‑251 (MQ=255)                                                                                                                                                                                                          gggAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGGCATCCCTTTACCGGCACACg                              >  1:410617/1‑249 (MQ=255)                                                                                                                                                                                                                ggAATGGCAGTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGcc                      >  1:128222/1‑251 (MQ=255)                                                                                                                                                                                                                   aTGGCAGTAAGGCAATGACACCGACAATCACCACCGCGCCGCTGACGATATACAGCGAGAAGGAAGATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTACTAAACTCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGATCGCCacc                   >  2:150033/1‑249 (MQ=255)                                                                                                                                                                                                                         gTAAGGCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGAccc             >  1:144788/1‑251 (MQ=255)                                                                                                                                                                                                                              gCAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGGACGCCAAACATGGCGGGCGGTATCCCCTTATCGGCACACGAGAGGGCCACCAGaacccccgc        >  2:118901/1‑243 (MQ=255)                                                                                                                                                                                                                               cAATGCCACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGGGTCATACCTTTAACGGCCCACCAGAGCGCCAACAgcaccctaa         >  2:160127/1‑241 (MQ=255)                                                                                                                                                                                                                                    ccACCGACAATCACCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCt  <  1:288945/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               CGCAAAACCAAAGTCAAAATACCAGGCTACGACGCTCCCCGCTGCGCGCAGAAAAGGTTTCGTCGCCTGCTGTACAGCAATTAAATCGAAGCGATCTACCGCACGCAGAATGCCGGTCGGCGTGGAGGAAGCCATTGAAGGAATGAGCGTGCAATAGAGCGCTGCCAGCCAAAAACT‑TTGGTCATCCAGACCTAATGAATGGGAAAGGAATGGCAGTAAGGCAATGCCACCGACAATCGCCACCGCGCCGCTGACGATATCCAGCGAGAAGGAAAATGAGACGACATTGCGGAATTGCTGCGGATTATTGTTGGTTAATGCTGGTGTTCCGTACTGAACCACCAGTTGCCATGTCTGAAACTTAATAAAATCGCTGATCGACTTGGCGTACGATTGCACAATCACCAGTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCT  >  NC_000913/3131839‑3132315

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