Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,625,044	A→G	Q50Q (CAA→CAG)	serB →	phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,625,044	0	A	G	95.5%	196.5 / ‑0.9	67	Q50Q (CAA→CAG)	serB	phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (34/30);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

GAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑ACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCA  >  NC_000913/4624802‑4625288                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        gAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGa                                                                                                                                                                                                                                                  <  1:123965/250‑1 (MQ=255)    aagttGAGACTTAATGTTGCCAGAAGCAATTGACACAAGGTAGCCTAATGCGTTATTTTCCCTTCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTTTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCa                                                                                                                                                                                                                                              <  2:155489/248‑1 (MQ=255)                         aGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTgc                                                                                                                                                                                                                         <  1:592751/251‑1 (MQ=255)                           aaGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTgcgg                                                                                                                                                                                                                       <  1:419664/251‑1 (MQ=255)                                aTGGATACAAGGTAGCCGCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCGTAATGCCTAACATTACCTGGTGCGACCTGCCGGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATg                                                                                                                                                                                                                  <  2:483642/251‑1 (MQ=255)                                 tGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATgg                                                                                                                                                                                                                 >  1:277461/1‑251 (MQ=255)                                        aaGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGAttgt                                                                                                                                                                                                          >  1:9239/1‑251 (MQ=255)                                        aaGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTAATTAAGTGGTGATGATGTGATGCCACTGGATGACCACGCAGGTCGTGGCGGCTGGCTGCTGTATGGTCGTGGACTGGATAAACAACGTCTGGCCCAATCCAAGAGCAAACTGGGTgcggcggagggggttgt                                                                                                                                                                                                          >  2:166030/1‑240 (MQ=255)                                             agcctgcggCGTTATTTACCCGGCTTCGAACGATTTTACAGGAGCCTTAAGGCCTAACATTGCCTGGTGTGACCAGCCGGAATATGTCTCTGTATGGCCGGGTCTGCCTCTTTCAGTAAGTGGTTATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGATGCTGTATGGTCGTGGGCTGGATAAAAAACGTCTGACCCAATACCAGAGCAAACTGGGGGCAGCGATGGTGATTGTTgca                                                                                                                                                                                                      <  2:226254/242‑2 (MQ=255)                                                 tCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGGGATTGTTGCCGACTg                                                                                                                                                                                                 >  2:583119/1‑251 (MQ=255)                                                  caggCGTTATTTGCCCTGCTTCGAACGAGTTTACAGGAGCCTTAATGCCTAACATTACCTGGTTCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGTTTGTTGCCGCCt                                                                                                                                                                                                  <  1:166030/246‑1 (MQ=255)                                                   aTGCGTAATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGAGGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCTGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGt                                                                                                                                                                                               <  2:164333/251‑1 (MQ=255)                                                     gCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTg                                                                                                                                                                                              <  1:154247/250‑1 (MQ=255)                                                          tttttCCCTGCTTCGGACGATTTTACAGAAGCCTTAATGACTAACATTACATGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTACGTGGTGATGAAGTGCTGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGGGGACTGGATACAATGCGTCTGACCCACTCCCAGAGATAACGGGGTGCGGCGGTGGTGATTGTGGCCGCCGGGTGCggggg                                                                                                                                                                                        >  2:401571/2‑247 (MQ=255)                                                               cccTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGAtt                                                                                                                                                                                   >  2:130177/1‑251 (MQ=255)                                                                ccTGCTTCGAACGATTTGTCAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTa                                                                                                                                                                                  <  2:277461/251‑1 (MQ=255)                                                                ccTGCATCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCAACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGAGTGCGGCGATGGTGCTTGTTGCCGCCTGGGGCGTGGgagaat                                                                                                                                                                                   >  2:567451/1‑248 (MQ=255)                                                                     ttCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATccgg                                                                                                                                                                             >  1:267977/1‑248 (MQ=255)                                                                          aCGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGAtt                                                                                                                                                                       <  1:130177/251‑1 (MQ=255)                                                                             cttttACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGAAGCGGCTGGCTGCTGTATGGACGGGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGGGCGGCGGGGGTGATTGTTGCCGCCTGGTGCGAGGACGAGTATccagggtttcgg                                                                                                                                                                     >  2:211941/2‑240 (MQ=255)                                                                             aTTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGGGCGGCGATGGTGATTGTTGCCGCCTGGGGCGGGGAAGATTAACAGG‑GGGTTCGt                                                                                                                                                                    >  2:37570/1‑251 (MQ=255)                                                                                  aCAGGAGCCCTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTgg                                                                                                                                                                <  2:216298/250‑1 (MQ=255)                                                                                    aGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAg                                                                                                                                                             <  1:46396/251‑1 (MQ=255)                                                                                       ggCCTAAATGCCTAACATTACCGGGTTCGACCGGCCTGAAGATGTCTCTTTATGGCCGGGTTTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGtt                                                                                                                                                          <  2:276005/250‑1 (MQ=255)                                                                                            tAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTtcact                                                                                                                                                      <  1:37570/250‑1 (MQ=255)                                                                                                          aacTGGTGCGACCTGCCTGAAGAAGGATCATTATGGCCGGGTCTGCCACTTTCATTAAGTGGTGATGAAGTGAAGCCACTGGATTACAACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGGCTGACCCAATACCAGAGCAAACTGGGTGCGGCTATTGTGATTGTTGCCGCCTGGTGCGTGTAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCTCTCACCTCACGGGCTac                                                                                                                                       <  2:597717/249‑1 (MQ=255)                                                                                                          aCCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTa                                                                                                                                        <  1:450965/250‑1 (MQ=255)                                                                                                                    ccctggctgAAGAGGTCTCTTTATGGCCGGGTCTGCTTCTTTCATTAAGTGGTGATGAAGAGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGAGAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCTGCGATGGTGATTGGTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGcc                                                                                                                            <  1:650252/250‑1 (MQ=255)                                                                                                                         cctgAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCatg                                                                                                                        >  1:53565/1‑248 (MQ=255)                                                                                                                                  gTCTCTTTATGGCCGGGTTTTCCTCCTGTATTAAGTGGTGAGGAAGGGATGACACGGGATTACCACGCAGGTCGTAGCGGCGGGCTGCGGTATGTGCGTGGGCGGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTTGTACGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCTCTCCCCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCaact                                                                                                              <  2:368476/251‑4 (MQ=255)                                                                                                                                      ctTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGAt                                                                                                          >  1:201023/1‑251 (MQ=255)                                                                                                                                          aTGGCCGGGTCTTCCTCTTTCATTAAGTGGTGATGGATTGATTCCACTGGATTACCACGCAGGTCGTTGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGGCCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATggcg                                                                                                      <  2:180288/251‑4 (MQ=255)                                                                                                                                          aTGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGTGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTc                                                                                                       <  1:665292/250‑1 (MQ=255)                                                                                                                                                gggTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACAGCACGGGCTacgc                                                                                                                                     <  1:527797/215‑3 (MQ=255)                                                                                                                                                 ggTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCt                                                                                               >  1:640237/1‑251 (MQ=255)                                                                                                                                                   tCTGCCTCTTTCATTAAGTGGTGATGCAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGCGCAAACTGGGTGCGGCGATGGAGATTGTTGCCGCCTGGTGCGTGGACGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCCCCGCACGGGCTacgc                                                                                                                                     >  2:458986/1‑210 (MQ=255)                                                                                                                                                   tCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTacgc                                                                                                                                     <  1:458986/212‑3 (MQ=255)                                                                                                                                                     tGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGGCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCTCTCACCGCACGGGATCC‑GCGCCTGGCCCATGAACCGCAACTGGATGTCGCCCCGCTgggg                                                                                           >  2:116378/1‑251 (MQ=255)                                                                                                                                                      gcatctTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGAAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGACTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCCCCTGGATCGTGGAAAATT‑CAAGGAGGATTCGTCTGGCAAGTTCACGCACCGCACGGGCTACGAC‑CCTGGACAATGAAGCGCACATGGATGTCGCCCCGCGgggg                                                                                           >  2:365558/4‑250 (MQ=255)                                                                                                                                                      gCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTgggg                                                                                           >  1:541493/1‑250 (MQ=255)                                                                                                                                                                    ttggagatgaAGTGATGCCACTGGATTAACACGCAGGTCGTAGCGGCTGGCTGCTGAATGGTCGTGGGCTGTATAAACAGCGTCTGAACCAATACCAGAGCAACCTGAGTGCGGCAGTGGGAATTGTTGCCGCCTGGTGCGTGGAAGATTATCATG‑AGAGTAGTCTGGCAGGTTCACTCACCGCACGGGATAC‑GCGCCCAACCCATGAAGCGCCACTGGATGGCGTCCCGCTggggtaacacccgaccctg                                                                            >  2:395847/6‑236 (MQ=255)                                                                                                                                                                    gtggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTg                                                                           >  1:91454/1‑251 (MQ=255)                                                                                                                                                                     tggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATACCGCACCTgc                                                                          >  1:300967/1‑251 (MQ=255)                                                                                                                                                                           tgaAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGACGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑AGATTCGTCTGGCAGATTCACTCACCGCCCGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCGGCACCTGCGCaagc                                                                    >  2:572740/1‑248 (MQ=255)                                                                                                                                                                             aaGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGAGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTTACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAAATGGATGTCGCCCCGCTGGGG‑AAACACCCGCCACTGCGCACGCCg                                                                  >  2:93958/1‑251 (MQ=255)                                                                                                                                                                              agggtggccccgggagtACCACGCAGGCCGTATCGGCGTGCTGCGGTATGGTCGTTGGCTTGAGAAACAGCGTCTGACCCAATCCCAGAGCAAACTGGGTGCGGCGAGGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCCCCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCCGAGCCCGCCgg                                                                 <  2:375592/235‑1 (MQ=255)                                                                                                                                                                                 gATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGACGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGtt                                                              >  1:575090/1‑251 (MQ=255)                                                                                                                                                                                  aTGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTTACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTAGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCTCGCTGGGG‑AAAATCCAGCCCCTGCGCACGCCgggctt                                                             >  1:13450/1‑248 (MQ=255)                                                                                                                                                                                               aCCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCCGAGCAAACTGGGTGCGGCGATGGTGCTTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGGTCACTCACCGCCCGGGCTCC‑GCGCCTTGCCCATGAAGCGCGACTGGAGGTCGCCCCGCTGGGG‑CAAATCCCGCACACGCGCACCCCGGGTGTGCTGGGGATggg                                                 >  2:544003/1‑249 (MQ=255)                                                                                                                                                                                                 ccctcAGGTCGTAGCGGCTGGCGTCGGTAGGTTCGGGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAATCTGGGTGCGTCGATTGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGACTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGacat                                              <  2:516705/247‑4 (MQ=255)                                                                                                                                                                                                 cACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGacat                                              >  1:516705/1‑248 (MQ=255)                                                                                                                                                                                                 cACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCATGGGCTGGATCTACAGCGTCTGAGCCAATACCCGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGGGGACGATTATCAGG‑TGATTCGTCTGGCAGGTTCATCCACCGCCCGGGGTAC‑GCGCCAGGCCCATGAAGCGCAACAGGCTGTCGCACCGCGGGGGCAAAATCC‑GCCCCTGCGCTCGCCGGGCTTgtcggcgttggacc                                               >  2:305300/1‑237 (MQ=255)                                                                                                                                                                                                  aCGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATg                                             <  1:583119/251‑1 (MQ=255)                                                                                                                                                                                                          ggTAGCGGATGGCTGCTTGATGTGCGTGGGCTGGATAAACAGCGGCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCTCCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGTATTCCAc                                     <  2:601993/250‑1 (MQ=255)                                                                                                                                                                                                                 gctggctgCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCCAACTGGGTGCGGCGATGGCGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATc                              >  1:330072/1‑251 (MQ=255)                                                                                                                                                                                                                   tggctgCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATcc                             >  1:341839/1‑250 (MQ=255)                                                                                                                                                                                                                       tgctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAg                                                                                                                                                             >  1:121263/1‑120 (MQ=255)                                                                                                                                                                                                                       tgctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAg                                                                                                                                                             <  2:121263/120‑1 (MQ=255)                                                                                                                                                                                                                       tgctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑CAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGat                                          >  1:222322/1‑232 (MQ=255)                                                                                                                                                                                                                        gctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGagt                        <  1:394999/250‑3 (MQ=255)                                                                                                                                                                                                                             aTGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATg                   <  1:15169/250‑1 (MQ=255)                                                                                                                                                                                                                                           gATAAACCGCGTCTGACCCAATAAAAGAGCAAACTGGGTGCGGCGATGGGGATTGTTGCCGCCTGGGGCATGGAAGATTATCAGG‑TGATTCGCCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAAATGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTGTGCTGGAGATGGACATGGATTCCACCGCCATCCAGATTGAATGTCTTGATGAAATTGc    >  2:269002/1‑251 (MQ=255)                                                                                                                                                                                                                                             tAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGcc   >  1:287203/1‑250 (MQ=255)                                                                                                                                                                                                                                             tAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑GGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGGGATGGACATGGATTCCCCCGCCATCCAGATTGAATGTATTGATGAAATTGccc  >  2:141575/1‑250 (MQ=255)                                                                                                                                                                                                                                              aaaCAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATg                                             <  1:89365/207‑1 (MQ=255)                                                                                                                                                                                                                                              aaaCAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                                    >  1:70579/1‑216 (MQ=255)                                                                                                                                                                                                                                              aaaCAGCGTCTGAACCAAGACCAGAGCAAACTGGGTGCGGCGATTGTGATTTGTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑GCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                                    <  2:70579/216‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        GAAAATTTGAGACTTAATGTTGCCAGAAGCAATGGATACAAGGTAGCCTCATGCGTTATTTTCCCTGCTTCGAACGATTTTACAGGAGCCTTAATGCCTAACATTACCTGGTGCGACCTGCCTGAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGG‑TGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTAC‑ACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGG‑AAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCA  >  NC_000913/4624802‑4625288

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