Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	495,680	C→T	F187F (TTC→TTT)	htpG →	chaperone protein HtpG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	495,680	0	C	T	97.9%	166.6 / NA	48	F187F (TTC→TTT)	htpG	chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (26/21);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAACCG  >  NC_000913/495437‑495919                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                cGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGa                                                                                                                                                                                                                                            <  2:117948/249‑1 (MQ=255)                 gACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTc                                                                                                                                                                                                                          <  2:348516/251‑1 (MQ=255)                                 gccagcTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATAc                                                                                                                                                                                                          >  1:400199/1‑251 (MQ=255)                                         aTCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGCAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTTCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGAcc                                                                                                                                                                                                   >  1:386054/1‑250 (MQ=255)                                           cGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGCCATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGCG‑TTTCTCGATGACTGGCGCGTGCCTTCCATCATCAGCAAATACTCCGACCat                                                                                                                                                                                                 >  2:75297/1‑250 (MQ=255)                                           cGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCat                                                                                                                                                                                                 >  1:160417/1‑250 (MQ=255)                                            ggTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCatat                                                                                                                                                                                               >  1:232510/1‑251 (MQ=255)                                                      gtgtTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCg                                                                                                                                                                                     <  2:472792/251‑1 (MQ=255)                                                        gtTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGt                                                                                                                                                                                   <  1:266958/251‑1 (MQ=255)                                                         ttGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGt                                                                                                                                                                                   >  1:358597/1‑250 (MQ=255)                                                          tGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTa                                                                                                                                                                                  <  1:172700/250‑1 (MQ=255)                                                            gTTTCTCCTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTagag                                                                                                                                                                               <  2:155225/251‑1 (MQ=255)                                                               tCTACTCTGCGTTTATCGTGTCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCAGCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACGGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCTGTAGAGAt                                                                                                                                                                             <  2:353706/250‑1 (MQ=255)                                                               tCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATTGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGAt                                                                                                                                                                             <  2:128238/250‑1 (MQ=255)                                                                  aCTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATc                                                                                                                                                                            <  1:131681/248‑1 (MQ=255)                                                                  aCTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGa                                                                                                                                                                          >  1:338977/1‑250 (MQ=255)                                                                                gTGGCCGACAACGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAgaga                                                                                                                                                            >  1:134449/1‑250 (MQ=255)                                                                                   gCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGTAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGAcc                                                                                                                                                                                                   <  1:104946/208‑1 (MQ=255)                                                                                     cGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGCGAATACACAGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACg                                                                                                                                                      >  1:475499/1‑251 (MQ=255)                                                                                                                   aGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATc                                                                                                                         >  2:533187/1‑250 (MQ=255)                                                                                                                      cGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCCTCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATcaac                                                                                                                      >  1:57131/1‑250 (MQ=255)                                                                                                                       gAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGGGAATACACCGTTGCCGGCATCATCAAAGAAGATCGTGGTACTGAAATCGCCCTGCATCTGCGTGAAGGCGAATACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCTCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATcaacaa                                                                                                                    <  1:494545/251‑1 (MQ=255)                                                                                                                         aaaaCCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGGGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAg                                                                                                                  <  2:528400/251‑1 (MQ=255)                                                                                                                                    aTGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCtg                                                                                                       >  2:413072/1‑251 (MQ=255)                                                                                                                                      ggCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCtgt                                                                                                      <  1:172805/250‑1 (MQ=255)                                                                                                                                        cgtcgtCTGGGAATCGGCTGGCGAAGGTGATTACACCGTTGCCGACATCAGCAAAGAAGATCGTGTTACTGAAATCACCCTGCGTCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCAGCATCAGCAAATACTCCGACCATATCGCTCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAATCCGTTATCTCCTGGGAGAAAATCAACAAAGCGCATGCGCTGTGGa                                                                                                   <  1:271564/251‑1 (MQ=255)                                                                                                                                          tcttctGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACt                                                                                                 >  2:447758/1‑251 (MQ=255)                                                                                                                                             tctGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTc                                                                                                <  2:305777/249‑1 (MQ=255)                                                                                                                                                gggAATCGGCTGGGGAAAGGGAATACACCGGTTCCGACATCAACAAAGAAGATCGTGGTCATGAAAACACCCTTAAGCTGCGTGAAGGCGAAGACGAGGTTCC‑CGATGACTGGCGCGTGCGGTCCATTATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGt                                                                                              <  2:415414/248‑1 (MQ=255)                                                                                                                                                gggAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGAACATATCGCGCTGCCGGTAGAGATCGAAAAACGCGATGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCCGGAGCTGTGGACTCGTAAc                                                                                           >  2:113398/1‑251 (MQ=255)                                                                                                                                                              gAAGGTGAATACACCGTTGCCGACATCATCAAAGAAAATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAAAGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGCGAAAATCAACAAAGAGCAGGCGCTGGGGACTCGTAAAAAGTCGGCAATccc                                                                             >  1:295905/1‑249 (MQ=255)                                                                                                                                                              gAAGGGGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGACAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCAc                                                                             >  1:110161/1‑251 (MQ=255)                                                                                                                                                                aGGTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCg                                                                           >  2:244599/1‑251 (MQ=255)                                                                                                                                                                  gTGAATACACCGTTGCCGACATCATCAAAGAAGATCGTGGTACTGAAATCCCCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGAt                                                                         <  2:444470/251‑1 (MQ=255)                                                                                                                                                                                           cAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCTATGAAGAGTACAAAGAGTTCTACAAac                                                <  2:366921/251‑1 (MQ=255)                                                                                                                                                                                           cAAAGAAAATCGTGGTACTGAAATCACCCTGCATCTGCGTGCAGGCGAAGACGAG‑TTTCTCGCTGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGCAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGCAGAGGACAAAGAGTTCTACAAac                                                >  2:77605/1‑251 (MQ=255)                                                                                                                                                                                                 agaTCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCg                                          <  1:447758/251‑1 (MQ=255)                                                                                                                                                                                                    tCGTGGTACTGAAATCACCCTGCATCTGCGTGAATGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGccc                                       >  2:65375/1‑251 (MQ=255)                                                                                                                                                                                                           aCTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACt                                 <  2:10827/250‑1 (MQ=255)                                                                                                                                                                                                                   cACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGACAATCACCAAAGCGCAGGCGCTGTGGACTCGTAAAAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACCCATCGCCCACGACTTTAATGATc                        >  1:501323/1‑251 (MQ=255)                                                                                                                                                                                                                   cACCCTGCATCGGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCAGCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATc                        <  2:310879/251‑1 (MQ=255)                                                                                                                                                                                                                     cccTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGAGTGCGTTCCATCATCAGCAAATACTACGACCATATCGCGCTGCCGGTAGAGATCGAAAAAAGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAATTCTACAAACACATCGCCCACGCCTTTAATGATCCg                      >  1:55083/1‑251 (MQ=255)                                                                                                                                                                                                                      ccTGCATCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGc                     >  2:133870/1‑251 (MQ=255)                                                                                                                                                                                                                            tCTGCGTGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCt               >  2:80340/1‑251 (MQ=255)                                                                                                                                                                                                                               gCGTGAAGGCGAAGACGAG‑TGTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCGCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGAcc                <  1:62638/247‑1 (MQ=255)                                                                                                                                                                                                                                 gTGAACGCGAAGACGAG‑TTTCTCGATGAATGGCGCGTGCGTTCCATCATCAGCAAATCCTCCGCCCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATATCCTGGGAGAAAATCACCAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGCAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGc          >  2:326204/1‑251 (MQ=255)                                                                                                                                                                                                                                  tGAAGGCGAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGcc         >  1:179450/1‑251 (MQ=255)                                                                                                                                                                                                                                         gAAGACGAG‑TTTCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAACCg  >  1:111671/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                CGAATCCCTGGGTTCTGACCAGGCGAAAGACAGCCAGCTGATCGGTCAGTTTGGTGTTGGTTTCTACTCTGCGTTTATCGTGGCCGACAAAGTGACCGTGCGTACTCGCGCGGCAGGCGAAAAACCAGAAAATGGCGTCTTCTGGGAATCGGCTGGCGAAGGTGAATACACCGTTGCCGACATCACCAAAGAAGATCGTGGTACTGAAATCACCCTGCATCTGCGTGAAGGCGAAGACGAG‑TTCCTCGATGACTGGCGCGTGCGTTCCATCATCAGCAAATACTCCGACCATATCGCGCTGCCGGTAGAGATCGAAAAACGCGAAGAGAAAGACGGCGAAACCGTTATCTCCTGGGAGAAAATCAACAAAGCGCAGGCGCTGTGGACTCGTAACAAGTCGGAAATCACCGATGAAGAGTACAAAGAGTTCTACAAACACATCGCCCACGACTTTAATGATCCGCTGACCTGGAGCCACAACCG  >  NC_000913/495437‑495919

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