Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,371,401	C→T	A126V (GCT→GTT)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,371,401	0	C	T	98.9%	307.6 / NA	93	A126V (GCT→GTT)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (44/48);  minor base A (1/0);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACA  >  NC_000913/4371163‑4371643                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    ccGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTgaaga                                                                                                                                                                                                                                         <  1:564808/251‑1 (MQ=255)  cGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGCAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTgaaga                                                                                                                                                                                                                                         >  1:224909/1‑250 (MQ=255)    aCCATCACCACAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACt                                                                                                                                                                                                                                      >  1:504335/1‑251 (MQ=255)    aCCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGACAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCACCCGTCCTGGCTCAGGCTATCATCACTGAAGGCCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGCATCGACAAAGCGGTTACCGTTGCCGTTGCAGAACt                                                                                                                                                                                                                                      >  2:482273/1‑251 (MQ=255)    aCCATCACCAAAGATGGTGGTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTACAGTTGAAGAACt                                                                                                                                                                                                                                      <  2:38913/251‑1 (MQ=255)     ccATCACCAAAGAGGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTg                                                                                                                                                                                                                                     <  2:421251/251‑1 (MQ=255)        tCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGaa                                                                                                                                                                                                                                   >  1:249892/1‑250 (MQ=255)              aaGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGCTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTg                                                                                                                                                                                                                            <  2:472788/251‑1 (MQ=255)                      gtTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGCAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGATGTCCGTAcc                                                                                                                                                                                                                    >  1:302575/1‑251 (MQ=255)                          cagTTGCTCGTGAAATAGACCTGGAAGACAAGTTCGAACATATGGGTGCGCCGCTGGTGAAAGACGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCTGGCTATCATCACTGAAGGGCTGAAAGCTGTTGCTGCGGGCATGAAACCGCTGGACCTGAAACGTGGTATCGACAAATCGGTTACCGTTGCAGTTGACGAACTGAAAGCGCTGTCCGTACCATGc                                                                                                                                                                                                                >  2:71496/3‑251 (MQ=255)                            gTTGCTCGTGAAATCGAACGGGAAGACAAGGTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGATGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTCACATGctc                                                                                                                                                                                                              <  2:81065/251‑1 (MQ=255)                                cTCGTGAAATCGAACTGGAAGACAAGTTCGAAACTATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGCAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAACCGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGAc                                                                                                                                                                                                          >  2:40943/1‑251 (MQ=255)                                     gAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGt                                                                                                                                                                                                                           >  1:263503/1‑229 (MQ=255)                                     gAAATCGAACAGGAAGACAAGTTCGAAAAGATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAATCATGCGACGCTTTAGGCGACGGTACCACCACTGCAATCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGCCAAAGCGGTTACCGTTTCAGTTGAAGAACTGAAAGCGATGTCCGTTCCATGCTCTGACTCTaa                                                                                                                                                                                                     <  2:485298/251‑1 (MQ=255)                                     gAAAGCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACGGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAgcgcggt                                                                                                                                                                                                                           <  2:263503/229‑4 (MQ=255)                                      aaaTCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTaa                                                                                                                                                                                                     >  2:207106/1‑250 (MQ=255)                                       aaTCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTaa                                                                                                                                                                                                     >  1:12734/1‑249 (MQ=255)                                          cGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCg                                                                                                                                                                                                 <  2:422923/250‑1 (MQ=255)                                          cGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGa                                                                                                                                                                                                <  1:102240/251‑1 (MQ=255)                                           gAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGAt                                                                                                                                                                                               >  2:395152/1‑251 (MQ=255)                                            aaCTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGAtt                                                                                                                                                                                              <  1:173452/251‑1 (MQ=255)                                                ggAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGcg                                                                                                                                                                                           <  2:443418/250‑2 (MQ=255)                                                ggAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTc                                                                                                                                                                                          <  1:135683/251‑1 (MQ=255)                                                       aaGTTCGAAAATAGGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTGGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTgg                                                                                                                                                                                   <  2:354778/251‑1 (MQ=255)                                                        aGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTAGCAGCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGt                                                                                                                                                                                  <  2:259739/251‑1 (MQ=255)                                                            cGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCa                                                                                                                                                                              <  1:129911/251‑1 (MQ=255)                                                                    tGGGTGCGCCGCTGGTGAAAGAAGTTGCCTCTACAGCAAACGACGATGCAGGCGACGGTCCCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGCAGGTCTGAAAGCTGTTGCTGCGGGCATGAAACCGATGGACCTGCAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGCCTCTAAAGCGATTGCTCCGGTTGGTACCCTCTCCGCt                                                                                                                                                                      >  2:526325/1‑251 (MQ=255)                                                                        tGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACt                                                                                                                                                                  >  1:104072/1‑251 (MQ=255)                                                                          cgcAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTcc                                                                                                                                                                <  1:112480/251‑1 (MQ=255)                                                                          cgcAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTcc                                                                                                                                                                <  1:500406/251‑1 (MQ=255)                                                                            cAGATGGTGAAAGAAGTTGACTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGACTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGAATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGAGATTGCTCAGGTTGGAACCATCTCCGCTAACTCcg                                                                                                                                                               >  2:430378/1‑250 (MQ=255)                                                                             aGATGGTGAAAGCAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGAGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGATAACTCcgac                                                                                                                                                             >  1:510880/1‑251 (MQ=255)                                                                               aTGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCcgacg                                                                                                                                                            >  2:436414/1‑250 (MQ=255)                                                                                    gAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCg                                                                                                                                                      <  1:94894/251‑1 (MQ=255)                                                                                     aaaGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGt                                                                                                                                                     <  1:68701/251‑1 (MQ=255)                                                                                      aagaagTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAAACGt                                                                                                                                                     >  1:322545/1‑250 (MQ=255)                                                                                      aagaagTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGGTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCCCCGCTAACTCCGACGAAACCGt                                                                                                                                                     >  2:370211/1‑250 (MQ=255)                                                                                      aagaagTTGCCTCTAAAGCAAACCACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTa                                                                                                                                                    <  1:348666/251‑1 (MQ=255)                                                                                       agaagTTGCCACTAAAGCAAACGACGCTGCGGGCGACGGTACCACCACTGCAACCTTACTGGCGCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTGGAAGAACTGAAAGCGCGGGCGGTACCATGCTCTGACTCTAAAGCGAGTGCTCAGGTTTGTCCCCTCTCCGCTAACTCCGACGAAACCGTAg                                                                                                                                                   <  1:376734/251‑1 (MQ=255)                                                                                              gCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGCCAAAGCGGTTACCGTTGCAGTTGAAGAACTGCAAGCGCTGTACGTACCATGCTCTGACTCTAAAGCGATTGATCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑t                                                                                                                                           >  1:230183/1‑251 (MQ=255)                                                                                              gCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑t                                                                                                                                           >  2:532347/1‑251 (MQ=255)                                                                                               ccTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑Tg                                                                                                                                          >  2:142212/1‑251 (MQ=255)                                                                                                      aGCAAACGACGCTGCAGGCGACGGTACCACCACTGCATCCGTACTGGCTCAGGCTATCATCACTGAAGGACTGAAAGCTGTTGCTGCGGGCAGGAACCCGAGGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCAGGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑Tgcacgctg                                                                                                                                   <  2:552518/251‑8 (MQ=255)                                                                                                      aGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTg                                                                                                                                   <  2:552540/251‑1 (MQ=255)                                                                                                       gCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGa                                                                                                                                  <  1:184939/251‑1 (MQ=255)                                                                                                        cAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACGGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCAGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCGCAGGTTGGTACCATCTCCGCTAACTCCGACGAATCCGttgggtaactgaaccctgaa                                                                                                                                  <  2:537191/251‑20 (MQ=255)                                                                                                          aaCGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGc                                                                                                                               >  1:497545/1‑251 (MQ=255)                                                                                                           acgacgCTGCCGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGAGCTGTCAGTCCCATGCTCTGACTCTAAAGTGATTGATACGGTAGGGAACAATTACGATAACTCCGACGAAACCGTACGGAAAC‑TGATCgaggaag                                                                                                                                >  2:310360/1‑243 (MQ=255)                                                                                                           acgacgCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCg                                                                                                                              >  2:310375/1‑251 (MQ=255)                                                                                                               cgCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATgg                                                                                                                          >  1:583454/1‑251 (MQ=255)                                                                                                                   gCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGAc                                                                                                                        <  1:160275/249‑1 (MQ=255)                                                                                                                        cGACGGTCCCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCGGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAA‑TGATCGCTGAAGCGATGGACAAAGTCg                                                                                                                 <  2:428962/251‑1 (MQ=255)                                                                                                                            ggTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTaa                                                                                                             <  1:207106/251‑1 (MQ=255)                                                                                                                            ggTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTa                                                                                                              >  2:474058/1‑250 (MQ=255)                                                                                                                             gTACCTCCACTGCAATCGTACTGGCCCAGGCTATCATCACTGAAGGTCTTAAAGCTGTTGCTGCGGGCTTGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAATC‑TGATCGCTGAAGCGATGGACAAAGTCGGTaaa                                                                                                            <  1:520519/251‑1 (MQ=255)                                                                                                                               accaccACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTaa                                                                                                             >  1:62947/1‑248 (MQ=255)                                                                                                                                    cacTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTGCCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCg                                                                                                     <  1:246296/251‑1 (MQ=255)                                                                                                                                       tGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGtt                                                                                                   >  2:124855/1‑250 (MQ=255)                                                                                                                                        gCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGtt                                                                                                   <  1:24478/249‑1 (MQ=255)                                                                                                                                               tACTGGCTCAGGCTATCATCAGTGAAGGTCTGAAAGCTGTTGCTGAGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGAGGCTACCGATGCAGGTGCAGAACTGAAAGCGCTGTCAGTACCATGCGCTGACTCTAAAGCGATCGCTCAGGCTGGTACCATCTCCGCTAATTCCGACGACACCGTAGG‑AAACATGATAGCAGAAGCGATGGACAAAGTCGGCAaagaggctgctatctacgtt                                                                                          >  2:464258/1‑235 (MQ=255)                                                                                                                                                    gCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGa                                                                                        >  1:367261/1‑248 (MQ=255)                                                                                                                                                      ccAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACTGTTGAAGACg                                                                                   <  2:46618/250‑1 (MQ=255)                                                                                                                                                        aGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGt                                                                                 >  1:246464/1‑251 (MQ=255)                                                                                                                                                             atcatcACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTAcc                                                                              <  1:434699/249‑1 (MQ=255)                                                                                                                                                                         ggTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGa                                                                <  1:124855/251‑1 (MQ=255)                                                                                                                                                                              gAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGCCAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACt                                                             <  2:419444/249‑1 (MQ=255)                                                                                                                                                                              gAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTgg                                                           >  1:147329/1‑251 (MQ=255)                                                                                                                                                                                 aGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACg                                                        <  1:263820/251‑1 (MQ=255)                                                                                                                                                                                  gCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGt                                                       <  1:395152/251‑1 (MQ=255)                                                                                                                                                                                       tgctgcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGtt                                                   <  2:412759/250‑1 (MQ=255)                                                                                                                                                                                        gctgcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTg                                                  <  1:265316/250‑1 (MQ=255)                                                                                                                                                                                           gcGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAgg                                              >  1:520815/1‑251 (MQ=255)                                                                                                                                                                                                 aTGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGc                                         >  1:99862/1‑250 (MQ=255)                                                                                                                                                                                                 aTGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGATGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGc                                         >  1:277247/1‑250 (MQ=255)                                                                                                                                                                                                    aaCCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGtt                                     <  2:246464/251‑1 (MQ=255)                                                                                                                                                                                                        cGATGGACCTGAAACGTGGTATCGACAAAGCGGATACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGAAAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGCCGGTACCGGTCTGCAGGACGAACTGGCCGTGGTTGAAGGTATGCAGTTCGAc                                 >  1:523582/1‑251 (MQ=255)                                                                                                                                                                                                         gATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTTGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAATAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGAcc                                <  1:257371/251‑1 (MQ=255)                                                                                                                                                                                                            ggACCTGAAACGGGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTg                             <  2:497545/251‑1 (MQ=255)                                                                                                                                                                                                                   aaaCGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCt                      >  2:500220/1‑251 (MQ=255)                                                                                                                                                                                                                   aaaCGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCt                      <  2:126698/251‑1 (MQ=255)                                                                                                                                                                                                                     aCGTGGTATCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGt                    >  1:195548/1‑251 (MQ=255)                                                                                                                                                                                                                             tCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAGCCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGGCGGTAAAGAAGGCGTTATCACCGt                                                                                           <  2:347005/172‑1 (MQ=255)                                                                                                                                                                                                                             tCGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTACAC‑TGATCGCTGAAGAGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGt                                                                                           >  1:347005/1‑172 (MQ=255)                                                                                                                                                                                                                              cGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTAc            <  1:264495/250‑1 (MQ=255)                                                                                                                                                                                                                              cGACAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAAACGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCAcc                                                                                             >  2:375079/1‑169 (MQ=255)                                                                                                                                                                                                                              cGACAAAGCGGTGACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCAcc                                                                                             <  1:375079/169‑1 (MQ=255)                                                                                                                                                                                                                                aCAAAGCGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTtc         <  1:84734/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGGTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCCCCGCTAACTCCGAAGAAACCGTAGGTACAC‑TGATCCCTGCAGCGCTGGACACAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGCCGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGCATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTtcatcacc   >  1:582704/1‑249 (MQ=255)                                                                                                                                                                                                                                       ggTTACCGTTGCAGTTGACGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTAAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTtca        >  2:235899/1‑245 (MQ=255)                                                                                                                                                                                                                                       ggTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaac   >  1:453223/1‑250 (MQ=255)                                                                                                                                                                                                                                       ggTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGGGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTtca        <  1:235899/245‑1 (MQ=255)                                                                                                                                                                                                                                        gTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTACCTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaaca  >  2:415382/1‑250 (MQ=255)                                                                                                                                                                                                                                        gTTACCGTTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATcaaca  >  2:339658/1‑250 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    CCGACCATCACCAAAGATGGTGTTTCCGTTGCTCGTGAAATCGAACTGGAAGACAAGTTCGAAAATATGGGTGCGCAGATGGTGAAAGAAGTTGCCTCTAAAGCAAACGACGCTGCAGGCGACGGTACCACCACTGCAACCGTACTGGCTCAGGCTATCATCACTGAAGGTCTGAAAGCTGTTGCTGCGGGCATGAACCCGATGGACCTGAAACGTGGTATCGACAAAGCGGTTACCGCTGCAGTTGAAGAACTGAAAGCGCTGTCCGTACCATGCTCTGACTCTAAAGCGATTGCTCAGGTTGGTACCATCTCCGCTAACTCCGACGAAACCGTAGGTAAAC‑TGATCGCTGAAGCGATGGACAAAGTCGGTAAAGAAGGCGTTATCACCGTTGAAGACGGTACCGGTCTGCAGGACGAACTGGACGTGGTTGAAGGTATGCAGTTCGACCGTGGCTACCTGTCTCCTTACTTCATCAACA  >  NC_000913/4371163‑4371643

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