Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,634,927	+A	intergenic (‑251/+438)	unknown ← / ← unknown	Permease of the drug/metabolite transporter (DMT) superfamily/LSU rRNA ## 23S rRNA, large subunit ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,634,925	1	.	A	100.0%	132.5 / NA	52	intergenic (‑249/+440)	unknown/unknown	Permease of the drug/metabolite transporter (DMT) superfamily/LSU rRNA ## 23S rRNA, large subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (30/22);  total (30/22)

CCCAGAGTAAATACCATTCCTAAATGCAATTTTTAAATATCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTAT‑AACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTTCTCTCGTTCTTT  >  1/1634747‑1635157                                                                                                                                                                                    |                                                                                                                                                                                                                                         ctcctgaccgacgataccgtaaaaTGCAATTTTTAAATATCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGc                                                                                                                                                                                                          <  1:158135/191‑1 (MQ=255) ctcctgaccgacgataccgtaaaaTGCAATTTTTAAATATCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCgg                                                                                                                                                                                       >  1:493680/22‑231 (MQ=39)               taccgtaaaaTGCAATTTTTAAATATCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTAc                                                                                                                                                     <  1:739448/244‑1 (MQ=255)                accgtaaaaTGCAATTTTTAAATATCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGccc                                                                                                                                                             <  1:14958/236‑1 (MQ=255)                 ccgtaaaaTGCAATTTTTAAATATCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCa                                                                                                                                                   >  1:410561/6‑251 (MQ=255)                                     atatCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCAc                                                                                                                               <  1:505821/251‑1 (MQ=255)                                         ccTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTAc                                                                                                                                                     >  1:580333/1‑225 (MQ=255)                                         ccTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACtt                                                                                                                             >  1:875140/1‑249 (MQ=255)                                                gTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATgg                                                                                                                                                                     <  1:253371/202‑1 (MQ=255)                                                   gATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTa                                                                                                                 >  1:543224/1‑251 (MQ=255)                                                    aTCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCacac                                                                                                                                                        >  1:655442/1‑211 (MQ=255)                                                    aTCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTAt                                                                                                                <  1:534188/251‑1 (MQ=255)                                                      cGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATgg                                                                                                              >  1:698282/1‑251 (MQ=255)                                                      cGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATgg                                                                                                              >  1:135433/1‑251 (MQ=255)                                                        cTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTGCCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATg                                                                                                                                                                      <  1:71476/193‑1 (MQ=255)                                                        cTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGt                                                                                                             <  1:173882/250‑1 (MQ=255)                                                              atgtcccTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAgg                                                                                                       <  1:864086/246‑1 (MQ=255)                                                              atgaCCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAggt                                                                                                      <  1:828490/251‑1 (MQ=255)                                                              atgaCCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAggt                                                                                                      <  1:297241/251‑1 (MQ=255)                                                              atgaCCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAgg                                                                                                       <  1:543302/250‑1 (MQ=255)                                                                  cccTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTATGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTa                                                                                                  >  1:835662/1‑251 (MQ=255)                                                                  cccTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAgg                                                                                                       >  1:52394/1‑246 (MQ=255)                                                                            gAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCAc                                                                                         <  1:793347/250‑1 (MQ=255)                                                                              aTTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGcc                                                                                                                                     <  1:515120/204‑1 (MQ=255)                                                                                ttAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCAc                                                                                         >  1:393878/1‑246 (MQ=255)                                                                                     gTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGCATGGTTTCAGGTGGTACCACCGCACTATcgccgc                                                                                >  1:131053/1‑250 (MQ=255)                                                                                      tCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTacc                                                                                                <  1:262240/233‑1 (MQ=255)                                                                                        ttACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCAt                                                                                                                                            <  1:826955/187‑1 (MQ=255)                                                                                          aCTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAatat                                                                          >  1:604548/1‑251 (MQ=255)                                                                                                     cTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAggt                                                                                                      <  1:34320/212‑1 (MQ=255)                                                                                                       tGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTAc                                                                                                                                                                               >  1:592874/1‑137 (MQ=255)                                                                                                       tGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTaccac                                                                                              <  1:882247/218‑1 (MQ=255)                                                                                                                  tCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCtctctctg                                                  <  1:126365/251‑3 (MQ=255)                                                                                                                            taataaAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAggt                                                                                                      >  1:727822/1‑189 (MQ=255)                                                                                                                               taaAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACt                                                                                                                                                       >  1:631918/1‑137 (MQ=255)                                                                                                                                 aaGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTacc                                                                                                >  1:160177/1‑190 (MQ=255)                                                                                                                                   gTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCGTTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTacc                                                                                                >  1:49533/1‑188 (MQ=255)                                                                                                                                   gTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATcc                                                          >  1:356047/1‑226 (MQ=255)                                                                                                                                   gTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGcgctttcgcttt                                 >  1:862751/1‑251 (MQ=255)                                                                                                                                                 aTTTATTATCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTc                          >  1:861181/1‑244 (MQ=255)                                                                                                                                                    tattatCTCTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGCATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCctct                                                       >  1:190116/1‑212 (MQ=255)                                                                                                                                                         tctcTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTc                                                                        >  1:554253/1‑190 (MQ=255)                                                                                                                                                         tctcTGAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTtctc           <  1:874561/251‑1 (MQ=37)                                                                                                                                                              gAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGAt                                                                                                                                                                 <  1:296779/96‑1 (MQ=255)                                                                                                                                                              gAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTCTCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTTCTCTCGtt      >  1:437185/1‑251 (MQ=21)                                                                                                                                                              gAGCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTcc                                                   >  1:639473/1‑206 (MQ=25)                                                                                                                                                                gCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTTCTCTCGtt      >  1:658941/1‑249 (MQ=21)                                                                                                                                                                gCTATACGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTTCTCTCGtt      >  1:215043/1‑249 (MQ=21)                                                                                                                                                                   ataCGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTAt                                                                                      >  1:546963/1‑166 (MQ=21)                                                                                                                                                                   ataCGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTTCTCTCGTTCttt  >  1:310331/1‑250 (MQ=17)                                                                                                                                                                   ataCGTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTTCTCTCGTTCttt  >  1:64173/1‑250 (MQ=17)                                                                                                                                                                       gTTGTCTAGCTATAAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTc                          <  1:596511/222‑1 (MQ=14)                                                                                                                                                                                    |                                                                                                                                                                                                                                         CCCAGAGTAAATACCATTCCTAAATGCAATTTTTAAATATCCTTTTGGTCGATCGCTAGTGATGACCCTAATTCGGAATTTAATGTCTTACTTTTATGGGCTTGAATTTACAATCCAAGTTCATAATAAAGTACCTGTGTCTGAATTTATTATCTCTGAGCTATACGTTGTCTAGCTAT‑AACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACACTACCATCGGCATTACGGCCTTTCACTTCTGAGTTCGGTATGGTTTCAGGTGGTACCACCGCACTATCGCCGCCAATATTCTTTGATGACTATCCTCTCTCCTTTCCGCGCTTTCGCTTTCTTATTCCGCTCTTCGCTTCTCTCGTTCTTT  >  1/1634747‑1635157

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