Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	651,297	0	A	G	26.4%	148.9 / 11.1	72	F187S (TTC→TCC)	citC	citrate lyase synthetase
Reads supporting (aligned to +/- strand):  ref base A (15/38);  new base G (19/0);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 1.58e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACA  >  NC_000913/651054‑651546                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                gTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAAGTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTga                                                                                                                                                                                                                                                          <  2:289048/251‑1 (MQ=255)         aCAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAAtctt                                                                                                                                                                                                                                                  <  1:595091/251‑1 (MQ=255)             aaGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATTATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGTGAAGCGTGAAGAATCTTCttt                                                                                                                                                                                                                                              <  2:408302/251‑1 (MQ=255)             aaGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCttt                                                                                                                                                                                                                                              <  2:227616/251‑1 (MQ=255)                ggTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGAc                                                                                                                                                                                                                                           <  2:529364/251‑1 (MQ=255)                 gTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACt                                                                                                                                                                                                                                          <  1:483578/251‑1 (MQ=255)                    cAGTACCGGCAAATCGGTGAGTTCCACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAAGTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTAAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATGCGGAGCCACGATGCACAGTCAGTCTGGGAATAGCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTaaa                                                                                                                                                                                                                                       <  2:609083/251‑1 (MQ=255)                         ccGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGCGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCCGCCGGTCTTTAT‑A‑GGGGAAGCGTGAAGAATCTTTTTTTGCTaaaaaaa                                                                                                                                                                                                                                   >  1:553202/1‑248 (MQ=255)                               aaGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACA‑CC‑AAATCCAGCCGGGCTTCAT‑A‑GGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAAATGCa                                                                                                                                                                                                                             >  1:530153/1‑250 (MQ=255)                                    gTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTGTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGGGGAATATCGGCGGTGGCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAG‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGt                                                                                                                                                                                                                      <  1:421553/251‑1 (MQ=255)                                         ttACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACA‑CC‑AAATCCAGCCGGGTTTCAT‑A‑AGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAAATGCCACC‑AGTCGCAc                                                                                                                                                                                                                 >  1:164995/1‑251 (MQ=255)                                               ccAGCGCGGGAGCGAGGTAATGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACC‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑cc                                                                                                                                                                                                          <  2:585754/251‑1 (MQ=255)                                                    gctggAGCTAGGTACTGACGGAAAATCTGCAGATCAATGTCGGTGTAACAAGGGTTAATGACGCTCGGTTCTTTAATGAAGTAGCAAGGGAACGTAGCTCTGGAGAGTATGTATTCGGAGCCACGATGCACAGTCAGGCGTGTAATATCGGCGGTGCCTGTTAACA‑CC‑AAATCCAGCCGGACGGCAG‑A‑GGTGAAACGTGAACAATCTTCttt                                                                                                                                                                                                                                              <  2:668811/209‑1 (MQ=255)                                                      gggAGCGAGGTCCTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTCATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTAACA‑CC‑AAAACCCGCCGGGCTTC‑T‑ATGGGGGAGCGGGAAGAATCTTCTTTTACTACAAACAGATGCAACC‑CGTCGCCCTGTG‑CCGCAGccc                                                                                                                                                                                                   >  2:663425/1‑250 (MQ=255)                                                          gCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTg                                                                                                                                                                                               <  2:330644/251‑1 (MQ=255)                                                          gCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTg                                                                                                                                                                                               <  1:544461/251‑1 (MQ=255)                                                               gTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCt                                                                                                                                                                                                   <  1:164489/242‑1 (MQ=255)                                                               gTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTTACA‑CC‑AAAACCAGCCGGTCTTCAT‑A‑GGGGGAGCGGGGAGAATCTTTTTTGTCTAAAAACAGATGCAACC‑AGTCGCACTGTG‑GCGCCGCCtatatcgga                                                                                                                                                                                           >  2:164489/1‑242 (MQ=255)                                                                  cTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAACA‑CC‑CAATCCAGCCGGTCTTTCT‑A‑GGGGGAGCGTGCAGAATCTTTTTTTAATAAAAACAGATGCAACC‑ACTCGCACGGTG‑GCGCCGCCTGTTTAATCAggt                                                                                                                                                                                       >  1:381575/1‑249 (MQ=255)                                                                                  aGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGt                                                                                                                                                                       <  2:627190/251‑1 (MQ=255)                                                                                    aTCAATTTCGGTGTAACAATGGTTAATGGCGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGGCAGGCGTGGAATATCGGCGGTGCCTTTTAACC‑CC‑AAATCCAGCCGGGCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAAGCAGATAACGGTGACCATTCGTaa                                                                                                                                                                     <  2:272886/251‑1 (MQ=255)                                                                                    aTCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAGCC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATGACGGTTACCATTCGTaa                                                                                                                                                                     <  1:504789/251‑1 (MQ=255)                                                                                           tCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑t                                                                                                                                                             <  2:278825/250‑1 (MQ=255)                                                                                           tCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑t                                                                                                                                                             <  2:239962/250‑1 (MQ=255)                                                                                                  aaCAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTtc                                                                                                                                                     <  2:530153/251‑1 (MQ=255)                                                                                                     aaTGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAA‑CC‑CAATCCAGCCGGGCTTCAT‑A‑GGGGGAGAGTGAAGAATCTTCTTTGACTAAAAAAAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGGCCATTCGTAAA‑GGGA‑TTgggcgatcgg                                                                                                                                                  >  1:271376/1‑243 (MQ=255)                                                                                                     aaTGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGc                                                                                                                                                         <  2:271376/244‑1 (MQ=255)                                                                                                         gTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAA‑CC‑AAAACCAGCCGGGCTTCAT‑T‑GGGGGAGCGGGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCCCCGTG‑CCGCAGCCCGTTGAATTAGAAAACGGTGACCATTCTTAAA‑AGGA‑TTGGCGTtcatccca                                                                                                                                               >  1:311540/1‑247 (MQ=255)                                                                                                          gtAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑t                                                                                                                                           <  1:505451/250‑1 (MQ=255)                                                                                                            aaTGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACC‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTGGGAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGc                                                                                                                                         <  1:550689/251‑1 (MQ=255)                                                                                                                aCGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTTACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCCGCCTGTTGAATCAGATAACGGGGACCATTCGTAAA‑AGGA‑TTGGCGTTCATCAC‑AA‑TGCAGcc                                                                                                                                     >  1:191246/1‑251 (MQ=255)                                                                                                                    tctGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGGGCCTTTTTACA‑CC‑AAAACCAGCCGGTTTTCCT‑A‑GGGGAAGGGTGAAAAAACTTCTTTTACTCAAAACAAATGCAACC‑CGTCGCACTGTG‑GCGCAGCCCGGTGAAACAGATAACCGCGACCATTCGGAAA‑GGGA‑TTGGCGTTTATCAC‑AA‑TGAAGCCAATc                                                                                                                                 >  1:149117/1‑251 (MQ=255)                                                                                                                              aaTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGATCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCGTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAAGCAGATAACGGTGACCATTCTTAAA‑GGGA‑GTGGCGTTCATCAC‑AA‑AGCAGCGAATCTTGGTCCCTg                                                                                                                       <  2:631134/251‑1 (MQ=255)                                                                                                                              aaTGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAA‑CC‑AAATCCAGCCGGCCTTCAT‑A‑GGGGGAGCGGGAAGAATATTTTTTGACTAAAAAAAGATGCAAAC‑ACACGCCCCGTG‑CCGCAGCCTGTTGAATCAGATAACAGGGACCATTCGTTAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TTCAGCCAATATTGTTCCCTg                                                                                                                       >  1:481735/1‑251 (MQ=255)                                                                                                                                  aaGTAGCAAGGGACCGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCTGGGCTTCATCA‑GGGGAA‑CGTGAAGAATCTTTTTTGACTAAAAACACATGCAACCCAG‑CGCACTGTG‑GCGCAGCCTGTTGCAACAGATAAAGGGGGCCATTCGCAAA‑GGGC‑TTGGCGGTCATCAC‑AA‑TGCCGCCAATCTTTGTCCCTGGagt                                                                                                                   >  1:442158/1‑249 (MQ=255)                                                                                                                                  aaGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACA‑CC‑AAAACCAGCCGGGCTTCAT‑T‑GGGGGAGCGGGGAGAATCTTCTTTTACTAAAAACCGATGCAAAC‑AGTCGCCCTGTG‑CCGCAGCCCGTTGAATCAAATAACGGGGACCATTCGTTAA‑GGGG‑TTGGGGCTCATCCCCAA‑T‑CAGCCAAACTTGTTTCCttggtg                                                                                                                   >  2:170228/1‑246 (MQ=255)                                                                                                                                  aaGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCATTTTAACA‑CCCAATTCCACCCGG‑CTTTCT‑A‑GGGGGAGCGTGAAGAATTTTCTTTGACTAAAAAAAGATGCAACC‑AGTAGCAATGTGGCCGC‑GCCCGTTGAATCAGATAACGGGGACCATTCGTAAA‑AGGAGTTGG‑GTTCAT‑AC‑AAATGCAGCCCAACTTGTTTCCTgggtg                                                                                                                   >  1:275659/1‑248 (MQ=255)                                                                                                                                       gCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTCTGACTAAAAAAAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAAAGGGGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGa                                                                                                               >  1:192948/1‑250 (MQ=255)                                                                                                                                         aaGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAAt                                                                                                            <  1:555967/251‑1 (MQ=255)                                                                                                                                              aaCGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTGGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGGTCCCTGGATTACGAAAtttttt                                                                                                       <  2:625422/251‑1 (MQ=255)                                                                                                                                                     ggcgGGAGATGAGGTATGCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGAt                                                                                                <  2:218526/250‑1 (MQ=255)                                                                                                                                                       cgGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACA‑CC‑AAATCCAGACGGGCTTCAT‑A‑GGGGGAGCGGGAAGAATCTTCTTTGACTAAAAAAAGATGCAACC‑AGTCCCCCTGTG‑CCGCCGCCTGTTGAATCAGATAAAGGGGACCATTCGTAAA‑GGGG‑ATGGCGGTCCTCCC‑AA‑TGCAGCCCATCTTGTTCCCTGGGTGAAGAAAtttttttcgcggttt                                                                                              >  2:97857/1‑242 (MQ=255)                                                                                                                                                           agaTGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACA‑CC‑AAAATCAGCCGGTCTTCAT‑A‑GGGGGAGCGTGAAGAATCTTCTTTGACTAAAAAAAGATGCAACC‑AGTCGCACTGTG‑GCGCAGCCTGTTGAATCAGATAAAGGTGACCATTTGTAAA‑AGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAAACTTGTTTCCTGGATGACGAAAATTTTTCAGCGAATCGGc                                                                                           >  1:366445/1‑250 (MQ=255)                                                                                                                                                            gatgatGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTAAA‑CC‑AACTACAGGCGGTCTTCAT‑A‑GGGGGAGGGAGAAGAATATTTTTTGCCTAAAAAAACATGCAACC‑AGTTGCACTGTG‑CCGCCGCCTGTTGAATCAAGTAAAGGTGACCCTTCGTAAA‑GGGA‑ATGGCGTTCATTAC‑AA‑AGCAACCAATATTGTTCACTGGATGACGAAATTTTTTTAGAGATTTGGc                                                                                           >  1:114585/1‑249 (MQ=255)                                                                                                                                                             atgatgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAACAACC‑AATTCC‑CCCGGCCCTACT‑A‑AGGGGAGCCTGGAAAATCTTTTTTTACTAAAAAAAGAAGCAACC‑AGGCGCACTGTG‑CCGCCCCCTGTTGAATCAGAAAACGGTGGCCATTCGTAAAAGGGA‑T‑GGCGTTTCTCAC‑AA‑TGCAGCCCATCTTGTTTCCCGGAAGACGCAATTTTTTTAGCGCTTCGGCAtt                                                                                        >  1:85600/1‑250 (MQ=255)                                                                                                                                                              tgatgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGGCTTTTAACA‑CC‑CCAACAAACCCGTCTTCAT‑A‑GGGGGAGCGGGGAGAATCTTCTTTTAATAAACACAGATGCCACC‑AGTCGCACTGTG‑CCGCCGCCTGTTGAATCAGATAACGGTGAACATTCGTAAA‑GGGA‑TTGGCGGTCATCAC‑AA‑TGCAGCCAAACTTGTTTCCTCGGTGGCGGAATTTTTTTAGCGAtttgtcata                                                                                        >  2:397505/1‑243 (MQ=255)                                                                                                                                                              tgatgTATTCGGAGCAACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACA‑CC‑CAATCCAGCCGGGCCTCCT‑A‑GGGGGAGCGGGAAGAAGCTTCTTTGACCAACAACAGCAGCAACC‑AGACGCACTGTG‑CCGCCGCCTGTTGAATCAGATAACGGGGATCATTCCTAAA‑GGGG‑TTGGCGGTCATCAC‑AA‑TGCAGCCAATCCTTTTCCCTGGATGCCGAAATTTTTTTAGCGATTTGGcct                                                                                         >  1:499938/1‑247 (MQ=255)                                                                                                                                                                atgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATAACGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTGGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TGGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTGCCCTGTATGACGAAATTTTTTCAGCGATTCGGCATAGCg                                                                                     <  2:669959/251‑1 (MQ=255)                                                                                                                                                                 tgTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTAACA‑CC‑AAAACCAGCCGGTCTTCCT‑A‑GGGGGAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AATCGCCCTGTG‑CCGCAGCCTGTTGAATCAGATAACGGGGACCCTTCGTAAA‑AGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTTTTCCCTGGATGACGAAATTTTTTTCGCGATTCGGCATAACGt                                                                                    >  1:475861/1‑251 (MQ=255)                                                                                                                                                                    aTTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGttt                                                                                  <  1:199609/250‑1 (MQ=255)                                                                                                                                                                           gCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACA‑CC‑AAATCCAGCCCGGCTTTCT‑A‑GGGGGAGCGGGAAGAATCTTCTTTGTCTAAAAAAAGATGAAAAC‑AGGCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGGGACCATTCGTAAA‑GGGG‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTTGATGACGAAATTTTTTTAGCGATTCGGCATATCGGTTCAATcgcg                                                                          >  1:501165/1‑251 (MQ=255)                                                                                                                                                                           gCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCCTTTTACA‑CC‑AAAACCAACCGGGCGTCCT‑A‑TGGGGAGGGGGGAGAATCTTTTTTTACTTAAAAAAAAAGCAAAC‑AAGCGCACTTTG‑CCGCAGCCTGTTGAAACAGAAAAAAGGGACCATTCGTAAA‑GGGG‑TTGGGGTTTCTAAC‑CA‑AGCAGCCCAACTTTTTCCCCTGATGAAGGAATTTTTTTTGCGATTCTGCATAGCggttttatcgcg                                                                          >  2:153886/1‑240 (MQ=255)                                                                                                                                                                                       aGTCAGGCGTGGAATATCGGCGGTGCCTTTTTACC‑CC‑AAATCCAGCCGGTCTTCAT‑A‑AGGGAAGCGTGAAGAATCTTCTTTGGCTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGAAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGG‑ATGGCGTTCATCAC‑AA‑AGCAGCCAATCTTGTTCCCTGGGTGACGAAATTTTTTTAGCGGTTTGGCCTTGCGTTTCAGTCGCGGGGGGCTGtttt                                                              >  1:572009/1‑251 (MQ=255)                                                                                                                                                                                       aGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGtttt                                                              <  1:97857/251‑1 (MQ=255)                                                                                                                                                                                          cAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGTCGTTCAGCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACTAAATTTTTTCAGCGAGTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTcc                                                            <  1:411175/250‑1 (MQ=255)                                                                                                                                                                                                        cggcggTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACg                                             <  1:663425/251‑1 (MQ=255)                                                                                                                                                                                                                ccTTTTAACG‑CC‑AAACCCAGCCGGTCTTCAT‑A‑GTGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGGTGAAGCAGATAACGGTAACCATGCGTAAA‑GGGA‑TTGGCGTTCAGCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCTCGTGGCGCTGTTTTCCAGCAGCTCCATCACGCCGGGTAc                                     <  2:443840/251‑1 (MQ=255)                                                                                                                                                                                                                  ttttacaa‑cC‑AAATCCAGCCGGTCTTCAT‑C‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCCACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGCTAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATCGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGTGTACGc                                   >  2:197673/8‑251 (MQ=255)                                                                                                                                                                                                                     tAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTgg                                >  2:128714/1‑251 (MQ=255)                                                                                                                                                                                                                          ccc‑AAAGCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTGGCACTGTG‑CCGGAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAg                            <  2:8938/250‑1 (MQ=255)                                                                                                                                                                                                                              aaaTCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGCTGAATCTTCTTTGACCAAAAACAGATGCAACC‑AGACGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGAGACCATTCGTAAT‑GGGC‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTACCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCCCGCCGggtccgcgggtccgcgg                         >  2:261047/1‑250 (MQ=255)                                                                                                                                                                                                                               aaTCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTg                        >  1:128955/1‑251 (MQ=255)                                                                                                                                                                                                                                    aGCCGGTCTTCAT‑A‑GGGTAAGCGTGAAGAATCTTCTTTGACTAACAACAGATACAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAACCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCAGCATAGCGTTTCAGTAGCATGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGaac                    >  1:329882/1‑249 (MQ=255)                                                                                                                                                                                                                                    aGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAACCAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCCGCACCATCACGCCGGGTACGCTGGTCAGCGTGGaaa                    >  2:395720/1‑250 (MQ=255)                                                                                                                                                                                                                                       cGGTCCTCAT‑A‑GGGGAAGCTTGAAGAATCTGCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGAGAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTCTTTGCCATCAGCACCAGCACGCCGGGTACGCTGGGCAGCGTGGAAAAAcc                <  2:114585/251‑1 (MQ=255)                                                                                                                                                                                                                                        ggTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCg               >  2:514015/1‑251 (MQ=255)                                                                                                                                                                                                                                            ttCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑CGTCGCACTGTG‑CCGCAGCCAGTTGAAACAGCTAACGGTGACAATTCGAACA‑GGGA‑ATGGCGTTCATCAC‑AA‑TGCAGCCCATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACt           >  1:640868/1‑251 (MQ=255)                                                                                                                                                                                                                                            ttCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGGTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCAc            <  2:553202/250‑1 (MQ=255)                                                                                                                                                                                                                                             tCAT‑A‑GGGAAAGCGTGCAGAGTCTTCTTTGACTACAAAAAGATGCCACG‑AGTCGCACTGTG‑CCGCCGCCTGTTGAATCAGATCACGGTGACCATTCGTCAA‑GGGA‑TTGGCGTTCATCAC‑AC‑TGCAGCCAATCTTGTTCCCTGGATGACGAACTTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCCCGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTg          >  2:629092/1‑251 (MQ=255)                                                                                                                                                                                                                                                      ggAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGAGTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCTCTGGCGGAAc   <  2:314153/251‑1 (MQ=255)                                                                                                                                                                                                                                                       gAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACa  <  1:581641/251‑1 (MQ=255)                                                                                                                                                                                                                                                       gAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACa  <  1:333008/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                GTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTTAATGAAGTAGCAAGGGAACGTAGCGCGGGAGATGATGTATTCGGAGCCACGATGCACAGTCAGGCGTGGAATATCGGCGGTGCCTTTTAACA‑CC‑AAATCCAGCCGGTCTTCAT‑A‑GGGGAAGCGTGAAGAATCTTCTTTGACTAAAAACAGATGCAACC‑AGTCGCACTGTG‑CCGCAGCCTGTTGAATCAGATAACGGTGACCATTCGTAAA‑GGGA‑TTGGCGTTCATCAC‑AA‑TGCAGCCAATCTTGTTCCCTGGATGACGAAATTTTTTCAGCGATTCGGCATAGCGTTTCAGTCGCGTGGCGCTGTTTTCCATCAGCACCATCACGCCGGGTACGCTGGTCAGCGTGGAAAAACCGCACTGGCGGAACA  >  NC_000913/651054‑651546

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