Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,259,283	G→A	H123Y (CAT→TAT)	tdcG ←	L‑serine deaminase III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,259,283	0	G	A	95.0%	195.4 / NA	60	H123Y (CAT→TAT)	tdcG	L‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/31);  minor base C (2/1);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 5.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCG  >  NC_000913/3259036‑3259528                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      aTTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCg                                                                                                                                                                                                                                                    <  1:340827/251‑1 (MQ=255)     ccGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATATCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGc                                                                                                                                                                                                                                                >  1:588656/1‑251 (MQ=255)       gCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGCGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCa                                                                                                                                                                                                                                               >  1:232002/1‑250 (MQ=255)       gCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTAAGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATCGCGAGGCa                                                                                                                                                                                                                                               >  2:446246/1‑250 (MQ=255)            gctgcGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGGTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACAtt                                                                                                                                                                                                                                          <  1:65578/250‑1 (MQ=255)               gctcAGCGCTAGCTCGTGGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCCCCTGCTGAGTGGAAATCGTAATGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCg                                                                                                                                                                                                                                       <  2:360076/247‑1 (MQ=255)               gcgcAGCGCTATCTCGTTGTGCATAATCAGCCCATATATACACAGGCCGTTGTACTCACACATTTTCTGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTCCTGCCGTTTCGTCATCGTTCGCCAGGCCGACGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTTTCCCTTCCAGGCAGTGATCCTCATACCGTGCTCATCGCGAgtctccatttc                                                                                                                                                                                                                                        <  2:283197/249‑11 (MQ=255)                      gcTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCAGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGGTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGaa                                                                                                                                                                                                                               <  1:482678/251‑1 (MQ=255)                      gcTAGCGCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCGCCTGCTGAGTGGAAATCGTAAGGTCCGGACGTTTCGACATCTTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGGCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACAGATCGGGATGGaa                                                                                                                                                                                                                               <  2:511338/251‑1 (MQ=255)                         aGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGa                                                                                                                                                                                                                             <  2:212870/250‑1 (MQ=255)                         aGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGGTCTCATAGCGAGGCAACATTTCGGGATTGAAGAt                                                                                                                                                                                                                            <  1:456162/251‑1 (MQ=255)                                   gCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACCCATTTTCAGCAGTTCACCTGGTGAGTGGAAATCGGAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTCCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACGTAATAGCTCTTCCTTTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCtt                                                                                                                                                                                                                   <  2:421393/250‑1 (MQ=255)                                     atcatcAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGc                                                                                                                                                                                                                >  2:140807/1‑251 (MQ=255)                                       catcaGTCAAGATATAGCCAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGGGGAAAGCGTAAGGTACGGACGTTTCGAAATCTTGCGTAAGGCCTAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCTGTGATCCTCATTCCGTTCTCATAGCGAGGCAACATTTCTGGATGGAAGATAATGTTCTTTTCTa                                                                                                                                                                                                              <  2:50324/251‑1 (MQ=255)                                       catcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCCCCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGAGCCGCATTCCTTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCGTTGCTa                                                                                                                                                                                                              <  2:253525/251‑1 (MQ=255)                                                    ctataCAGGCCGTTATAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAACTCGTCAGGTACGGACGTTTCGACATCGTGCGACAGGCCGGAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGGCCCTTCCAGGCAGTGATCCGCATCCCGTTCTCAGAGCGAGGCAACATTTCGGGGTGGAAGATAATGTTCTTTACTACAGGAAAATcaac                                                                                                                                                                                                 >  2:235372/5‑251 (MQ=255)                                                         cAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAatat                                                                                                                                                                                            >  1:232727/1‑251 (MQ=255)                                                           ggCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATgc                                                                                                                                                                                          <  1:457141/251‑1 (MQ=255)                                                           ggCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGGAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTGCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATgc                                                                                                                                                                                          <  1:577510/251‑1 (MQ=255)                                                                          acaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACt                                                                                                                                                                           <  1:459145/251‑1 (MQ=255)                                                                                 cagcagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAgcc                                                                                                                                                                    <  1:235352/251‑1 (MQ=255)                                                                                        tCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCt                                                                                                                                                             >  2:474360/1‑251 (MQ=255)                                                                                         cACCTGCTGAGTGTAAATCGTAAGGTCCGGCCGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCATATGCCACTGGCAGCCGCCCGCTg                                                                                                                                                            <  2:232002/251‑1 (MQ=255)                                                                                          aCCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCGCTTCCTGTCCCTTCCAGGCAGTGATCCGCAGTCCGTTCTCATAGCGAGGCAACAGGTCGGGATGGAAGATAATGTTCGTTGCGACAGGAAAATCAACAATATGCGCACCAGATGCCACTAGCAGCCGCCtgctgc                                                                                                                                                           <  2:492620/251‑1 (MQ=255)                                                                                           ccttctgAGTGGAAATCGTAAGGTACGGATGTTTCGACACCGTGCGACAGGCCGAAGTGTTCTTCTCCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCTCTGGCAGCCGCCCGCTgcg                                                                                                                                                          <  2:121257/247‑1 (MQ=255)                                                                                             tgctgAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGcgcg                                                                                                                                                        >  1:185972/1‑251 (MQ=255)                                                                                                tgAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCa                                                                                                                                                                             <  2:384074/227‑1 (MQ=255)                                                                                                tgAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCa                                                                                                                                                                             >  1:384074/1‑227 (MQ=255)                                                                                                     ggAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACACGCCGCAGACAGAGTAATAGGGTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGCGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGCTGCCACTGGCAGCCGCCCGCTGCGCGTTACTAAc                                                                                                                                                >  1:358595/1‑251 (MQ=255)                                                                                                            gTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCtata                                                                                                                                          >  1:520120/1‑250 (MQ=255)                                                                                                            gTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGAGCGTTACTAACTCtata                                                                                                                                          >  2:227423/1‑250 (MQ=255)                                                                                                                  tACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGc                                                                                                                                     >  2:114359/1‑249 (MQ=255)                                                                                                                       aCGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGAtc                                                                                                                              >  1:360336/1‑251 (MQ=255)                                                                                                                        cGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACGAACTCTATAAATGCAGGGAtct                                                                                                                             <  1:140807/251‑1 (MQ=255)                                                                                                                           ttCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCtcat                                                                                                                          >  1:286907/1‑251 (MQ=255)                                                                                                                             cGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCCTAGCGAGGCAACATTTCGGGATGGAAGCTAATGTACTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCtcatc                                                                                                                         >  1:567112/1‑250 (MQ=255)                                                                                                                                        gACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCAGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGATACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATc                                                                                                             >  2:18639/1‑251 (MQ=255)                                                                                                                                              ccGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCgg                                                                                                       <  1:329808/251‑1 (MQ=255)                                                                                                                                                gAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCGCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGAc                                                                                                     <  1:492317/251‑1 (MQ=255)                                                                                                                                                             tcGACAATAAACCCGCCGCCGACAGAGTAATAGGTTGTACTTAATAGCTCTTCCTGTCCCTTACAGGCAGTGATCCGTATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAAGGTTCTTTGTTACAGGAAAATCAACAATATGCGCACCAGATGCCACGGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGc                                                                                           <  1:261289/248‑1 (MQ=255)                                                                                                                                                               gACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTGCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTc                                                                                      <  1:284016/251‑1 (MQ=255)                                                                                                                                                                              gccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTAc                                                                                                                                                    <  1:172307/174‑1 (MQ=255)                                                                                                                                                                              gccgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTAc                                                                                                                                                    >  2:172307/1‑174 (MQ=255)                                                                                                                                                                                cgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGAc                                                                      >  2:529308/1‑250 (MQ=255)                                                                                                                                                                                cgACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGAc                                                                      >  2:15518/1‑250 (MQ=255)                                                                                                                                                                                 gACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACAt                                                                    >  1:513668/1‑251 (MQ=255)                                                                                                                                                                                     gagTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGt                                                                <  1:140683/251‑1 (MQ=255)                                                                                                                                                                                        tAATAGGATTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGGCCCATGATGATGGCGACATCCGTGGc                                                             >  1:510813/1‑251 (MQ=255)                                                                                                                                                                                                     ttAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATCTGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCTGGGCTTGGCCTTTGccc                                                >  2:489584/1‑251 (MQ=255)                                                                                                                                                                                                          aGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGCTGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTc                                             >  1:302477/1‑249 (MQ=255)                                                                                                                                                                                                              ctTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCGTTGCCCGTCAGTg                                         <  1:186153/249‑1 (MQ=255)                                                                                                                                                                                                                      cccTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGAcc                                <  1:127177/250‑1 (MQ=255)                                                                                                                                                                                                                       ccTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCg                              >  1:162431/1‑251 (MQ=255)                                                                                                                                                                                                                        cTTCCAGGCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGt                             >  2:434459/1‑251 (MQ=255)                                                                                                                                                                                                                              ggcagtattCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCGTTGCCCGTCAGTGACAACGACCCGTACAGAt                       <  2:232727/243‑1 (MQ=255)                                                                                                                                                                                                                              ggCAGTGATCCGCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGAt                       <  1:499549/251‑1 (MQ=255)                                                                                                                                                                                                                                         gCATTCCGTTCTCATCGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGAtct                                                                                                                             >  1:454190/1‑138 (MQ=255)                                                                                                                                                                                                                                         gCATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGAtct                                                                                                                             <  2:454190/138‑1 (MQ=255)                                                                                                                                                                                                                                          cATTCCGTTCTCATAGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAatat           <  1:488160/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ctcATAGCGAGCCACCATTACGGGATGGAAGATAATGTTCTTTGCTAAAGGAAAATCAACAATATGCGCACCAGATGCAACTGGCAGCCGCACGCTGCGCGCTACTAACACTATAAATGCAGGGATATCATCAATGACAACCTCCTGCGGACTGTTTCATGCCCGCCACATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCCGTGACAACGACCCGTACAGATCGACCCCAATATGCCtccccc  >  2:525786/1‑247 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGCATTCCGTTCTCATGGCGAGGCAACATTTCGGGATGGAAGATAATGTTCTTTGCTACAGGAAAATCAACAATATGCGCACCAGATGCCACTGGCAGCCGCCCGCTGCGCGTTACTAACTCTATAAATGCAGGGATCTCATCAATGACAACATCCTGCGGACTGTTTCCTGCCAGTCCCATGATGATGGCGACATCCGTGGCATGGCCTTTGCCCGTCAGTGACAACGACCCGTACAGATCGACCACAATATGGCTCGTCG  >  NC_000913/3259036‑3259528

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