Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,445,924	C→T	Q605* (CAG→TAG)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,445,924	0	C	T	98.4%	198.6 / NA	64	Q605* (CAG→TAG)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (27/36);  minor base G (0/1);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGACCCGCG  >  NC_000913/4445680‑4446170                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        gCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGAGCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCg                                                                                                                                                                                                                                                   <  2:31668/251‑1 (MQ=255)       cAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGa                                                                                                                                                                                                                                             <  1:125918/251‑1 (MQ=255)        aGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGa                                                                                                                                                                                                                                             >  1:185014/1‑250 (MQ=255)        aGTGCCGAAGGTAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTGGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGaa                                                                                                                                                                                                                                            <  2:70521/251‑1 (MQ=255)               aaGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGAc                                                                                                                                                                                                                                     <  1:109107/251‑1 (MQ=255)               aaGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGAc                                                                                                                                                                                                                                     <  1:245955/251‑1 (MQ=255)               aaGTGAAAGGAGAGCTGGGGGGAAAAGATCTCAATCGTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCGTGGCGGGACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGGGGAAGCGCCACAACTACGGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCTTGGCGCAGGTTCGCGGGGAAGGCGATATCAAATCCGCCGATTAGATCGCCGGGAAACTCGAc                                                                                                                                                                                                                                     <  2:420031/251‑1 (MQ=255)                   gAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTa                                                                                                                                                                                                                                  <  1:350095/250‑1 (MQ=255)                       ggCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGt                                                                                                                                                                                                                             <  1:125566/251‑1 (MQ=255)                        cgatcgtGGGGGTAAAAGATCTCAAGCGTGATGCCACCATCAACGCGCCGGGGCTGGAGAACGCGCGGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGGACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGtt                                                                                                                                                                                                                            <  2:221628/245‑1 (MQ=255)                         cGAGCTGGGGGTAAAAGACCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGAGATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTg                                                                                                                                                                                                                           <  2:216344/251‑1 (MQ=255)                              tGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGtt                                                                                                                                                                                                                            >  1:221628/1‑245 (MQ=255)                                ggggTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAAt                                                                                                                                                                                                                    <  1:395734/251‑1 (MQ=255)                                   gTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTc                                                                                                                                                                                                                 >  1:132643/1‑251 (MQ=255)                                        aGATCTCAATTTTGATGCCACCACCAACGCGCCGGGGCTGGAGAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCGTGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCGGCGCTGGCAGGAACCTTCCGTGGCGCAGGTTCGCGTGGAAGTCGATAGCAAATCCACCGATTAGATCGCCGGGAAACCCGACGTACGCGTTGAGCAAATGTCGCAGc                                                                                                                                                                                                            <  2:474526/251‑1 (MQ=255)                                           gctcaatcgggacgcccccaTCAACGCGCCGGGGGTGGATAACGCGATGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGAACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATAGCACCGCGCGCGGCCTGCGCTGGCAGGAACTTGCCGTGGCGCAGGGTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCgg                                                                                                                                                                                                         <  2:246192/234‑1 (MQ=255)                                                    tGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTaaa                                                                                                                                                                                                  <  1:20892/249‑1 (MQ=255)                                                           accaTCAACGCGCCGGGGCGGGATAACGCGCGGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCt                                                                                                                                                                                         <  2:185014/251‑1 (MQ=255)                                                                 aaCGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAc                                                                                                                                                                                   >  1:92464/1‑251 (MQ=255)                                                                   cgcgCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGc                                                                                                                                                                                                            >  1:25045/1‑224 (MQ=255)                                                                            gCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCa                                                                                                                                                                       <  1:407573/251‑1 (MQ=255)                                                                                  tcacgagcTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGGTAAAGGACGCGGCACGGTGGAAGCGCCACAACTACTAGCAGATATCACCGCGCGCCGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCCCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCCGCCGGCTGGAAATATCAACCTCG‑GACACCTGAATGCCAAAGGc                                                                                                                                                                  >  2:222994/7‑250 (MQ=255)                                                                                   aaCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGc                                                                                                                                                                  >  1:175338/1‑249 (MQ=255)                                                                                   aaCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGc                                                                                                                                                                  >  1:198540/1‑249 (MQ=255)                                                                                     cgcgCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAAACTCGTCAC‑CCTGAATGCCAAAGGCAGcc                                                                                                                                                              >  2:459399/1‑250 (MQ=255)                                                                                           gCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGaaacgc                                                                                                                                                        >  2:66516/1‑248 (MQ=255)                                                                                             cGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAagc                                                                                                                                                        >  1:427256/1‑249 (MQ=255)                                                                                              ggggCTTGGCGTTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATAGCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGGACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCACGGCAC‑CCTGAATGCCAAAGGCAGCGAAAagcagc                                                                                                                                                     <  2:235893/251‑1 (MQ=255)                                                                                               gggctgaggggaaCAGCGAAAGGGTGGGGGAAAGTACGCGGCACGGGGGAAGCGCCGCAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTGCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACGCGACGTACGCGTTGAGCAAATGTCGCAGCCGGATGTAAATATCAACCGCGTCAC‑CCTGAATTCCAAAGGCAGCGAAAagcagct                                                                                                                                                    <  2:311910/239‑2 (MQ=255)                                                                                                ggCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAAGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCAc                                                                                                                                                   >  1:440664/1‑251 (MQ=255)                                                                                                 gCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACg                                                                                                                                                  >  2:29110/1‑251 (MQ=255)                                                                                                       ggggaCAGCGAAAGGGCTGGGGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTAc                                                                                                                                            <  2:414345/247‑1 (MQ=255)                                                                                                         gTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGCTATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTAc                                                                                                                                            >  2:45297/1‑249 (MQ=255)                                                                                                             aGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGc                                                                                                                                      >  1:104660/1‑251 (MQ=255)                                                                                                                gAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTg                                                                                                                                                                              <  2:325500/208‑1 (MQ=255)                                                                                                                gAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTg                                                                                                                                                                              >  1:325500/1‑208 (MQ=255)                                                                                                                 aaaGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGAt                                                                                                                                  >  1:166871/1‑251 (MQ=255)                                                                                                                   aGGGCTGGTGAAATTACGCGGCACGGTGGAAGCGCCACCACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCCAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTACATATCAACCTCGTCAC‑CCTGCATGCCAACGGCAGCGAAAAGCAGCCCGTGCTACTGTTGCGGTTTc                                                                                                                                >  2:17952/1‑251 (MQ=255)                                                                                                                   aGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTc                                                                                                                                >  1:69739/1‑251 (MQ=255)                                                                                                                    gggCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATAGCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCa                                                                                                                               <  2:175338/251‑1 (MQ=255)                                                                                                                         ggTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAggg                                                                                                                            <  1:459399/249‑1 (MQ=255)                                                                                                                          gggAAAGTACGCGGCAGGGGGGAAGCGCCACACCTACTGGCGGACATCACCGCGCGCGGCCTGCGCTGTGAGGAACTTTCCGTGGCGCAGGTTCTCGGGGAAGGCGATAACAAATCCACCGATGAGATCGCCGGGAAACTCGACGTGCGCGTGGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCGGAATGCGAAAGGCGGCGAAAAGCAGCACGAGCTACAGTTGCGGAGTCAGGGGGa                                                                                                                         <  2:451778/249‑1 (MQ=255)                                                                                                                             aaaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGcc                                                                                                                      >  1:122494/1‑251 (MQ=255)                                                                                                                                              ggAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCtt                                                                                                      >  2:87397/1‑250 (MQ=255)                                                                                                                                                        aaaCTACTGGCAGATAGCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCTTTGAATGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATGCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAgg                                                                                           <  2:464630/250‑1 (MQ=255)                                                                                                                                                             actGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGtt                                                                                      >  1:349288/1‑251 (MQ=255)                                                                                                                                                                  cAGATAACCCCGCGCGCGGCCGGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGAt                                                                                 <  2:349288/251‑1 (MQ=255)                                                                                                                                                                       atCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCaa                                                                            <  1:476415/251‑1 (MQ=255)                                                                                                                                                                        acACCGCGCGCGGCCTGCGCGGGCAGGAACTGTCCGTGGCGCAGGTTCGCGTGGAAGGCGATAGCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCaaa                                                                           <  2:359389/250‑1 (MQ=255)                                                                                                                                                                         cACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACTTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAaag                                                                          <  1:263097/251‑1 (MQ=255)                                                                                                                                                                                  cgGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTTCAGTTGCGGATTCAGGGGGAGGCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCt                                                                 <  1:216103/251‑1 (MQ=255)                                                                                                                                                                                      ctgagcTGGCAGGAACTTTCCGTGGCGAAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTAGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGAAAAGAAGAACGCTggt                                                              >  2:212447/5‑249 (MQ=255)                                                                                                                                                                                      cTGCGCTGGCAGGCACTTTCCGTGGCGCAGGTTCGCGGGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAAAATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGATACAGTTGCGGAATCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATAGCAAAGAAGAACGCTgggg                                                             >  2:163992/1‑249 (MQ=255)                                                                                                                                                                                              gCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATTCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGAGTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAActc                                                     <  1:438748/251‑1 (MQ=255)                                                                                                                                                                                                   aaCTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGc                                                >  2:61232/1‑251 (MQ=255)                                                                                                                                                                                                      tttCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAAt                                             >  1:333540/1‑251 (MQ=255)                                                                                                                                                                                                                    gTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGAc                               <  1:90471/251‑1 (MQ=255)                                                                                                                                                                                                                    gTTCGCGTGGAAGGCGATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGAc                               <  1:262716/251‑1 (MQ=255)                                                                                                                                                                                                                                 gcgttaacAAATCCACCGATTAGAACGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTTGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGATCTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCg                  <  2:80781/244‑1 (MQ=255)                                                                                                                                                                                                                                   gATATCAAATCCACCGATTAGATCGCCGGGAAACTCGACGTCCGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGCGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCACTTACAGACGCCGGTTGGCCCGTg                >  1:378153/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaaTCCACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCt          <  2:92464/251‑1 (MQ=255)                                                                                                                                                                                                                                             ccACCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGAcc      <  1:87397/251‑1 (MQ=255)                                                                                                                                                                                                                                               aCCGATTAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGACCcg    <  1:29110/251‑1 (MQ=255)                                                                                                                                                                                                                                                 cGATTAGATCGCCGGGAAACTCGACGTACGCGCTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCGGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGGCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGACCcgcg  <  2:122494/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        GCCAAACAGTGCCGAAGTGAAAGGCGAGCTGGGGGTAAAAGATCTCAATCTTGATGCCACCATCAACGCGCCGGGGCTGGATAACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAC‑CCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAGCCTGTCTCCGGGCAGCTTAATCTGGCAGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGACCCGCG  >  NC_000913/4445680‑4446170

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