Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	262,762	Δ1 bp	intergenic (+6/‑109)	proA → / → thrW	glutamate‑5‑semialdehyde dehydrogenase/tRNA‑Thr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	262,762	0	A	.	96.9%	265.5 / ‑1.1	65	intergenic (+6/‑109)	proA/thrW	glutamate‑5‑semialdehyde dehydrogenase/tRNA‑Thr
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (16/47);  total (18/47)
Fisher's exact test for biased strand distribution p-value = 7.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCTCCCACCAAAAAATTATCTTAATGTAACAGC  >  NC_000913/262525‑263012                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           cGATGCGATCCTGACCCGCGATATTCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGTCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCaaa                                                                                                                                                                                                                                              <  1:414846/251‑1 (MQ=255)     gCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCCCG‑CGC‑GTGGCCCAATGGGGCTGGAAGCCCTGGCCCACT‑ACAAGTGGG‑TCGGCATTGGTGGTTAACCCATTCGGGCGTAAATAAAAC‑GGGTGGTGCCAAAGt                                                                                                                                                                                                                                           >  1:226806/1‑250 (MQ=255)       gATCCTTACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTTGCGGTAAGCACAC‑A‑AGAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGc                                                                                                                                                                                                                                        <  1:40592/251‑1 (MQ=255)              aCCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCCCG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TAGGCATTGGTGATTACACCATTCGTGCGTAAATAAAAC‑GGGGGGTGCAAAAGTAGCCAttttg                                                                                                                                                                                                                                 >  1:119232/1‑249 (MQ=255)              aCCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGa                                                                                                                                                                                                                                 <  1:2400/251‑1 (MQ=255)                       atGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACGCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑AGCACTTACAATTGGA‑TCGGCTTTGGTGATTACAGCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAg                                                                                                                                                                                                                         <  1:201708/250‑1 (MQ=255)                            cAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTATCCATTTGATTCACAAGGCCAtt                                                                                                                                                                                                                   <  2:70554/251‑1 (MQ=255)                            cAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCAtt                                                                                                                                                                                                                   <  1:295426/251‑1 (MQ=255)                               cGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGAc                                                                                                                                                                                                                <  1:365133/251‑1 (MQ=255)                               cGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAAATCCACG‑CGC‑GTGGCCCAAAGGGGGTGGGAGCACCG‑ACCACTTCCAAATGGAATCGGC‑TTGGTGGTTAAACCATTCGTGCGGAAATAAA‑CCCGGTGGTGCAAAAATAACCACTTTATTCACCAGGGCGTTGAc                                                                                                                                                                                                                >  2:305637/1‑251 (MQ=255)                                   cAGCGTTTTGTTAACGAAGTGGCTTCGTCCGCTGTTTAAGTTAACGCCTATACGCGTTTTAACGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGACCAC‑A‑AAAACACCACG‑‑GCAGGGGGCCAATGGGGGAGGCAGGACCG‑GCCACCTACAAATGGG‑TCGGCATCGGGGGTTAAACCATTCTTGCGTAAATAAA‑CACGGTGAAGCAAACATAGCCAttttgttctacaagcccttatacgcc                                                                                                                                                                                                            >  2:183034/1‑228 (MQ=255)                                    aGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATc                                                                                                                                                                                                           <  2:10702/251‑1 (MQ=255)                                         ttttttAACGAAGTGGATTCGTACGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGccc                                                                                                                                                                                                       <  2:97840/246‑1 (MQ=255)                                         tttGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGccc                                                                                                                                                                                                       <  1:305637/250‑1 (MQ=255)                                         tttGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGGGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGTAGTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGccc                                                                                                                                                                                                       <  1:174131/250‑1 (MQ=255)                                          ttGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGGTGGAAGCACTG‑ACCACTTACAACGGGA‑GAGGCATTGGGGGTTTTACCATTCGCGCGTAAATAAA‑CCGGGGGATGAAAAAAAAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCgg                                                                                                                                                                                                     >  1:229674/1‑251 (MQ=255)                                                 cGAAGTGGATTCGGCCGCTGTTTACGGTAACGCCTCTACGCGGTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCAGTTTATTCACAAGGCCATTGACGCATCGCCCGGTTAGttt                                                                                                                                                                                              <  2:366546/251‑1 (MQ=255)                                                  gAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAACACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGG‑TCGGCATTGGGGATTACACCATTCGGGCGTAAATAAA‑CCGGGGGATGCAAAAGTAGCCATTTTATTCACAAGGCCATTGACGCATCGCCCGGTTAGtttt                                                                                                                                                                                             >  1:95837/1‑251 (MQ=255)                                                         aTTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACCC‑A‑AAAACCCCCCG‑CGC‑GGGGCCCAATGGGGCTGGAAGCACTG‑AGCACTTACAAGAGGA‑TCGGCATTGGTGATTTCACCATTCGGGCGTAAATAAA‑CCGGGGGAGGCAAAAGTAGCCATTTGATTCACAAGGGCATTGACGCATCGCCCGGGTAGTTTTTACCTTg                                                                                                                                                                                      >  1:187996/1‑251 (MQ=255)                                                                 gCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGt                                                                                                                                                                              <  2:226806/251‑1 (MQ=255)                                                                 gCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGCCA‑CACCCCA‑AAAAATCCACC‑AGC‑GTGGGCCACTGGGGCTGGGAGAACTG‑ACCCCTTACAACTGGG‑TGGGGATTGGGGGTTACACCCTTCTTGCGTCAATAAA‑CCGGGTGATGCCAAAGAAGCCATTTTATTCACAAGGACATTGACGCGTCCCCCGGGTTGtttttaacttttcccccgc                                                                                                                                                                              >  2:170078/1‑236 (MQ=255)                                                                   tGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACCC‑A‑AAAACACCCCG‑CGC‑GTGGCCCCATGGGGCTGGAAGCACTG‑ACCACTTACAAGGGGA‑TCGGCATTGGTGTTTACACCCTTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGACATTTGATTCACAAGGCCATTGACGCATCGCCCGGGTAGTTTTAACCTTGTCCACCGTgt                                                                                                                                                                            >  1:77719/1‑250 (MQ=255)                                                                   tGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAACAC‑C‑ACAAAACTCCACGTCG‑‑GTGGCCCCCAGGGGCTGGGAGCACTG‑ACCACTTACAAGGGGA‑TCGGCATTGGGGATTACACCATTCGTGCGTAAATAAA‑CCGGGGGATGCAAAAGTAGCCCCTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCtttttccccgtg                                                                                                                                                                             >  1:231381/1‑240 (MQ=255)                                                                        aCGTTAACGCCTCTCCGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAACACAC‑A‑AAAACTCCCCG‑CGC‑CTGGCCCCCTGGGGCTGGGAGCACCG‑ACCCCTTTCCAATGGG‑TCGGCATTTGTGATTTCACCCTTCCTTCCTCCATATA‑CCGGGTGATTCACAACAAACCCTTTGATTCACCAAGCCCTTGACGCCTCCCCCGGTTCGTTTTTACCTTTTTCACCGTGCTTccc                                                                                                                                                                       >  1:351174/1‑249 (MQ=255)                                                                            tAACGCCTCTACGCGTTTGACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑GCGGCATTGGTGATTACACCATTGGTGCGTAAATAAA‑CCGGGTTATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGTTCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGtt                                                                                                                                                                    <  1:34967/250‑1 (MQ=255)                                                                             aaCGCCGCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTTGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACGCCACG‑CGC‑GTGGCCCAATGGCGCTGGAAGCACTG‑ACCACTTACAAGTGTA‑TCGGCATTGGTGATTACACCATTCGTGCGGAAATAAA‑CCGGTGGGTGCAAAAGGAGCCATTTGATTCACAAGGCCATTGACGCATCCCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTc                                                                                                                                                                   <  1:417737/250‑1 (MQ=255)                                                                                      aCGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACa                                                                                                                                                            <  1:220695/248‑1 (MQ=255)                                                                                        gcgTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTcc                                                                                                                                                       <  2:284150/251‑1 (MQ=255)                                                                                         cgTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCt                                                                                                                                                      <  1:14422/251‑1 (MQ=255)                                                                                          gtgttACCGACGGCTGCCAGTTTGGTCGGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACGCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑GCGGCATTGGTGATTACCCCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTATCCATTTGATTCACAAGGCCATTGACGCAGCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCtt                                                                                                                                                     <  2:374538/248‑1 (MQ=255)                                                                                             ttACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAATCTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATGCACAAGGCCATTGACTCATCGCCCGGTTATTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACAGGTCCTTTc                                                                                                                                                   <  2:372333/250‑1 (MQ=255)                                                                                                 cGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGc                                                                                                                                              <  1:346316/251‑1 (MQ=255)                                                                                                  gccggcggCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTTATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTGCAGGGcc                                                                                                                                             <  1:123386/249‑1 (MQ=255)                                                                                                    cggcggCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCa                                                                                                                                                                        <  1:318631/222‑1 (MQ=255)                                                                                                            aGTTTGGTCTGGGTGCGGAAGTTGCGGTTAGCACAC‑A‑AAAACTCCACT‑CGC‑GTGGCCCAAGGGGGCTGGATGCACCG‑ACCACTTACAAGTGGA‑TCGGCATGGGTAATTACACCATTCGTTCGTAAATAAA‑CCTGGTGATGCAAAAGTATCCATTTGAGTCACAAATCCAGTTACGCATCTCCCGTTTATTTTTAACCTTGTCCTCCGGGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTc                                                                                                                                   <  1:379967/251‑1 (MQ=255)                                                                                                             gTTTGGCCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCa                                                                                                                                  <  1:199059/251‑1 (MQ=255)                                                                                                                           cGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGagcagc                                                                                                                    <  1:48656/251‑1 (MQ=255)                                                                                                                                  ggCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCg                                                                                                             >  2:300036/1‑251 (MQ=255)                                                                                                                                  ggCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATGCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTGTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCg                                                                                                             <  1:183034/251‑1 (MQ=255)                                                                                                                                      gTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTTATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTTGTAGAGCAGCGCATTCGTAAt                                                                                                         <  1:300036/251‑1 (MQ=255)                                                                                                                                      gTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCAGTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTGTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGAGATAGCTCAGTTGGTAGAGCAGCGCATTCGTAAt                                                                                                         <  2:318941/251‑1 (MQ=255)                                                                                                                                              cc‑a‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGt                                                                                                 <  1:26686/250‑1 (MQ=255)                                                                                                                                                 a‑aaaaTTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGAGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTACACCGTGATTCACGTTCGTGACCCTGTCCTTTCAGGGCCGATATCTCTCCGTTGGTAGAGCAGCGCATTCGTAATTCGAAGGTcc                                                                                               >  2:375280/1‑250 (MQ=255)                                                                                                                                                   aaaaCTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCAcc                                                                                                                                                                                <  1:179845/169‑1 (MQ=255)                                                                                                                                                   aaaaCTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCAcc                                                                                                                                                                                >  2:179845/1‑169 (MQ=255)                                                                                                                                                      aCTCCAAG‑CGA‑GTGGCCCCATGGGGCTGGACGCACTG‑ACCACATACAAGTGGC‑TCGGCATTGGTGATTACACCCTTCGTGAGTAAATAAA‑CCGGGTGATGAAAAAGAAGCCATTTGATTCAAAAGGCCCTTGACGCATCGCCCGGTTAGTTTTACCCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAgg                                                                                           >  2:396782/1‑251 (MQ=255)                                                                                                                                                          cACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCTTAGGTTCg                                                                                       <  1:350269/251‑1 (MQ=255)                                                                                                                                                               cgc‑gTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGGAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTATAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTc                                                                                   <  1:253654/251‑1 (MQ=255)                                                                                                                                                               cgc‑gTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGGAAATAAA‑CCGGGTGATGCAAAAGTGGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTATCCTTGTCCACCGTGATTCACGTTGGTGAACATGTGCGTTCAGTGCCGATAGAGCTCAGGTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTc                                                                                   <  1:226728/251‑1 (MQ=255)                                                                                                                                                                         cAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACGCCTATTATCg                                                                          <  2:181373/251‑1 (MQ=255)                                                                                                                                                                                   ggAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCAtt                                                                  <  1:111767/249‑1 (MQ=255)                                                                                                                                                                                    gAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATGAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTaaa                                                               <  2:194513/251‑1 (MQ=255)                                                                                                                                                                                     aaGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTaaaa                                                              <  2:6381/251‑1 (MQ=255)                                                                                                                                                                                         aCTG‑ACCACTTACAAGGGGA‑TCGGCATTGGTGATTACACCATTCGTGCGGAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGGCCTTGCAGGGCCGATATAGCTCAGGTGGGAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCaa                                                          <  2:181564/251‑1 (MQ=39)                                                                                                                                                                                               ccACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGt                                                     >  2:131997/1‑251 (MQ=255)                                                                                                                                                                                                cacctACAAGTGGA‑TCGGCATTGGTGATTACCCCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCAGTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCGATTATCGGCACCATTAAAATCAAATTGtt                                                    <  2:429737/247‑1 (MQ=255)                                                                                                                                                                                                        aGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTa                                               >  1:324819/1‑248 (MQ=255)                                                                                                                                                                                                         gTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAACTAAA‑CAGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCACATTGTTACGTa                                               >  1:383565/1‑247 (MQ=255)                                                                                                                                                                                                                       ggTGATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCTCCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCTcc                              <  1:33368/251‑1 (MQ=255)                                                                                                                                                                                                                          gATTACACCATTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCTCccac                           <  1:116665/251‑1 (MQ=255)                                                                                                                                                                                                                              acacCAGTCGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGGCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAAGTGTTACGTAAGATCTTATCATTCTCCCACCaaa                       <  2:294823/251‑1 (MQ=255)                                                                                                                                                                                                                                cccATTCGTGCGTAAATAAA‑CCTGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGAGTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTTTAGCAGCGCATTCGTAATGCGAAGGTCTTAGGTTCTACTCCTATTATCGGCACCATTAAAGTCAAATTGTTACGTAAGATCTTATCATTCTCCCACCaaaaa                     <  1:411643/250‑1 (MQ=255)                                                                                                                                                                                                                                      cGTGCGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATGCTCCCACCAAAAAATTATc               <  2:58474/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCGTAAATAAA‑CCGGGGGATGCAAAAGTAGCCATTTTATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCTCCCACCAAAAAATTATCTTa            <  2:241036/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGTAAATAAA‑CCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCTCCCCCCAAAAAATTATCTTaa           <  1:152155/251‑1 (MQ=255)                                                                                                                                                                                                                                                  aaaCCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCTCCCACCAAAAAATTATCTTAATGTAACa    <  1:170078/251‑1 (MQ=255)                                                                                                                                                                                                                                                    aCCGGGTGATGCAAAAGTAGCCATTTGATGCACAAGGCCATTGACGCATCGCCCGGTGAGTTTTAACCGTGTCCACCGTGAGTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCGCCCACCAAAAAATTATCTTAATGTAACAGc  <  2:324819/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           CGATGCGATCCTGACCCGCGATATGCGCAACGCCCAGCGTTTTGTTAACGAAGTGGATTCGTCCGCTGTTTACGTTAACGCCTCTACGCGTTTTACCGACGGCGGCCAGTTTGGTCTGGGTGCGGAAGTGGCGGTAAGCACAC‑A‑AAAACTCCACG‑CGC‑GTGGCCCAATGGGGCTGGAAGCACTG‑ACCACTTACAAGTGGA‑TCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACGCATCGCCCGGTTAGTTTTAACCTTGTCCACCGTGATTCACGTTCGTGAACATGTCCTTTCAGGGCCGATATAGCTCAGTTGGTAGAGCAGCGCATTCGTAATGCGAAGGTCGTAGGTTCGACTCCTATTATCGGCACCATTAAAATCAAATTGTTACGTAAGATCTTATCATTCTCCCACCAAAAAATTATCTTAATGTAACAGC  >  NC_000913/262525‑263012

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