Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,190	0	A	G	31.0%	46.9 / 28.5	64	intergenic (+396/+246)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (23/17);  new base G (11/7);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAC  >  NC_000913/4295892‑4296481                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                       ggattAGCGGGATTTTGTAGGTCGGATCAGGCGTTTATGCCGCATCCGACATCAACGCCGTAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATAAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAAc                                                                                                                                                                                                                                                                                                      <  1:1136893/298‑1 (MQ=255) aaTTTTGCGTAATTTTGTAGGTCTGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAAc                                                                                                                                                                                                                                                                                                      <  2:331665/300‑1 (MQ=21)       gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGGTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCACCCCCTGATCGCG‑CGCTTAAagccttt                                                                                                                                                                                                                                                                                                >  2:841828/1‑297 (MQ=255)             ttttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTATGTCGGATAATGCGTTAATGCCGCATCCGACATTAACGCCAGAG‑GCGACGCTTACCGCGTCttgac                                                                                                                                                                                                                                                                                            >  1:90873/1‑296 (MQ=255)               tggactgatggATAACGCGTATATGCCTCATGCTACATGCA‑GCCTGGTAGCGAGCCTTGTCGGGTCTTACCCTGACTCAGCCAGCCAGCTCAATTGCCTGTGTTAGCGTGCTTATGTAGGTCTGATAAG‑GCGTTTACGCCGCCTCCGACATCAACG‑CCTGTTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑AAGACGCGTTAAGCGTCGCAT‑CAgg                                                                                                                                                                                                                                                                                        <  1:888908/291‑1 (MQ=2)                    ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGACGCGTCTCATCAGACCTACGCCAGACAACGCAATAGACCGATCTAGCGGTATTTTGTAGGTCGGATAAGGGGGTGATCCCGACTCAGCCCTCAACGCATGAT‑GCGGCACTGAGCGCGTCTTAT‑TAGGCCTAcc                                                                                                                                                                                                                                                                                  >  1:801710/1‑300 (MQ=11)                         gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGATGCTTAACGCGTCTTAT‑CACGCCTaacatacg                                                                                                                                                                                                                                                                             >  2:708127/1‑294 (MQ=25)                           tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTACGCCCGCCAGCGCAATAGCATACTTTAGCGTTATTTTGTAGGTCGGAAAAGGCTGTTCTGCCACATCCGACATCAACGCCTGAT‑TCGACTGTTCACGCGTCTTTT‑CAGGCCTACGCCCGAcc                                                                                                                                                                                                                                                                           >  1:1069723/1‑300 (MQ=11)                            aaGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTACGGCGTTTATGCCGCATCCGACATCAACGCCTGCT‑GCGCCCCTTAACGCGTCTTAT‑CAGGCCTCCGCCAGAccg                                                                                                                                                                                                                                                                          >  2:134030/1‑299 (MQ=255)                                      tGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGCCATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGGTTATGCCGACTCCGAACTAAACGACTGAT‑GCGACGCTTCAAGCTTCTTTT‑CAGGCCTACGCCTGCCAGCTCAATAGcc                                                                                                                                                                                                                                                                >  2:796160/1‑301 (MQ=11)                                               cGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCCGCGCATTAGCCTGATTTCTCTTGATTTTTTAGGTCGAATAAGTCGTTTATGCGGCTTCCGGATTAAACGCCTGAT‑TCGACGAGTGACGCGCTTAAT‑CCGGCATCCGCATGACAgcgctagagattctcatcac                                                                                                                                                                                                                                                       >  1:381738/1‑285 (MQ=11)                                                ctggcgtaggCCTGAT‑AAGACGCGGCAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTATGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTACCGTTATTTTGTAGGTCGAATAAGGCGTTTACGCCGAATCCGACATCCAAGCCGGAG‑GCACCGCCTTTCGAG‑CTTGTACAGCCCTACGCCAGACAGGTCAAtcgacgtactcttac                                                                                                                                                                                                                                                      >  1:609483/10‑287 (MQ=11)                                                  cATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGCTGCGATGCTTGCCGCGTCTTATCAGCCCTACGCCACACAGCGCAATAGCCTGATCTAGCGTGCTTTTGTAGGTCGTATAAGACGTTGACGCCGCTTCCGACATACACGCCTACT‑GCACCGCTTAACTCCTGTCAT‑CAGGCCTACGTCATACAgcgtcacatcctcagtacggctg                                                                                                                                                                                                                                                    >  1:275548/1‑281 (MQ=11)                                                    tCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGTTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCT‑GCGACGCTTAACGCGTCTTAT‑CAGCCCTACGCCAGCCAGCGAAATAGCCTGATTTAgtggaat                                                                                                                                                                                                                                                  >  2:1122267/1‑295 (MQ=255)                                                     cAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGCCACCGCAATAGGCTGATTTCGCGTGGtt                                                                                                                                                                                                                                                 >  2:853487/1‑301 (MQ=21)                                                      aaCGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAttt                                                                                                                                                                                                                                                >  2:969710/1‑301 (MQ=34)                                                                 gtgagtcTTAACGCGTCTGTTCAGGCCTACGACAGCCAGCGCAATAGCCTGATTTAGCGGGTTTTTGTAGGTCGGATAAG‑GCTTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCg                                                                                                                                                                                                                                       <  2:844918/294‑1 (MQ=255)                                                                 gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAATG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCAAGCGCAACAGCCTGTCTTAGCGTGACTTTGTAGGTCGGATAAGGCGTTTATGCCTCATCCTACATCAGAGCTTGAT‑GCGCGGCTTGACGCGGCTTAT‑CAGGCCTGAGCTCGGCAGCGCAATAGCCTCAtatcccgcgactctgtagattgg                                                                                                                                                                                                                                      >  1:185127/1‑279 (MQ=11)                                                                   gACGATTTACGCGGCTATGCAGGCTTAGGCCAGACAGCGAAATAGCCTAATTTAGCGTGATTTTGTCGGTCGGATAAG‑GCGTTTACGCCGCATCAGACATCAACG‑CCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAt                                                                                                                                                                                                                                    <  2:165490/301‑1 (MQ=11)                                                                    acgacgaacgctgcgtATCAGTCCTCAGCCAGACAGCGCAATAGCCTAATTTAGCGTAACTTTTTAGGGCGGATAGG‑GCGTTTCCGCCGCCTCCGACATCGACG‑GCTGATGCGACGCTTTCCGTGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTTATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATa                                                                                                                                                                                                                                   <  2:1087391/286‑1 (MQ=11)                                                                     cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCTGATa                                                                                                                                                                                                                                   >  2:366216/1‑300 (MQ=25)                                                                       gatgcggggtgttctCAGGCCAACACCCGACAGCGCCACAGCCCGATTTAGAGGTCTTTTGTAGCTCGGAATAG‑GCGGGTACGCGGCATCCGCCATCAGCG‑CCCGATGCAACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAgg                                                                                                                                                                                                                                <  2:145023/287‑1 (MQ=11)                                                                          gatccttCTTAGCAGGCCTCCGCCGGATAGCTCATTTGGTTGATTTTGCGTGATTTTTTAGGTCGGTTAGA‑GCCGTTACGCGCCAGCCTTCCT‑AACGCCCTGTTGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGACCGCGGAATGGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCCGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                                                                                                                                                                                              <  1:738325/294‑1 (MQ=11)                                                                          aaCGCGTCTTATCAGGACTACACCCTACAGCGCAATAGCCTGTGTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCAGCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                                                                                                                                                                                              <  2:801931/300‑1 (MQ=255)                                                                           aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTCGGTCGGCAAAggag                                                                                                                                                                                                                              >  2:738534/1‑297 (MQ=18)                                                                            cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGACTGTT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTCCGCCCAACAGCGCAATATCCTGCTTTCGCGTGAtttttgagtgagtatcaggctttt                                                                                                                                                                                                                           >  2:411296/1‑281 (MQ=11)                                                                            cgcgGCTTATCAGGCCTAGTCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGttt                                                                                                                                                                                                                           <  2:426176/301‑1 (MQ=17)                                                                               gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATACCCTGATTTAGCGTGATTTTGTAGGTCGGGTAAGGCGTTTacg                                                                                                                                                                                                                        >  2:1171096/1‑299 (MQ=17)                                                                                      cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCTACGCTTAACGCGTCTTAC‑CAGGCCTACCACAGACAGCGCAATAGCCGTATTACGCGTCGTTTTGTTGGTCGGACACGGCGtttcagccccctc                                                                                                                                                                                                                 >  1:527873/1‑291 (MQ=11)                                                                                        gtgcaacgccaaaagacgaactAGCCTGATTTAGCGGTGTTTGAAAGGGCGGATAAG‑GCCTTTACGCGGCACCGCACATCAACG‑CCTGATGCGACACTTGCCGCGTCTTATCAGGCCTAGGCCAGACAGCGCAATAGCCTAATTTAGCGTGATGTTGTAGCTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                                                                                               <  2:275604/280‑1 (MQ=11)                                                                                             cctgcataccggtggagaccaggatgttgggTGATTTTGTGGGTCGGATAAG‑GCCTTTACGCCGAATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCGTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATc                                                                                                                                                                                                          <  2:609659/271‑1 (MQ=11)                                                                                               gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGTATGATTTCGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GTGACGCTTAACGCGGCTTAT‑TACGCCTACTCCCACCATCGCAAAAGCACGATTTAGCGTCATGTCGTTGGTCGGCTAACGCGTTTAcgccgcttgctgcaccac                                                                                                                                                                                                        >  1:156710/1‑289 (MQ=11)                                                                                                ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCTCATCCGACACCAACGCCTGCT‑GCGCCGCTTACCGCGTCTTAT‑CAGCCCTACGCCCGCCCGCGCAAGAGCCTGCTCTAGCGTCATTTTGTCGCTCGGaacgccctgttagacctcatcatccacccac                                                                                                                                                                                                       >  1:561287/1‑271 (MQ=11)                                                                                                   gACAGCTAAATGGCCTTATTTGGCGTTATTTTGTGGGTCTGATAAG‑TCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                    <  2:1101919/301‑1 (MQ=11)                                                                                                    aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCCCAACACCCTGAATTCGCGGTAGTTTTTGGGCCGAATAATGCGTTTCTGCCGCATCTGACATCaatcccc                                                                                                                                                                                                   >  1:1049734/1‑296 (MQ=11)                                                                                                                 cTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCttg                                                                                                                                                                                      >  2:554313/1‑300 (MQ=18)                                                                                                                       tAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGACGCGt                                                                                                                                                                                >  2:224483/1‑301 (MQ=21)                                                                                                                        agggtggttgtggaggcgcgacaaG‑GCGTTTCCACCGCATCCGCAATCAGCG‑CCTGATGAGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                               <  2:1107416/279‑1 (MQ=11)                                                                                                                           gTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTCCGCCGCATCCGACATCCATGCCTGATGCGACGCTTGCcgcggcgtt                                                                                                                                                                            >  2:592463/1‑296 (MQ=255)                                                                                                                                cgttgatgTCGGATGCG‑GCATAAACGCCTTATCCGACATCAACG‑CCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTAC‑CAGGCCTACGCCAGACCGCGCAATAGCCTGATTTCGCTTGATTTTGTACGTCGGATCCGGAGTTTAGTCCGGATACGATATCCACAACTGAGGCGTCACTGATcgcggcgtactcgg                                                                                                                                                                       >  1:942925/8‑291 (MQ=11)                                                                                                                                 ttGTAGGTGGGATAAG‑GTGTTTGTGCCGCATCAAACAGCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCACACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                                                                       <  1:773231/300‑1 (MQ=11)                                                                                                                                 ttGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCCTCCTACAGCAATGCCTGTCGCGACCCTTGCCGCGTCTTTTCAGGc                                                                                                                                                                      >  1:696258/1‑301 (MQ=11)                                                                                                                                 ttGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCAGCATCCGACATCAATGCCTGATCCGACGCTTGCCGCGGCTTATCCGGc                                                                                                                                                                      >  1:707249/1‑301 (MQ=255)                                                                                                                                  tGTAAGTCGGAGAAGTGCG‑TTACGTCGCAGCCGACATCAGTG‑CACGGTGGGACTCGTGCCTGGTCTTATTAGGCTTACGCCAGACAGCGCAATAGGCTGTTTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGGCATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                      <  1:1170781/300‑1 (MQ=1)                                                                                                                                   gTAGGGCGGAAAAG‑GCTGTTACGCCGCATCCGACATCAACG‑CCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                     <  2:1002330/300‑1 (MQ=255)                                                                                                                                               aG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACTCCACACAGCGCAATAGCTTGCTTTAGCGTGTTTTTGTAGGTCGGCTAAGGCTTTTACGCCTCCACCGACATAACCGCCTGACGCGACGCTTGCCGCGTCTTAACAGGcccatcttcaccgtt                                                                                                                                                        >  1:935883/1‑288 (MQ=11)                                                                                                                                                                                                                                                                   ggCGTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                      >  2:1130557/1‑301 (MQ=2)                                                                                                                                                                                                                                                                   cctcgtcggccctctctcgagctCAACGCCTAGT‑GCGCCGTTTGCCCCGTCTTGT‑CAGGCGTACGCCAGACACCGCAATCGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAGTGCGTTTATGCCGCATCCGACATCAACGCCTGAGGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                      <  1:364432/279‑1 (MQ=11)                                                                                                                                                                                                                                                                    gcgtgcacgCCGCAGCCGACCTCCATGCCCTAT‑GCGACGCTCACCGCGTTTCAC‑CCGGACTATCTCAGCCAGCGAAATAGCCTGAGTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                     <  2:736508/293‑1 (MQ=11)                                                                                                                                                                                                                                                                    gCGTTTATGCCCCATCCGACATCACCGCCTGAT‑GCGACGGTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                     <  2:762810/301‑1 (MQ=255)                                                                                                                                                                                                                                                                      gTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTTATGCGACGCTTCCCGCGACTTAACAGGACTATCTTAACAGTTGGtt                                   >  2:1059977/1‑301 (MQ=11)                                                                                                                                                                                                                                                                      gTTTACGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                   >  2:261809/1‑301 (MQ=12)                                                                                                                                                                                                                                                                              gcgctcccacacaaccGCCTGAT‑GCGCCCCTTGCCGCGTCTTAT‑CAGGCCTACGACAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                           <  2:140722/286‑1 (MQ=11)                                                                                                                                                                                                                                                                                  cccGACATCCACGCCTGAT‑GCGACGCTTGCCTCGTCTTAT‑CAGGCCTACGCCAGACAGCGAAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                       <  2:936132/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                   ccGACATCAACGCCTGAT‑GCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCCGCGCAATAGCCTGCTTTAGCGATATTTTTTAGGTCGGACAAGGCGTTTACCCCGCATCCGACATCAATGACCGGTGCGCCGCCTGCCGATTTTTATAATGCCCACCTTAACAGTCTCGTAATTTTTcgcgtc                      >  1:140694/1‑296 (MQ=11)                                                                                                                                                                                                                                                                                    cgagtccaccgccgact‑GCGTCGCTTGCCTCGTCTTAT‑CAGGCCTACGCCAGACAGCTCAATAGCCTTTTTTAGCGTTATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                     <  1:853254/285‑1 (MQ=11)                                                                                                                                                                                                                                                                                    cgaattCAACGCCTGAT‑GCGTCGCGTTCCGCTTCTTAG‑CAGGCCTACGCCGGACAGGCCAATATGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCGTCCGACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                     <  2:366715/296‑1 (MQ=14)                                                                                                                                                                                                                                                                                       ggcgtaggCCTGAT‑AAGACGCGGCAAGCGTCGCAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCAACCGACCTCAATGCCTGATGCGACGCTTGCCGCGTCGTATCAGGCCTATCTTAACCGTGTGTTAATTTTTCAAACttttt                  >  2:884699/8‑300 (MQ=11)                                                                                                                                                                                                                                                                                          cAACGCCTGAT‑GCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCATACATGCCTGATGCGACGCTTGCCGCGTCGTCTCAGCCCTACCTTAACCGTTGGTTAATTTTTCAAACTCTTcaa                >  1:391255/1‑299 (MQ=11)                                                                                                                                                                                                                                                                                          cAACGCCTGAT‑GCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgggatcggacgagcacacgtctgg               >  2:364514/1‑279 (MQ=11)                                                                                                                                                                                                                                                                                                cTGAT‑GCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAATCTCTTCacac              >  2:131199/1‑296 (MQ=255)                                                                                                                                                                                                                                                                                                  gtt‑GCCACGCTTGTCGCGTCTTAT‑AAGGCCTACGCCAGAGAGCGCTATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa        <  1:366134/298‑1 (MQ=14)                                                                                                                                                                                                                                                                                                      gCGACGCTTGCCGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACACGGTATCa     >  2:1093843/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                        aaCGCTTGCCGCGGCTTAG‑CAGGCCTACGCCAGACAGCGCAATCTCCTGGTTTAGCGTGATGTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc  <  2:185169/300‑1 (MQ=14)                                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                       GATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAG‑GCGTTTACGCCGCATCCGACATCAACG‑CCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAT‑GCGACGCTTAACGCGTCTTAT‑CAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAC  >  NC_000913/4295892‑4296481

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑