Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,176,115	0	G	A	22.4%	134.4 / 16.2	77	A58T (GCT→ACT)	tufB	translation elongation factor Tu 2
Reads supporting (aligned to +/- strand):  ref base G (18/41);  new base A (17/0);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.19e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAATGCGACAT  >  NC_000913/4175870‑4176362                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           ttCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGACAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGGGCTGCCCGCGCACTCGACCAGCAAGCT‑AA‑AGCGCCGGGAGAAAAAGCTCgg                                                                                                                                                                                                                                                       >  2:91811/1‑250 (MQ=255)       aaGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CCCGCCGGAAGAAAAAACTCgcgggatt                                                                                                                                                                                                                                                 >  1:102228/1‑244 (MQ=255)       aaGTCGGGCATGTTGCCGGGTTGAGGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGGTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATc                                                                                                                                                                                                                                                 <  2:507502/251‑1 (MQ=255)          tCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCCCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTgggataacc                                                                                                                                                                                                                                              >  1:435369/1‑244 (MQ=255)              gCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCa                                                                                                                                                                                                                                         <  2:402838/251‑1 (MQ=255)                  ggtGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAacac                                                                                                                                                                                                                                     <  2:7496/249‑1 (MQ=255)                  gTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCCCGCATTCGCCCAGATCGAT‑AA‑AGCGCCCGAAGAAAAAACCCGTGGTATC‑CCCATCCacac                                                                                                                                                                                                                                     >  2:79388/1‑251 (MQ=21)                     gCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGGGCTGCTCGGGCATTCGACCCGATCGAT‑AA‑AGCCCCGGAAGGAAAAACTCGTGGTATT‑CCCCTCAacacatc                                                                                                                                                                                                                                  >  2:167973/1‑248 (MQ=255)                      ccTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGCGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCCCCCGCCTTCGACCACATCGAT‑AA‑AGCGCCGGAAGAAAAAACTCGTGGTATT‑ACCATCCACACTtcc                                                                                                                                                                                                                                 >  2:354814/1‑250 (MQ=18)                        tGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGCGGGACAATCGATGTCTAAAGAAACGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGAGCTTCTCGCGCCTTCGGCCCGATCGAT‑AA‑CGCGCCCGGAGAAAAAACTCGTGGTATC‑ACCCTCAACCCTtctcc                                                                                                                                                                                                                               >  1:449063/1‑250 (MQ=255)                          gTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCCTTCGCCCAGATCGGT‑AA‑AGCGCCGGCAGAAAAAGCTCGTGGTATC‑AACATCAACACTTCTCACg                                                                                                                                                                                                                             >  1:172639/1‑251 (MQ=255)                          gTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCCCATTCGACCAGATCGAT‑AA‑ACCGCCGGGAGAAAAAGCTCGTGGTATC‑ACCATCAACACCTCTCACg                                                                                                                                                                                                                             >  1:394093/1‑251 (MQ=255)                                gtggtgATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAAt                                                                                                                                                                                                                       <  1:519412/251‑1 (MQ=255)                                                         gtCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCCCGCCATCGGCCAGATCGAT‑AA‑CCCGCCCGAACAAAAAACTCCCGGTATCCACC‑TCCACACTTCTCACGTTTGATACGACACCCCCACCCgcccctcccc                                                                                                                                                                                              >  2:151152/1‑242 (MQ=17)                                                             tagaggtgcaagcggtGTCTAAAGAAAAGTTTGAACGTACAAATCCGCACGTTAACGTCGGTACGATCGGCCACGTTGACCATGGGAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CTCGCCAGAATGAAATGCTCGGTGTATC‑ATCCTCAACACTTCTCACGTTGACTACGACACCCCGACCCGTCTCTACGca                                                                                                                                                                                             <  2:524219/233‑1 (MQ=14)                                                              agagGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTTCTCGCGCCTTCGACCAGAACGAT‑AA‑CGCGCCGGAAGAAAAAACTCGTGGTATC‑ACCCTCAACACTTCTCCCGTTGAATACGAAACCCCCACCCGTCACTCCGCacacc                                                                                                                                                                                         >  1:293058/1‑250 (MQ=17)                                                                    aCAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAACTCGTGGGATC‑ACCATCAACCCTTCTCCCGTTGAATACGAAACCCCCACCCCTCACTACCCACACCTAGAC‑t                                                                                                                                                                                  >  1:58770/1‑251 (MQ=14)                                                                     cAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAACTCGTGGGATC‑ACCATCAACACTTCCCACGTTTAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑Tg                                                                                                                                                                                 >  1:182105/1‑251 (MQ=255)                                                                     cAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGGGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGGAGAAAAAACTCGTGGGATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑Tg                                                                                                                                                                                 >  1:404664/1‑251 (MQ=255)                                                                         cGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCCCGCGCCTTCGACCAGATCGAT‑AA‑CGCGCCGGGAGAAAAAACTCGTGGGATA‑CCCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGACCTGccc                                                                                                                                                                              >  1:516539/1‑251 (MQ=14)                                                                           aTGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCCTTCGACCCAATCGATGAA‑C‑CCCCGGAAGAAAAAAATCGTGGGATC‑ACCCTCAACACTTCTCCCGCTGTATACGACAACCCGACCCCGCACCACGCCCACGTAGAA‑CGCCC‑ggg                                                                                                                                                                          >  2:382725/1‑251 (MQ=14)                                                                                 aaaGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATGCGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCAc                                                                                                                                                                      <  1:500698/249‑1 (MQ=255)                                                                                    gAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAACCAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGGGCTGCTCGCGCATTTGGCCCGAACGGT‑AAACCCCCCGGA‑GACAAAACTCGTGGGGTC‑ACCCTCAACCCTTCTCCCGCTGGATACGGCCCCCCCACCCGCCCCTACGcacccctataccaccccgggggcgccgc                                                                                                                                                                   >  2:131889/1‑226 (MQ=255)                                                                                        aGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTAt                                                                                                                                                             <  1:111163/251‑1 (MQ=255)                                                                                           ttGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCAATCGACCAGCTCGAT‑AA‑CCCGCCGGGAGAAAAAACCCGGGGTATC‑CCCACCAACACATCTCACGTTGAATTCGACACCCCGACCCGTCCCTAAGCCCACCTAGAC‑TGCCC‑CGGGCACGCCCACCATGtt                                                                                                                                                          >  1:61860/1‑251 (MQ=14)                                                                                                              cGTTAACGTCGGTACTATCGGCCACTTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATTTTAAAAACATGATCACCGGtg                                                                                                                                       <  2:512436/251‑1 (MQ=255)                                                                                                              cGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGtg                                                                                                                                       <  2:143281/251‑1 (MQ=40)                                                                                                               gTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGtgc                                                                                                                                      <  1:86743/251‑1 (MQ=40)                                                                                                                   aCGTCGGTACTAACGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCTCCGGTGCTgcg                                                                                                                                  <  2:58770/251‑1 (MQ=255)                                                                                                                       cGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCCTTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAACTCGTGGTATC‑ACCATCAACCCTTCTCACGTTGGATACGACACCCCGCCCCCTCCCCACGCACACGTAGAC‑CGCCC‑CGGGCACACCGCCTCTGTTAAAAACACGATCCCCGGGGCTGCGCAgg                                                                                                                              >  2:68412/1‑250 (MQ=14)                                                                                                                              aTCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGAc                                                                                                                         <  2:435369/249‑1 (MQ=40)                                                                                                                                cGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGgc                                                                                                                      <  1:91811/250‑1 (MQ=40)                                                                                                                                cGGCCACGTTGACCAGGGTAAAACAACGCGGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGAACAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGgc                                                                                                                      <  2:427215/250‑1 (MQ=255)                                                                                                                                   ccACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGa                                                                                                                  <  1:79388/251‑1 (MQ=40)                                                                                                                                       gTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCt                                                                                                              <  2:3538/251‑1 (MQ=40)                                                                                                                                         tGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGAACGCT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATA‑ACCCTCAAAACATCTCACGTTGAATAAGAAACCCCGACCCGTCCCTACGCACCCGTGGAC‑AGCCC‑GGGGGACCCCGGCCATGTTAAAAACATGACCAACGGGGGTGCGCAGCTGGGCGGCGCGATcccgg                                                                                                            >  2:2640/1‑248 (MQ=14)                                                                                                                                          gACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTgg                                                                                                            <  1:101420/250‑1 (MQ=40)                                                                                                                                             cATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCCGATAGAT‑AA‑CGCGCCGGAAGAAAAAACTCGTGGTATC‑ACCATCAACACTTCTCCCGTTTAATACGACACCCCCACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGCCTATGTTAAAAACATGATCACCGGGGCTGCGCAGATGGACGGCGCGATCCTGGTAGt                                                                                                        >  1:350950/1‑251 (MQ=14)                                                                                                                                              aTGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGtt                                                                                                       <  2:102228/251‑1 (MQ=39)                                                                                                                                                     aacaacGCTGACCGCTGAAATCACTACCGTACCGGCTAAAACCTACGGCGGTGCTGCTCGCGCATGCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTAGC‑ACCATCGACACTGCTCACGTGGAATACGACGCCCCGACCCGTCGCTACGCACACGTAGAC‑TGCCC‑GGGGCAAGCCGACTATGTTAAAAACATGATCACCGGGGCTGCGCATATAGACGGCGCGACCCTTGTAGTTTCTGCGa                                                                                                <  2:356335/251‑1 (MQ=11)                                                                                                                                                      acaacGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGAc                                                                                               <  2:147131/251‑1 (MQ=38)                                                                                                                                                      acaacGCTGACCCCCGCAATCACTACCGGACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGACGGAAGAAAAAGCTCGTGGTAAC‑ACCATCAACCCGTCTCACCTTGAATACGACACCCCGGTCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATTTTAAAAACATGATCTCCGGTGCTGCTCAGATCGACGGGGCGATCCTGGTAGTTGCTGCGAc                                                                                               <  2:92591/251‑1 (MQ=11)                                                                                                                                                        aacGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCCGAACGAT‑AA‑AACGCCCGAAGAAAAAACTCGTGGGATT‑ACCCTCAACACTTTTCCCGTTGAATACGACACCCCGACCCGCCACCACGCACACGTCGAC‑TGCCC‑GGGGCCCGCCGACCCTGTTAAAAACATGAACACCCGGGCTGCGCAGATTGACGGCGCGATCCTGGTAGCTGCcgcgcgtg                                                                                             >  2:9303/1‑247 (MQ=11)                                                                                                                                                              gACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGc                                                                                        <  1:198252/250‑1 (MQ=35)                                                                                                                                                              gACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCGCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACgg                                                                                         <  2:376817/249‑1 (MQ=25)                                                                                                                                                                 cGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCCCGCATTCGGCCAGATCGAT‑AA‑CGCGCCGGCAGAAAAAGCTCGTGGGATC‑ACCCACAACACCTCTCACGTTGAATACGACACCCCGCCCCGTCACTACGCACACGTAGAC‑CGCCCCGGGGC‑CGCCGACTATGTTTAAAACATGAACCCCCGGGCTGCGCCGCTGGACGGCGCGGTCCTGGTAGGTGCCGCCGCTGACGGCCCgc                                                                                    >  2:259953/1‑250 (MQ=11)                                                                                                                                                                     gCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTTTC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGc                                                                                 <  1:524204/250‑1 (MQ=25)                                                                                                                                                                                cGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTg                                                                     <  2:404664/251‑1 (MQ=33)                                                                                                                                                                                   aCTGGCTAAAACCTACGGCGGTGCTGCTCGCGCAGTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCGCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCGGCGACTGACGGCCCGATGCCGCAGACTCGTGAGc                                                                  <  2:87496/251‑1 (MQ=18)                                                                                                                                                                                      ggCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGGAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGca                                                                 <  2:93176/249‑1 (MQ=18)                                                                                                                                                                                          aaaaCCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATcc                                                            <  1:2640/250‑1 (MQ=33)                                                                                                                                                                                             aCCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGGTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTTCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCctgct                                                        <  2:449063/251‑1 (MQ=18)                                                                                                                                                                                                  cggcggTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtc                                                   <  1:131889/251‑1 (MQ=33)                                                                                                                                                                                                   ggcggTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTgggccg                                                  >  1:18311/1‑248 (MQ=18)                                                                                                                                                                                                      ggTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGgtcgtc                                                >  1:188330/1‑250 (MQ=33)                                                                                                                                                                                                            gctCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAgg                                         >  2:147082/1‑251 (MQ=33)                                                                                                                                                                                                                 cgcATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTc                                    <  2:293058/251‑1 (MQ=31)                                                                                                                                                                                                                  gcATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTcc                                   <  1:320888/251‑1 (MQ=31)                                                                                                                                                                                                                    aTTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                                 >  1:43973/1‑251 (MQ=31)                                                                                                                                                                                                                    aTTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGt                                 <  1:225152/251‑1 (MQ=31)                                                                                                                                                                                                                      tCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTGCCGTa                                <  1:167973/250‑1 (MQ=14)                                                                                                                                                                                                                          ccAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTAcatc                            >  2:45005/1‑250 (MQ=31)                                                                                                                                                                                                                              aTCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCCTCgt                       >  2:333524/1‑251 (MQ=14)                                                                                                                                                                                                                              aTCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCgt                       <  2:182105/251‑1 (MQ=31)                                                                                                                                                                                                                               tCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCgtg                      >  1:26575/1‑251 (MQ=31)                                                                                                                                                                                                                                 gAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGtt                    >  1:43796/1‑251 (MQ=31)                                                                                                                                                                                                                                 gAT‑AA‑CGCGCCGGAAGAAAAACCTCGTGGTATC‑ACCCTCCACACTTCACCCGTTGAATACGACACACCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCAGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGACAGGGAGGCGTTCCGTACATCATCGTGtt                    >  2:266932/1‑251 (MQ=11)                                                                                                                                                                                                                                  aT‑AA‑CGCGACGGAAGAAAAAGCCCGGGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTc                   <  2:164222/251‑1 (MQ=14)                                                                                                                                                                                                                                   gttaatgccccGGAATAAAAAGCTCGTGGTATA‑ACCATCAACACTTCTCACGATGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑GGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCGGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCCCATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGtt                    <  2:302643/242‑1 (MQ=11)                                                                                                                                                                                                                                     aa‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACCCCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCt                 <  2:11479/251‑1 (MQ=14)                                                                                                                                                                                                                                     aa‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAg                                                                            >  1:523814/1‑192 (MQ=32)                                                                                                                                                                                                                                     aa‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTGGCTGCGACTGACGGCCCGATGCCGCAg                                                                            <  2:523814/192‑1 (MQ=17)                                                                                                                                                                                                                                            ccGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACaa           <  1:147082/251‑1 (MQ=31)                                                                                                                                                                                                                                               gaagaaAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAATg        <  1:151152/251‑1 (MQ=31)                                                                                                                                                                                                                                                     aaaGCTCGTGGTCTC‑ACCATCAACACTTCGCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGCC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATACTGGTAGTTGCTGCGCCTGACGGCACGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTGGTCAGGTAGGCGGTCCGTTCCTCATCGGGTTCCTGAACAAATGCGACAt  >  2:305226/1‑251 (MQ=11)                                                                                                                                                                                                                                                     aaaGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAATGCGACAt  >  2:6420/1‑251 (MQ=31)                                                                                                                                                                                                                                                     aaaGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAATGCGACAt  >  1:497523/1‑251 (MQ=31)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           TTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATCCGTGTCTTAGAGGGACAATCGATGTCTAAAGAAAAGTTTGAACGTACAAAACCGCACGTTAACGTCGGTACTATCGGCCACGTTGACCATGGTAAAACAACGCTGACCGCTGCAATCACTACCGTACTGGCTAAAACCTACGGCGGTGCTGCTCGCGCATTCGACCAGATCGAT‑AA‑CGCGCCGGAAGAAAAAGCTCGTGGTATC‑ACCATCAACACTTCTCACGTTGAATACGACACCCCGACCCGTCACTACGCACACGTAGAC‑TGCCC‑GGGGCACGCCGACTATGTTAAAAACATGATCACCGGTGCTGCGCAGATGGACGGCGCGATCCTGGTAGTTGCTGCGACTGACGGCCCGATGCCGCAGACTCGTGAGCACATCCTGCTGGGTCGTCAGGTAGGCGTTCCGTACATCATCGTGTTCCTGAACAAATGCGACAT  >  NC_000913/4175870‑4176362

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