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	98.3%	181.1 / NA	58	P26P (CCC→CCG)	ileS	isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (34/23);  minor base A (0/1);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAG  >  NC_000913/22222‑22711                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    aCAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGaa                                                                                                                                                                                                                                                 <  2:260503/251‑1 (MQ=255)    agtagtGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGc                                                                                                                                                                                                                                              <  1:290335/251‑1 (MQ=255)      gaatcctgcTTAAAAAAATACGCAGTAAGCAGCGATTTGCGTAGCTGGACGAAATGAAAGCGCAAAGTGCACGTGATGAAATAACCGCCCGCGAATTTGTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACAGGGAACCGAGAATCTGATGAGTGACTATAAATCATCCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCAGGAATGCTg                                                                                                                                                                                                                                            <  2:286780/246‑1 (MQ=255)           ctgcGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGcgcg                                                                                                                                                                                                                                       >  1:493253/1‑251 (MQ=255)                   aaaaaaACGCAATGAGCAGCGGTTTGCGCCGCTGAACGAACTGAAAGCGCAGATTGCGCGTGATGGATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGAGTTGGACtg                                                                                                                                                                                                                               >  2:565166/1‑251 (MQ=255)                   aaaaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtg                                                                                                                                                                                                                               <  1:276564/251‑1 (MQ=255)                       aTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatga                                                                                                                                                                                                                           <  1:354637/251‑1 (MQ=255)                         aCGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAAACCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACtgatgatg                                                                                                                                                                                                                         >  2:397569/1‑251 (MQ=255)                                gAGCAGCGATTTGCGTCGCTGGACGCACTGAAAGCGACGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTCTGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGAGCGAACCGGGACTGCTGGCGCGTTGGACTGCTGATGATCTGt                                                                                                                                                                                                                   >  2:407826/1‑250 (MQ=255)                                           tgagTCGCTGGCCGAACTTAAAGCGCAGATTGCGCGTTATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCatcatc                                                                                                                                                                                                        <  2:18065/247‑1 (MQ=255)                                                cgatGGACGAACTGAAAGCGCAGATTGCGCGAGATGAATTAACCGCCCGCGAATTTTTTAGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGACTTTGCCGAAAACAGGGTTCCCGATGCGTGTCTATCTCGCCACGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATTATCTTTACGGCATCATCCGTGCg                                                                                                                                                                                                  >  2:462136/4‑251 (MQ=255)                                                   tGGACGAACTGTAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGACTTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGCGAATCTGATGAGTGCCTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCt                                                                                                                                                                                               >  2:325732/1‑251 (MQ=255)                                                    ggACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGc                                                                                                                                                                                                <  1:39631/249‑1 (MQ=255)                                                     gACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTa                                                                                                                                                                                              >  1:196088/1‑250 (MQ=255)                                                      aCGAACTGAAAGCGCTGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGTCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAAGCAACCCTGAATTTGCCGGAAATAGGGTTACCGATGCGGGGCGACCCCGCCAAGCGGGAACCGGGAACGCGGGCGCGTGGGACTGAGGATGATCTGTACGGCATCACTCGTGCGGCTaaa                                                                                                                                                                                            <  2:614913/251‑1 (MQ=255)                                                             gAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGCATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGc                                                                                                                                                                                      >  1:32051/1‑250 (MQ=255)                                                             gAAAGCGCAGATTGCGCGTGATGAATTAACAGCCCGCGAATTTTTTGGGCTAACAACACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCATGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGc                                                                                                                                                                                      >  1:557811/1‑250 (MQ=255)                                                              aaaGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCa                                                                                                                                                                                     >  1:112342/1‑250 (MQ=255)                                                               aaGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCaa                                                                                                                                                                                    >  1:554911/1‑250 (MQ=255)                                                                  cgcAGAGTGCTCGTTATGAATTAACCGCCCGCGAATTTGTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAACCACGTAACCGAGAATCTGATGAGGGACTATAAAGCAACCCTGAATATGCCGGAAACAGGGTTCCCGAAGCGTGGCGATCGCTCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGAGGATGATCGGAACGTCATGATCGGTGCGGCGAAAAAAGTCaaaaaa                                                                                                                                                                                <  2:625416/251‑1 (MQ=255)                                                                    cAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAAc                                                                                                                                                                              >  1:93980/1‑251 (MQ=255)                                                                      gATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCt                                                                                                                                                                            <  2:597050/251‑1 (MQ=255)                                                                      gATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCt                                                                                                                                                                            >  1:597050/1‑251 (MQ=255)                                                                         tGCGCGTGACGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGATTAAGCCTGTTATGTAATCAAACCGAAACACGGCACCGGGGATCTGATGAGTGACTATGAATCAACCCTGAATATGCCGGCACAAGGGTTCACGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGCACGCTGGCGCAATGGACTGATGATGATCAGTACGGCATCATCCGTGGGGATAACAAAGGCACAAACACCTTc                                                                                                                                                                          >  1:611131/1‑250 (MQ=255)                                                                         ggcgcGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCTGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGAGCTCGCCAAGCGCGAACCGGGAATGCTGGCGCTTTGGACTGATGATGATCTGTACGGCATCAGCCGTGCGGCTAAGAAAGGCAAAAAAACCTTCa                                                                                                                                                                         <  2:64003/250‑1 (MQ=255)                                                                           cgcgTGATGAATTAACCGCCCGCGACTTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCAtt                                                                                                                                                                       >  2:35319/1‑251 (MQ=255)                                                                             ggTGATGAATTAACCGCCCGCGAATGTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCt                                                                                                                                                                     <  2:521275/250‑1 (MQ=255)                                                                                   gaATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAA‑TCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCatgat                                                                                                                                                              >  1:532381/1‑251 (MQ=255)                                                                                      ttAACCGCCCGCGAATTTTTGGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGCAACACGGAACCGAGACTCTGATGAGTGACTATAAATCAACACTGAATTTGACGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGCCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATgg                                                                                                                                                            >  2:29104/1‑251 (MQ=255)                                                                                           cgcccgcGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCACTCTGCATGATGGCcctc                                                                                                                                                       >  1:549789/1‑251 (MQ=255)                                                                                              ccgcGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCtt                                                                                                                                                    >  1:457221/1‑251 (MQ=255)                                                                                                              cTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCAAATGGCAGCAt                                                                                                                                    >  2:306190/1‑251 (MQ=255)                                                                                                               tAACAAAGCCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAACTCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGGACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCAtt                                                                                                                                   >  1:445425/1‑251 (MQ=255)                                                                                                                     aaCCGGCTTAAGCCTGTTATGTAATCAAACCTAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCTAATGGCAGCATTCATAtt                                                                                                                             <  1:537890/251‑1 (MQ=255)                                                                                                                      aCCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTg                                                                                                                            >  1:347936/1‑251 (MQ=255)                                                                                                                        cGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGAGGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGt                                                                                                                          <  1:365701/251‑1 (MQ=255)                                                                                                                           cTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCAc                                                                                                                       >  1:218403/1‑251 (MQ=255)                                                                                                                              aaGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCCTTCATATTGGTCACTCg                                                                                                                    >  2:489866/1‑251 (MQ=255)                                                                                                                                  cTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGtt                                                                                                                 >  2:202824/1‑250 (MQ=255)                                                                                                                                                      aacaCGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACAtt                                                                                            <  2:240790/248‑1 (MQ=255)                                                                                                                                                        caCGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTAt                                                                                          >  1:556934/2‑251 (MQ=255)                                                                                                                                                                 gagaATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTc                                                                                 >  1:109447/1‑251 (MQ=255)                                                                                                                                                                            gaGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGAGCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGTCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCttt                                                                      <  2:32051/251‑1 (MQ=255)                                                                                                                                                                               tGACTATAAATCAACCCTGAATTTGAAGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGTTGGCGCGGTGGACTGATGCTGATCTGTACGTCATCATCCGTGCGGCTAAAAAACGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGTTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAAGGCTTTCCg                                                                   >  1:381897/1‑251 (MQ=255)                                                                                                                                                                                         tCAACCCTGCATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCg                                                                   >  1:532941/1‑241 (MQ=255)                                                                                                                                                                                         tCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCg                                                                   <  2:532941/241‑1 (MQ=255)                                                                                                                                                                                                  aaTTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTAGCCAAGCGCGAACCGTGAATGATGGCGCGTTCGACTGATGATGATCTGTACGGCATCATCCGTGCGTCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTAtgt                                                <  1:515485/251‑1 (MQ=255)                                                                                                                                                                                                          ggAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCctggct                                        >  1:399481/1‑251 (MQ=255)                                                                                                                                                                                                            aaaCAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGAATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGCCATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCctggctgg                                      >  1:300669/1‑251 (MQ=255)                                                                                                                                                                                                                  ggTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGc                                <  2:112342/251‑1 (MQ=255)                                                                                                                                                                                                                     tCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGGTGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCAc                             >  1:22921/1‑251 (MQ=255)                                                                                                                                                                                                                     tCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCAc                             >  1:59314/1‑251 (MQ=255)                                                                                                                                                                                                                       ccGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCg                                                     >  2:13705/1‑225 (MQ=255)                                                                                                                                                                                                                       ccGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCg                                                     <  1:13705/225‑1 (MQ=255)                                                                                                                                                                                                                                        cGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCt            <  2:150341/249‑1 (MQ=255)                                                                                                                                                                                                                                        cGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGAGCGAGCTGa          <  1:489866/251‑1 (MQ=255)                                                                                                                                                                                                                                               ggcgAACCGGGAATGCTGGCGCGTTGGACTGATGAGGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCGTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTAGGACTCGCCGTATGCGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGa   <  2:259232/250‑1 (MQ=255)                                                                                                                                                                                                                                                gcgAACCGGGAATGCTGGCGCGTGGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATGCTGCATGATGGCCCGCCTTATGCGAATGGCAGCATGCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTAGCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCGCGGTCTGCCGAGCGAGCTGAAATTCGAg  <  2:317450/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    ACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAGATTCTGAAAGACATTATCGTGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAG  >  NC_000913/22222‑22711

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