Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	22,468	C→G	P26P (CCC→CCG)	ileS →	isoleucine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	22,468	0	C	G	97.5%	124.0 / ‑6.8	40	P26P (CCC→CCG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (24/15);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTG  >  NC_000913/22230‑22702                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            ggctgcGTAAAAAAATACGCAATGAGCAGCTATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGgcg                                                                                                                                                                                                                                 <  1:375852/250‑1 (MQ=255)  gctgcGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAATCGCAGATTGCGCGTGATGAATTAACCGCCCTCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGgcg                                                                                                                                                                                                                                 <  1:353916/250‑1 (MQ=255)             aaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGCATCTGATGAGTGACTATCAATCAACCCTGAATTTGCCGGAAAAAGGGTTCCCGATGAGTGGCGATCTCGCCACGCGCGAACAGGGAATGCTGGCGCGTTGGACtgat                                                                                                                                                                                                                     >  1:226762/1‑251 (MQ=255)                tACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGACAGCGCCGATTGCGCGTGATGACTTCACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGCCTCTAAATCAACCCTGAATTTGCCGGAAACAGGGTCCCCGAGGCGTGGCGATCCCGCCAAACGCGAACCGGGAATGCTGGCGCGGTGGACagatgat                                                                                                                                                                                                                  >  2:180126/1‑251 (MQ=255)                  cGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCCGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGCACCGGGAATGCTGGCGCGTTGGACtgatgatgg                                                                                                                                                                                                                >  1:135052/1‑250 (MQ=255)                     aaTGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCt                                                                                                                                                                                                             <  1:230387/251‑1 (MQ=255)                            agcGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGc                                                                                                                                                                                                      >  1:32881/1‑251 (MQ=255)                                    gCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGcatcagcc                                                                                                                                                                                               <  1:36023/250‑4 (MQ=255)                                     cGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCgg                                                                                                                                                                                             >  1:158613/1‑250 (MQ=255)                                       ccGCTGGACGAACTTAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGc                                                                                                                                                                                           <  2:12548/250‑1 (MQ=255)                                                      aaaGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTGGTGGGCTAACAAAACCGGCTTAAGACTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTTATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCaa                                                                                                                                                                            <  1:37063/251‑1 (MQ=255)                                                        aGCGCAGATTGCGCGAGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACAGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAAAAGGGTTCCCGATGCGAGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGACGCGTTGGAGTGAGGATGATCTGTACGGCATCAGCCGCGCGGCGAAAAAGGGCaac                                                                                                                                                                           >  2:278912/1‑249 (MQ=255)                                                             aGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCACTCTCGCCAAGCGCGAACCGGGAATGCTGGAGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGAGGCTAAAAAAGGCAAAAAAAcc                                                                                                                                                                     >  2:363425/1‑251 (MQ=255)                                                               aTTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑tt                                                                                                                                                                  >  1:371085/1‑251 (MQ=255)                                                                          tgaATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCatg                                                                                                                                                       >  1:86392/1‑251 (MQ=255)                                                                                 aCCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGccc                                                                                                                                                <  1:227914/251‑1 (MQ=255)                                                                                     cccgaGAATTTTTTGGGCTAACAAAAAAGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGCTACCGATTCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGAGCGTTGGACTGCTGATGATCGGTACGTCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTATCCATGATGGccatcat                                                                                                                                            >  2:270173/1‑246 (MQ=255)                                                                                       cgcgAATTTTTTGGGCTAACAAAACCGTCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAAGCTGATGAGTGACTATAAATCAACCCGGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCTAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGt                                                                                                                                                                                             <  2:216779/201‑1 (MQ=255)                                                                                       cgcgAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGt                                                                                                                                                                                             >  1:216779/1‑201 (MQ=255)                                                                                               ttttGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGt                                                                                                                                                                                                           <  1:269147/179‑1 (MQ=255)                                                                                               ttttGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACAGGGAATGCTGGCGCGTTGGACTGATGATGATCTGt                                                                                                                                                                                                           >  2:269147/1‑179 (MQ=255)                                                                                                           aaaaacagtttAAGCCTGTTATGAAATCAAACCGAAACAAGGAACCGCGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCAGGAAACAGGGTTCCCGGTGAGTGGCGATCTCGCAAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGCAGATGATCTGTACGGCCTCATCCGTGCGGCTAAAAAAGG‑AACACAACCATTCATTCTGCATGATGGCCCTCATTCTGCGAATGGCAGCATTCat                                                                                                                       >  2:337960/10‑250 (MQ=255)                                                                                                            aaaCCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCctat                                                                                                                     >  1:248516/1‑251 (MQ=255)                                                                                                               ccGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCAGATGCGTGGCGATCTCGCAAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGATAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGAAGGCCCTCATTAAGCGAATGGCAGCATTCATATTgg                                                                                                                  >  1:23433/1‑251 (MQ=255)                                                                                                                      aaGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTc                                                                                                            <  1:363425/250‑1 (MQ=255)                                                                                                                      aaGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCg                                                                                                           >  1:94353/1‑251 (MQ=255)                                                                                                                            gTTATGTAATCAATCCGTAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACATGGTTCCCTATGCGTGGCGATCTCTCCAAGCGCGAACCGGGAATGCTGTCTCGTTGGACTGATTATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTaa                                                                                                      <  1:265645/250‑1 (MQ=255)                                                                                                                              tATGTAATCAAACCGAAACACGGAACCGAGCATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAACGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACa                                                                                                    >  1:1/1‑250 (MQ=255)                                                                                                                                tGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACa                                                                                                    >  1:308637/1‑248 (MQ=255)                                                                                                                                 gTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCat                                                                                                                       >  1:79663/1‑228 (MQ=255)                                                                                                                                 gTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAAGTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCat                                                                                                                       <  2:79663/228‑1 (MQ=255)                                                                                                                                      cAAACCGAAACACGGAACCGAGACGCTGATTAGTGACTATAAATCAACCCTGAAGTGGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGTTAAAAAAGTCAAAAAAACC‑TGCATTCTGCATGATGTTCCTCATTATGCGAATGGCAGCATTCATATAGGTCACTCGGTTAACAAGATTCTGa                                                                                           <  2:192566/251‑1 (MQ=255)                                                                                                                                                   ggAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCACGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGt                                                                              >  1:369210/1‑251 (MQ=255)                                                                                                                                                         gagaATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTc                                                                        <  1:270173/251‑1 (MQ=255)                                                                                                                                                                       tGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATATCACCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGCACGGCATCATCCGTGCGGCTAAAAAAGGCAAACAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGTTTCTGAAAGACATTATCGCGACGTCAAAAGGGCTTTCCg                                                          >  2:294902/1‑251 (MQ=255)                                                                                                                                                                                   aaCCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGc                                              >  1:134700/1‑251 (MQ=255)                                                                                                                                                                                      ccTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCACAAAAACC‑TTCATTCTGCATGATGACCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAATTCCAAAGGGCTTTCCGGTTATGACTCGCCGt                                           >  2:152190/1‑251 (MQ=255)                                                                                                                                                                                          aaTTTTCCGGAAACCAGGTTCCCGATTCGTGGCGATCTCGCCAAGCGCGACCCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTAtgt                                       <  1:129045/251‑1 (MQ=255)                                                                                                                                                                                                                  aTGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCt               >  1:377527/1‑251 (MQ=255)                                                                                                                                                                                                                                cgacAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTACCTGGCTGGGACTGCCACGGTCTGCCGATCGAgatg  >  2:294104/4‑247 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            TGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACC‑TTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTG  >  NC_000913/22230‑22702

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: ‑