Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,337,723	0	A	C	28.6%	103.9 / 13.4	74	G566G (GGT→GGG)	gyrA	DNA gyrase subunit A
Reads supporting (aligned to +/- strand):  ref base A (38/12);  new base C (1/19);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 3.45e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCGCGATCT  >  NC_000913/2337474‑2337970                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          cGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTAc                                                                                                                                                                                                                                                          <  2:7063/251‑1 (MQ=255)    cATGAAGACTTTAACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGAC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGATGGAGAAGCACAGAATATGGTCGTGAGTG‑TCCGCCACC‑AGCAGT‑CGGTCG‑ATAAAACCTTCTTCTTT‑AATACGGGCGGCAGATTTACCtt                                                                                                                                                                                                                                                       >  1:545222/1‑251 (MQ=255)      tGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑c                                                                                                                                                                                                                                                    >  1:7745/1‑251 (MQ=255)       gAAGACTGTCACGCCGTCTTCAAACTCGGTCACTGGCAGGATCGCAGGGATACGTGCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGTCCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCTACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑GTCGCCACC‑AGCATT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATGCGTGCGTCAGATTTCCCTTT‑c                                                                                                                                                                                                                                                    <  2:669884/250‑1 (MQ=255)          gACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAGGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑Gc                                                                                                                                                                                                                                               >  1:585385/1‑251 (MQ=255)          gACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTACGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGCAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATCCGTGCGACTGCTTTACCTTT‑CCC‑Gc                                                                                                                                                                                                                                               >  1:106947/1‑251 (MQ=255)           aCTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑Gcc                                                                                                                                                                                                                                              >  2:258194/1‑251 (MQ=255)             tttCACGCCTTCTTCAAAATCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCATGCTCCAGCGGCAGCAGGATGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑CGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCac                                                                                                                                                                                                                                            >  2:329229/1‑251 (MQ=255)               tCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCCGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCacg                                                                                                                                                                                                                                           >  1:324612/1‑250 (MQ=255)                        cttcAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCt                                                                                                                                                                                                                                 >  1:8707/1‑251 (MQ=255)                           cAAACTCGGTCACGGGCAGGATCGCAGTGATCCGTCCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGCCCGCGCGCGCCACGAGGGGCTTCCGGA‑AACTAATAAACTTTCATCGAATAGACGCGACCAGGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCCCC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGGGGCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCg                                                                                                                                                                                                                              <  1:17044/251‑1 (MQ=255)                                        gggCAGGAGCGCAGGGATGCGTTCGTCCGCCTCCGTCGCCAGAAGGTTGACAATCGGACGCCCGCGCGCGCCACGAGTGGCTTCCGGA‑AATTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCAAAGAATAGGGTCGTGAGTG‑TCCCCCACC‑AGCAGT‑CGGTCG‑ATAAATTCTTCTTCTTA‑AATCCGTGCGGAAGATTTACTTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCg                                                                                                                                                                                                                 <  1:520715/250‑1 (MQ=255)                                         ggCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCTCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGCGACGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑CATACGTTCGGCCGATTTCCCTTT‑CCC‑GCCCCGACGCTGCGCTTCGTATTCAGAAATCgg                                                                                                                                                                                                                >  2:92807/1‑251 (MQ=255)                                                        attcGTTCGCCCGGCCCCAGCGGCAGCAGGTTGACGACCGGACGACCGCGCCCGCCACGGGTGGCT‑CCGGCCAACTGATAAACTTTCATCGAATAGACGCGACCCCGGCGGGAGAAGCACAGAATATGGCCGTGAGTT‑TTCCCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATCCGTGCGGCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTa                                                                                                                                                                                                 <  1:252846/248‑1 (MQ=255)                                                        aTACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGCCGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGCCGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGCCGCTTCGCTTCTTATTCACAAAGCGGATGCTACTTAACGt                                                                                                                                                                                                  >  2:217088/1‑250 (MQ=255)                                                         tACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATTCGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCCGAAAGCGGCTGATACTTAACGt                                                                                                                                                                                                  >  2:665760/1‑249 (MQ=255)                                                                  ccTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGC‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                                                                                                                        >  2:552836/1‑250 (MQ=255)                                                                    tGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                                                                                                                                                        >  2:194272/1‑248 (MQ=255)                                                                      ccccAGCGGCGGCAGGTTGACGAGCGGACGACCGCGCGCGCCACGAGTGGCTTCGGGC‑ACGTGATAAACTTTCATCGAAGGTACGCGCCCCCGGCTGGAGAAGCAAAAAATATGTTCGTG‑GTGGTCCGCCACC‑ATCAGT‑CGGTCT‑AGAAGGTCTTCTTCTTT‑AAAACGGGCGTCTAATTTCCTTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTgagaga                                                                                                                                                                                   <  1:661675/249‑1 (MQ=255)                                                                            cggcatcaggtggtctaTCGGACGACCGCGCCCCCCCCGAGGGGCTTCCGGC‑ATCTGATAAACTTTCATCAAAAAGACGCCACACCGGCTGGATAAACCCAGAATATGGTGGTGAGTG‑TTCGCCCCC‑AGCAGC‑CGGCCA‑ATAAGTTCTTTTTCTTT‑ATTTCGGGCGCCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGGAAGCGTCTCATACTTAACGTAACCCGGGTGAGAGAGCGTCa                                                                                                                                                                             <  1:236515/235‑1 (MQ=255)                                                                            cGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCTTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCATATTTACCtt                                                                                                                                                                                                                                                       <  2:39491/179‑1 (MQ=255)                                                                            cGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCtt                                                                                                                                                                                                                                                       >  1:39491/1‑179 (MQ=255)                                                                                    ggTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACcac                                                                                                                                                                      >  1:287288/1‑250 (MQ=255)                                                                                           gATCGGACGACCGGGCGGGCCAGGGGGGGCTGCCGGA‑AACGGATAAACTTTCAAAGAAAAGACGCGCACACGGCGGGAGAAGCACAGAAGAGGGGCGTGAGTG‑TTCGCCACC‑AGCAGG‑CGGTCA‑AAAAAGTCTTCTTTTTT‑AATACGTGCGGAAGATTAACCTTG‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                                                                                                                                              <  1:92807/251‑1 (MQ=255)                                                                                               ggACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCCACGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATACAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTgggg                                                                                                                                                          >  1:46377/1‑250 (MQ=255)                                                                                                       gcgcgcgcCACGGGTGGCTGCCGGC‑AACTAATAAACTTTCATCGAATGGACGCCACACCGGCGGGAGAAGCAAAGAATTGGGTCGGGGGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTTTTT‑AATACGGGCGGCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGAt                                                                                                                                                  <  1:516438/251‑1 (MQ=255)                                                                                                       gcgcgcgcCACGAGTGGCTTCCGGC‑AACTGATACACTTTCATCGAATAGACGCGACCACGGCTGGCGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACATTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATc                                                                                                                                                      >  1:246988/1‑247 (MQ=255)                                                                                                            cgcCACGGGCGGCTTCCGGC‑AACTTATACAGTTTCATCGAAGAGACGCCCCCCCGGCTGGAGAAGCAAATAAGATGGTCGTGATTT‑GTCGCCCCC‑GGCAGT‑CGTTCG‑AAAAATTCTCCTCCTTT‑AAGACGTTCGGCAGATTCACCTGC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCtt                                                                                                                                               <  2:566145/249‑1 (MQ=255)                                                                                                                ccggggggtttCCGGA‑AACTAATAAATTTTCATGGAATAGCCGCGCCCCCGGGGGAAGAAGAAAAAAATATGGTCGTGAGTG‑TTCGCCCCC‑AGCAGT‑CGGTCA‑ATAAAGTTTTCTTTTTT‑AATCCGTGCGGCAGATTTACCTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGt                                                                                                                                         <  2:276416/242‑1 (MQ=255)                                                                                                                 cgggTGGTTTCCGGC‑AACTGAAAAACTTTCATCGAATAGACGCGACCGCGGCGGGAGAAGCAAAGAATATGGTCGTGAGTG‑TTCCCCACC‑AGCAGC‑GGGTCG‑ATAAAGTCTTCTTCTTT‑AAGACGTGCGGCATATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGtt                                                                                                                                        <  1:552836/248‑1 (MQ=255)                                                                                                                       gCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAACGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGCCGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCCCGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTc                                                                                                                                  >  2:100179/1‑251 (MQ=255)                                                                                                                           cccggaaacgtaaaaaCTTTATTCGAAGAGCCGCGCCCCCGGCTGGAAAACCACAGAATATGGTCGGGATTG‑T‑CCCCACCCAGCTGT‑GGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgc                                                                                                                               <  2:171536/238‑1 (MQ=255)                                                                                                                            cGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATAAGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCGTCCAGGTTGATGTCTgcgc                                                                                                                             >  2:456385/1‑251 (MQ=255)                                                                                                                             ggC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCt                                                                                                                            >  1:364735/1‑251 (MQ=255)                                                                                                                                  attGATAAACTTTCATCGAATAGACGCGACCACGGCGGGAGAACCACAAAATATGGTGGTGGGTG‑TCCCCCCCC‑AGCAGT‑GGGTCA‑ATAAATTCTTCTTCTTT‑AATACGGGCGGAAGATTTACCTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                                                                                                                        <  1:258194/249‑1 (MQ=255)                                                                                                                                  acggATAACCTTTCAGAGAATAGCCGCGACAACGGCGGGAGAAGCAAAGAATATGGTCGTGATTG‑TGCGCCACC‑AGCATT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                                                                                                                        <  2:420986/248‑1 (MQ=255)                                                                                                                                  aCTGATAACTTTTCAGCGAGAAGACGCGCCCGCGGTTGAAAAAGCACAGAATAGGGTCGTGGGTG‑TTCGCCACC‑AGCAGC‑CGGCCG‑ATAAATCCTTCTTCTTG‑AATAGGTGCGGCAGATTTCCCTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTg                                                                                                                        <  1:598081/251‑1 (MQ=255)                                                                                                                                    tGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGACTATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGCTCTTCCAGGTTGATGTCTGCGCTGTTGGc                                                                                                                      >  1:246014/1‑251 (MQ=255)                                                                                                                                    tGATAAACTTTAAGCGAATAGCCGCGACCACGGTGGGAAAAGCACAGAATTGGGTCGTGAGTG‑TCCGCCCCC‑AGAAGT‑CGTCCG‑ATAAAGTCTTCTTCTTT‑ATTCCGTGCGTCAGTTTTACTTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTgg                                                                                                                       <  1:452764/250‑1 (MQ=255)                                                                                                                                          aCTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGAt                                                                                                                >  1:76247/1‑251 (MQ=255)                                                                                                                                             ttCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTc                                                                                                             >  1:434740/1‑251 (MQ=255)                                                                                                                                               cATCGAATAGACGCGACCACGGCTGGAGAAGCACAGCATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTATTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTc                                                                                                             >  2:362040/1‑249 (MQ=255)                                                                                                                                                        gACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGtt                                                                                                  >  1:110760/1‑251 (MQ=255)                                                                                                                                                              accccggtgggaaaaGCACAGAATATGGTCTGGAGTG‑TGCCCCTCC‑AGCGGT‑CGGTCG‑AAAAAGTCTTCTTCTTT‑ATCCCGTGCGGCAGATTTTCTTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATGTCAGTACGACGTTTGTCAc                                                                                            <  1:269719/238‑1 (MQ=255)                                                                                                                                                              aCCACGGCTGGAGAAGCACAGAATATGGTCGGGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTGGTATTCAGCAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGAGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGCTTTCCGTACGACGTCTGTCAc                                                                                            >  1:318766/1‑251 (MQ=255)                                                                                                                                                                cACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGCATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGCGAGCGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCg                                                                                          >  1:689810/1‑251 (MQ=255)                                                                                                                                                                 ccggcggGAGAAGCAAAGAATAGGGTCGGGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGAAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGa                                                                                         <  1:194272/250‑1 (MQ=255)                                                                                                                                                                      tGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGGGATTTCAGTACGACGTTTGTCACCGAACTGt                                                                                    >  1:411599/1‑251 (MQ=255)                                                                                                                                                                       ggAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGCTCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGCTTTCAGTACGACGTTTGTCACCGAACTGtt                                                                                   >  2:405187/1‑251 (MQ=255)                                                                                                                                                                       gaagaAGCAAAGAATGTGGTCGGGAGTG‑TTCGCCACC‑AGCGGT‑CGGCCG‑ATAAAGTCTTCTTCTGT‑AATACGTGCGGCATATTTCCCTCT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCTTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGt                                                                                    <  1:604168/248‑1 (MQ=255)                                                                                                                                                                                    aatggtgcGTGAGTT‑TCCGCCACC‑AGCAGT‑CG‑TCGTAAAAAGTCTTCTTCGTT‑AATACGGGCGGAAGATTTTCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACTACGTTTGTCACCGAACTGTTCACGAAccagctc                                                                      <  2:258214/244‑1 (MQ=255)                                                                                                                                                                                      atGGTCGTGAGTG‑TTCGCCGCC‑AG‑AGTTAGTTGG‑AAAAAGTCTTCT‑CTTTGAATACGTGGGGAAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGAGGTCTGCGCTGTTGGCGGTGATTTCAGGACGACGTGTGTCACCGAACTGTTCACGAAccagctcca                                                                    <  2:105237/251‑1 (MQ=255)                                                                                                                                                                                      atGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAAccagctccc                                                                    >  1:222998/1‑250 (MQ=255)                                                                                                                                                                                         gttgtggtgttgcccccccc‑aGCAGT‑CGCTCG‑CTAATGTCTTTTTCTTT‑ATTACGTAGGGCAGGTTACCCTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGGTTGTCACCGAGCTGTTCACGAAccagctccagc                                                                  <  2:490661/233‑1 (MQ=255)                                                                                                                                                                                            gTGATTG‑TTCGCCCCC‑AGCAGT‑CGGTCG‑TAAAAGTCTTCTTTTTT‑AATCCGTGCGGAAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCtt                                                              <  1:534951/251‑1 (MQ=255)                                                                                                                                                                                            gTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGGTTTCAGTACGGCGTTTGTCACCGAACTGTTCACGAACCAGCCCCAGCTCtt                                                              >  2:396343/1‑251 (MQ=255)                                                                                                                                                                                                    gtccgccccacgcagT‑CGGCCA‑AAAAAGTTTTCTTCTTA‑AAGTCGTGGTGCAGTTTTACCTTC‑CCC‑GCCACGACGCCGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCGGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACgg                                                         <  2:365579/239‑1 (MQ=255)                                                                                                                                                                                                       cccacc‑aGCAGC‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGGGCGGCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCAc                                                    <  1:163642/250‑1 (MQ=255)                                                                                                                                                                                                         cacc‑aGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCCGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTTACCGAACTGTTCACGAACCAGCTCCAGCACTTCACGGTTCACtt                                                  >  2:179015/1‑251 (MQ=255)                                                                                                                                                                                                                 ggt‑cggtca‑AAAATGTCTTCTTTTTT‑AATGCGTGCGGTAGATT‑CCCTTTACCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAg                                           <  2:253538/250‑1 (MQ=255)                                                                                                                                                                                                                          caaTAAAGTTTTCTTCTTT‑AATGCGGGCGGCAGATTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCGTCACGGATCACTTCCATCAGACGATCg                                    <  2:364735/249‑1 (MQ=255)                                                                                                                                                                                                                              aaaGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCCTCAGACGATCGgcgc                                >  2:445081/1‑251 (MQ=255)                                                                                                                                                                                                                                    ttcttcttT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGACCACTTCCACCAGACGCTCGGCGCTACCa                           >  1:483261/1‑250 (MQ=255)                                                                                                                                                                                                                                    catcttttcttataccgggggaaGATTTTCCTTT‑CCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTAcc                            <  2:7745/229‑1 (MQ=255)                                                                                                                                                                                                                                         cttT‑AATACGTGCGGCAGATTTACCTTT‑CAC‑GCCACGACGCTGAGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCAAGGTTGATGTCTGCGCTGTTGGCGGGGATTTAAGTACGCCGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAAt                      >  2:295833/1‑250 (MQ=255)                                                                                                                                                                                                                                              aaTACGTGCGGCAGTTTTACCTTT‑CCC‑GACACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTATGCGCTGTTGGCGGTGATTTCAGTACGCCGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTAAAGGGTCCCTTCTCTCAGACGAACAGCGATACCCAGATTACGcc                 >  2:281183/1‑250 (MQ=255)                                                                                                                                                                                                                                               aTACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCa                 >  2:495540/1‑250 (MQ=255)                                                                                                                                                                                                                                                    ggCGGCAGTTTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGtt           <  2:208618/250‑1 (MQ=255)                                                                                                                                                                                                                                                      gggaaGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTcc         <  1:306607/247‑1 (MQ=255)                                                                                                                                                                                                                                                      ggccatttttACCTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTGCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTcc         <  2:222998/244‑1 (MQ=255)                                                                                                                                                                                                                                                       ggCAGATTTACCTTT‑CCC‑GCCAAGACGCTGCGCTTCGTATTCAGAAAGAGGCTGATACTTAACAAAGCCCTGGGGAGAGAGCGTCACGACCCCATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGAGATTTCAGTACGACGTTTGGCACCGCACTGTTAACGACCCAGCTCCAGCTCTTCACGGATCAATTCCATCAGACGATAGGCGCTTGCAAGAATACCCAACAGTTCcg        >  2:512085/1‑251 (MQ=255)                                                                                                                                                                                                                                                          aGATTTTCCTTC‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCGCGa     <  1:101794/251‑1 (MQ=255)                                                                                                                                                                                                                                                            aTTTCCCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTGCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCGCGATc   <  1:263946/251‑1 (MQ=255)                                                                                                                                                                                                                                                             tttACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACATAACGTAGCCCTGGTGCGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACAAGCTCCAGCTCTTCCCGGCTCACTTCCATCAGACGATCGGCGCTACCAAGCATACGCAACAGTTCCGCGctct  >  1:693933/1‑251 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                          CGCCATGAAGACTTTCACGCCTTCTTCAAACTCGGTCACTGGCAGGATCGCAGTGATACGTTCGTCCTGCTCCAGCGGCAGCAGGTTGACGATCGGACGACCGCGCGCGCCACGAGTGGCTTCCGGC‑AACTGATAAACTTTCATCGAATAGACGCGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTG‑TTCGCCACC‑AGCAGT‑CGGTCG‑ATAAAGTCTTCTTCTTT‑AATACGTGCGGCAGATTTACCTTT‑CCC‑GCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCGCTGTTGGCGGTGATTTCAGTACGACGTTTGTCACCGAACTGTTCACGAACCAGCTCCAGCTCTTCACGGATCACTTCCATCAGACGATCGGCGCTACCAAGAATACGCAACAGTTCCGCGATCT  >  NC_000913/2337474‑2337970

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