Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,310,770	A→G	H32H (CAT→CAC)	yciI ←	protein YciI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,310,770	0	A	G	94.5%	177.8 / 0.9	55	H32H (CAT→CAC)	yciI	protein YciI
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (25/27);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

CCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCT  >  NC_000913/1310548‑1311013                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    ccttaaatatggaatcgtttactaGAACACTTTCTTAAATGGTTTCCCTGAAACGTGGTCATAGACTCCTGCCGCGACGTACGTGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAAtgtg                                                                                                                                                                                                                         <  2:361694/228‑1 (MQ=255)       attttGAGTCTCTAGTCAGAACTCTTTCTTAAATGGTTTCAGTTAAACGTTTTCATAGCCTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATCCAAATTCAGCAATAACGGTGGAGCCAGTAAATCACGCAGCCCCTGGATCGTTACTGTCTACCGCGGTCATTGTACCAGCTGTCAACAACCGACCTGCATCGTGCAGTATCTGTAAACGTGCTAAATGTTCCGGAc                                                                                                                                                                                                                   <  2:558068/248‑1 (MQ=255)           tGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAAc                                                                                                                                                                                                               >  1:784689/1‑251 (MQ=255)             aaTCTTTATTCAGAACACTTACTTAAAGGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACgg                                                                                                                                                                                                             <  2:330609/251‑1 (MQ=255)                   tttccataaCACGTTCTTAAAGGGTTTCACTGAAACGTGTTCATAGACGCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGaaa                                                                                                                                                                                                          <  2:833079/241‑1 (MQ=255)                     ttCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGt                                                                                                                                                                                                     >  2:28688/1‑251 (MQ=255)                          aaCACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTc                                                                                                                                                                                                >  1:565997/1‑251 (MQ=255)                             acTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAg                                                                                                                                                                                             <  1:251709/251‑1 (MQ=255)                                ttCTTAAAGGGTTCCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGc                                                                                                                                                                                           <  1:85719/250‑1 (MQ=255)                                     aaaaggTTTCACTGAAACTGGTTCTAAGCCCCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAg                                                                                                                                                                                     <  2:801294/247‑1 (MQ=255)                                      aaTGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGc                                                                                                                                                                                    >  1:471296/1‑251 (MQ=255)                                      aaTGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCACAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGGCAACAAACAACCTTCATCGTGCAGTAACTGTAACCGTGATAAATATACCGGACGAACGGAAAGGCGTTTTTCGAGGCTATccgc                                                                                                                                                                                    >  2:338462/1‑248 (MQ=255)                                         ggTTTCACTGAAACTTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCttt                                                                                                                                                                                 <  1:855644/251‑1 (MQ=255)                                           tttCACTGAAACGTGATCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCCGGCCTGACATGCTTCCAGCGCTTCAAATTCACCAATCACGGTTGCGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACAAGCTGTCAAAAAACGACCTTAATCGTGCAGTAACAGAAAACGTGATACCTGTGCCAGCCGAACTGAAAGTCGTTTTTAGAGGCTATCCGATTTGt                                                                                                                                                                               >  2:41834/1‑251 (MQ=255)                                            ttCACTGAAAAGTGTTCATAGACACCTGCCGCTACGTACGTGACAGCATCGGCCCAGGCCTGAACTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTAACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCCGTACCTGTAAACGTGCTCAATATGCCGTACGAACCGAAAGGCTTTTTTGGCAGCTATCAGCTTTggc                                                                                                                                                                              >  2:280716/1‑249 (MQ=255)                                               aCTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTg                                                                                                                                                                           >  1:680795/1‑251 (MQ=255)                                                  gAAACGGGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGa                                                                                                                                                                          <  1:705516/249‑1 (MQ=255)                                                   aaaCGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCg                                                                                                                                                                       >  1:701303/1‑251 (MQ=255)                                                       gtgtTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTaa                                                                                                                                                                    >  2:332069/1‑250 (MQ=255)                                                                   tCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATacaac                                                                                                                                                        >  1:376883/1‑250 (MQ=255)                                                                           cTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGACAGTAAATCCCGCAGACCCTGGATCGTTACTGTCTACCGCAGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGAAGAACGGAAAGGCGTTTTTCGAGGCTATCAGATATGTCTTGTGCGTAAATAACATAAAACACGGGcac                                                                                                                                                >  1:36403/1‑250 (MQ=255)                                                                                    ggTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATGGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGt                                                                                                                                      <  2:784689/251‑1 (MQ=255)                                                                                             cctcccAGGCCTGAGCTGCTTCCAGCGCTTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAg                                                                                                                             <  2:584337/248‑1 (MQ=255)                                                                                                cccAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAg                                                                                                                             >  1:584337/1‑248 (MQ=255)                                                                                                   aGGCCTGAGCTGCTTACAGCGATTCAAATTCAGCAATAACGGTTGAGCAAGTAAATCCCGCAGCCCCGGGATCGTTACTGTCTACCGCTGGCATGGGACCAGCTGTCAACACACGACCTTCATCGCGCAGAAACTGTGAACGTGCTAAATGTGCCGGCTGACCGGAAAGACGTTTTTCGAGGCTATCAGCTTTGTCTTGAACGTAAATAACATACACCCCGGGACCAAATCCTTGGTCGGGAAAGTTGTa                                                                                                                        >  2:236431/1‑250 (MQ=255)                                                                                                          aGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGt                                                                                                                >  1:551696/1‑251 (MQ=255)                                                                                                           actgataccAGCGATTCAAATACAGCAATAATGGTTGAGCTAATAAAACCCGCAGCCCCGGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGAGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGta                                                                                                               <  2:617013/244‑2 (MQ=255)                                                                                                            ctgctTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTACAGTACGtt                                                                                                               >  1:657503/1‑250 (MQ=255)                                                                                                              attttcAGCGATGCAAATTCAGCAATAACGGTTGAGTCAGTAAATCCCGAAGCCCCTGGATCGTTACTGTCTTCCGCTGGCATTGGACCAGCTGTCAACAAACGACCGTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGTGGCTATCAGCTTTGTCTTGAGCGTAAATAGCATACAACTCGGGCACAACTCCTTGTTCGGGAAAGTTGTAAGGTGCGTTAt                                                                                                             <  2:792570/245‑1 (MQ=255)                                                                                                                  ccAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAAt                                                                                                                                                                                                                            <  1:305271/135‑1 (MQ=255)                                                                                                                  ccAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAAt                                                                                                                                                                                                                            >  2:305271/1‑135 (MQ=255)                                                                                                                      cGATTCAAATTCATCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAATGCGTTTTGCGAGGCTATCATCTTTGTCTTGAGCGTAAATAACATACAACACGGGCAGAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAggg                                                                                                    <  1:36534/251‑1 (MQ=255)                                                                                                                      cGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAggg                                                                                                    >  1:822488/1‑251 (MQ=255)                                                                                                                        aTTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCgg                                                                                                                                                                                                                     <  2:607977/136‑1 (MQ=255)                                                                                                                        aTTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCgg                                                                                                                                                                                                                     >  1:607977/1‑136 (MQ=255)                                                                                                                           cAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTAGATCGTTACTGTCGACCGCTGGCATTGGACCAGCTGTCAACAAACGTCCTTCATCGTGCATTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTGTTTCGAGGCTAGCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAg                                                                                               <  1:549136/251‑1 (MQ=255)                                                                                                                             aagacagtaaTAACGGTTGAGCCAGTAAACCCCGCAGCCCCTGGATCGTTACTGGCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGGTGTAAAGTACGTTATTTGAAAGGGCGAAGAt                                                                                             <  2:657503/243‑1 (MQ=255)                                                                                                                             aaTTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACACACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCgc                                                                                                 >  1:318864/1‑246 (MQ=255)                                                                                                                                     aaTAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACg                                                                                     <  1:23773/251‑1 (MQ=255)                                                                                                                                          cGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCAGCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAg                                                                                <  2:44810/251‑1 (MQ=255)                                                                                                                                            gTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAgat                                                                              <  2:701303/251‑1 (MQ=255)                                                                                                                                                     gTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTTAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCttt                                                                     >  1:302123/1‑251 (MQ=255)                                                                                                                                                               cAGCCCCTGGATAGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAATCGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTCTTCGCGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTGGTTCGGGTCAGTTAAAAAGTACGTTATTTGAAAGGGCGATGATCTGCGACGGAAAGGTGATGTGTTTGTTAAgacc                                                            >  1:798491/1‑247 (MQ=255)                                                                                                                                                                 gCCCCTGGATCGTTACTGGCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCAGCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACATCACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATg                                                         <  2:741847/251‑1 (MQ=255)                                                                                                                                                                  ccccTGGATCGTTACTGTCTCCCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGc                                                        <  2:279344/251‑1 (MQ=255)                                                                                                                                                                       tgATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGAACGTTATTTGAAAGGGCGAAGTTCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATaa                                                    >  2:355785/2‑250 (MQ=255)                                                                                                                                                                           cGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCTAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTa                                                >  1:482837/1‑250 (MQ=255)                                                                                                                                                                                  gTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTa                                          <  1:28688/249‑1 (MQ=255)                                                                                                                                                                                  gTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGGTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTa                                          <  2:388379/249‑1 (MQ=255)                                                                                                                                                                                         gCTGGCCTTGGACCAGCTGTCAACAAACGACCTTCATCGTGCCGTAACTGTACACGTGCTAAATGTGCCGGACGAACGTAAAGGCGTTTTTCGGGGCTATCCGCTTTGTCTTGAGCGTAAATAACATACGACACGGGCACAACTCCTTGTTCGGGAAAGTTGTCAAGTACGTTATTTGAGAGGGCGAAGATCTGCAACGGAAAACTGATGTATTTGTTAAGGCCATGCATAAAGTAAGAGAAATTACGcc                                  >  2:774197/1‑250 (MQ=255)                                                                                                                                                                                                                   aCGACCTTCATCGGGCAGTAACTGTAAACGGGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCActgatcctgat       <  1:401921/251‑1 (MQ=255)                                                                                                                                                                                                                   aCGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGTGCGAAGATCTGCAACGGAAAGATGATGTCTCTGTTAAGGCCATGCATAAAATAAGGGTAATTACGCCAAAAATGACATTTTCActgatcc           <  1:22097/247‑1 (MQ=255)                                                                                                                                                                                                                   aCGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGTGCGAAGATCTGCAACGGAAAGATGATGTCTCTGTTAAGGCCATGCATAAAATAAGGGTAATTACGCAAAAAATGACATTTTCActgatcc           >  2:22097/1‑247 (MQ=255)                                                                                                                                                                                                                   aCGAACTTCATCGTGCAGTAACTGTAAACGTTCTAAATGTGCCGGAAGCACGGAAAGGCGTTTTGCGAGGCTATCAGCTTTGTCTTTAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTGGTAAAGTACGTTATTTGAAAGGGCGAATATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGGAAGGGTAATTACGCCAAAAATGACAATTTCActgatcctgat       <  2:80511/251‑1 (MQ=255)                                                                                                                                                                                                                    cGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGTTAATTACGCCAAAAATGACATTTTCActgatcctgatc      <  1:331630/251‑1 (MQ=255)                                                                                                                                                                                                                          tcatcgTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCt  <  1:332069/249‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    CCTTAAATATTGAATCTCTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCCTGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCT  >  NC_000913/1310548‑1311013

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