Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,061,935	G→A	S46N (AGC→AAC)	torD →	trimethylamine‑N‑oxide reductase‑specific chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,061,935	0	G	A	91.8%	192.8 / NA	74	S46N (AGC→AAC)	torD	trimethylamine‑N‑oxide reductase‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (39/28);  minor base C (0/6);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 7.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGGCG  >  NC_000913/1061691‑1062179                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      ggtggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCCCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGaa                                                                                                                                                                                                                                                   <  1:757019/249‑1 (MQ=255)  gtggAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAAcc‑g                                                                                                                                                                                                                                               >  1:204647/1‑251 (MQ=255)  ggtgAAATTGAGAAGGTCAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAAcc‑g                                                                                                                                                                                                                                               <  2:769904/248‑1 (MQ=255)         ttGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCCGATGGCTGAATGG‑ATTTCGTTG‑CTGGAAAACGAACc‑gccgcacc                                                                                                                                                                                                                                        >  2:476809/1‑248 (MQ=255)                  aCAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGcggcgg                                                                                                                                                                                                                               >  1:460054/1‑251 (MQ=255)                    aaCGGAACAGTGGAGCAGGTGACGACGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAACAACGAACC‑GCCGCTCACTGCGGCGGTg                                                                                                                                                                                                                             >  1:755063/1‑251 (MQ=255)                             ggggAGAAGGTGAGGGCTTTTAGCGGCCCGGTGGGGATGGTGGGGCGGTGCGAATATGT‑TCCCGC‑GGCG‑AGGGGAAAATCAGGCCCAC‑GCTGACAGCAC‑AACAGATTGCCTTTTTTTCC‑GCC‑TGGCTAGCGCAGTTGTTCGCCCTTGAGTTGGACGATAACCAACTGACAAAAACCGCCAGTGCGCAGATGGCGGAATGT‑TTTTCGTTG‑CGGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTg                                                                                                                                                                                                                    <  1:683696/249‑1 (MQ=255)                                gggcatgTGAAGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAA‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGATAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGGACA‑CCCGCTCACCGCGGCGGGGAActggcgtgg                                                                                                                                                                                                                  >  2:130506/7‑242 (MQ=255)                                 aGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAAAAACTGACGAAAACCGCCAGTGCGCAGATGGCTGAAGGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAAAGAGCTGGaaa                                                                                                                                                                                                                <  1:563833/251‑1 (MQ=255)                                  ggaggggACGGCGTTTAAGGCCCCCGTGGAGGTGGTGGCGCAGGGCGAATATGT‑CCCCGC‑GTCGAAGGGAAAA‑TCATGCCCAG‑GCTGACAGCAC‑AACAGTTTGCCTGTGTTTAC‑GCC‑TGGCAAGCGCAGTTTTTCTCCCGTGGGCTGGACGAGGAACAACTGACGAAAATCCCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGaaaa                                                                                                                                                                                                               <  1:619730/245‑1 (MQ=255)                                   cAGGTGACGGCGGTTAACGGCCCCGTGGAGATGGTGGCGAGGTGCAAATATGT‑CCCCGC‑TGCCCAGGTGAAA‑TCATGACCAC‑TCTGACA‑CACCAACAGATTCCCTGT‑TTTACCGCC‑GGGCTGGCGCTTTTGTTCTCCCGGGGGCGGGACGATGAACAACTGAGGCAAAGCGCCAGTGCGCAGGTGGCTGAAGGG‑TTTCCTTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAAc                                                                                                                                                                                                              <  1:513968/251‑1 (MQ=255)                                       ggacgggtgttacgggcccgggggggaTGGTGGGGCAGTGCGAATATGT‑CCCCGC‑GTCGCAGGTGAAA‑GCATGCCCAC‑GCTGACAGCAC‑AACAGATTGCCGGTGTTGAC‑GCG‑TGGCTAGCGCAGTTGTTCTCCCGGGAGCTGGACGAGGACCAACTGCCGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTGG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTa                                                                                                                                                                                                          <  2:51911/226‑1 (MQ=255)                                               tttAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGGGCTGGAAAACCGTATTGCCAcc                                                                                                                                                                                                  >  1:372566/1‑250 (MQ=255)                                                   aCGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGGACCAGCTGGAAAACCGTATTGCCACGCTgg                                                                                                                                                                                              >  2:144715/1‑250 (MQ=255)                                                     ggCCCCGTGGAGGTGGGGGCGCAGTGCGAATATGT‑TCCCGC‑GCCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTTTTTTTAC‑GCT‑GGGCTAGCGCAGTTGTTCTCCCGTGGGTGGGACGATAAACAACTGACGAAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACa                                                                                                                                                                                            <  1:62600/251‑1 (MQ=255)                                                          cGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAAcc‑gccg                                                                                                                                                                                                                                            >  1:114067/1‑198 (MQ=255)                                                                 ggagggtgggccggggcgattttgtCCCGCCGTCGCGGGGGAAA‑TCATGCCCCC‑GCTGACAGCAA‑AAAAGGTTGCTTGTGTTACC‑GCC‑GGGCTGGGGCAGTTGTGCCCCCGGGAGATGGACGACGACCAACAGACAAAAACGGCCTGTCCGCGGATGGCGTAATG‑‑TTTTTGTTGGCTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGAc                                                                                                                                                                                   <  2:453358/226‑1 (MQ=255)                                                                  tggtggCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCATATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATg                                                                                                                                                                               >  1:503330/1‑251 (MQ=255)                                                                    gtggCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGACTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTATCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACGGGCTGGAAAAACGTATTTCCACGCTGACAGTACGGGACGATGcc                                                                                                                                                                             >  2:364155/1‑251 (MQ=255)                                                                     gggCGCAGGGCGAAT‑TGTGTCCCGC‑GTCGCAGGTGAAA‑TCAGGACCCC‑GCTGACAGCAC‑AACATAGTGCCTGTGTTTAC‑GCC‑TGGCTGGCGCAGTTGTTCCCCCGTGAGCTGGAGGATGAACAACTGACGCAAATCGCCGGGGCGCAGATGGTTGAATGG‑TTTTCGTTG‑CTAAAAAACAAGCC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGccc                                                                                                                                                                            <  1:95180/250‑1 (MQ=255)                                                                         gcAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTcc                                                                                                                                                                        >  1:422628/1‑250 (MQ=255)                                                                          gagggCGAATATGT‑CCCCGC‑GTCGGAGGGAAAA‑ACAT‑ACCACCGCGGAAAACAC‑AAAAGATTGCCTGGGTTT‑C‑GCCCGGGCTGGCGCTGTTGTCCCCCCGGGGGCGGGACGAGGAACAACTGACGAAAACCGCCAGGGCGCAGAGGGCGAAATGG‑TTTGCGTTG‑CTAAAAACCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTg                                                                                                                                                                       <  2:422628/247‑1 (MQ=255)                                                                                  atatGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGc                                                                                                                                                               >  1:791850/1‑251 (MQ=255)                                                                                    atGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGCAAACCCTATTGCCCCACAGCCAGTACGTGACGATGCCCGTCTGGAACCGGccc                                                                                                                                                             >  2:437385/1‑250 (MQ=255)                                                                                              c‑gTCGCGGGTGAAA‑TCATTCCCGC‑GCTGCCAGCAC‑AATAGATTGCTTGTGTTTCC‑GCC‑TGGCACGCGCTGTTGTTCTCCGGGGAGCTGGACGAGGAACAACTGACGAAAACCGCCAGGGCGCAGGTGGCTGAAGGG‑TTTTGGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACtttt                                                                                                                                                    <  2:372566/251‑1 (MQ=255)                                                                                              c‑gTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGGACGTGCTGGAAAACCGGATTGCCCCGCTTAACGTACGTGAAGCGGCCCGTCTTGAACCGGGCGCGGGCtttt                                                                                                                                                    >  2:669449/1‑251 (MQ=255)                                                                                                gTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTg                                                                                                                                                   >  1:42561/1‑251 (MQ=255)                                                                                                                           tGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGacaa                                                                                                                           >  1:828951/1‑250 (MQ=255)                                                                                                                            gACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGGTTCTGATGGCCGacaaa                                                                                                                          >  1:43089/1‑250 (MQ=255)                                                                                                                            gACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGGGGACGGGGGGGAAAACAGGCTTGCCACGCTGACCGTACGAGAAGAAGACCGTCTGGAACTGGCCGCGGGCTTTTTCGACCTGTTTCtgattaactgcaaa                                                                                                                          >  2:737405/1‑240 (MQ=255)                                                                                                                                gCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGTCAAACAAgc                                                                                                                     >  1:659974/1‑251 (MQ=255)                                                                                                                                agcacaaaagttggctggtgtTTAC‑GCC‑GGGCTAGCGCAGTGTTTCTCCCGTGAGTGGGACGATGAACAACTGCCGCAAATCGCCGGGGCGCAGAGGGCGGAAGGT‑TTTTCGTGG‑CAGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAg                                                                                                                      <  1:437385/234‑1 (MQ=255)                                                                                                                                 cac‑aACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAgcg                                                                                                                    >  1:192695/1‑251 (MQ=255)                                                                                                                                 cac‑aACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGATGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGGGAACGAGCTGGCAAACCCGATTGCCACCCTGACCGTACGTGCCGCAGCCCGTCTGGACCCGGCCGCGGACTTTTTTGGCCTGtttttgtataccgccaaacacacc                                                                                                                    >  2:351873/1‑230 (MQ=255)                                                                                                                                 cac‑aACAGATTGCATGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGACCC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGCAAACCGTATTGCCACGCTGACAGTACGTGAAGATGCCCGACTGGAACTGGCCGCGGGCTTTTGCGGCCTGTTTCTGATGACCGACAAACATgcg                                                                                                                    >  1:451029/1‑251 (MQ=255)                                                                                                                                  caaaacaGATTGCTTGTGTTTAC‑GCC‑GGGCTAGCGCATTTTTTCCCCCGTGGGCTGGACGATGAACAACTGACGCAAACCGCAAGTGCGAAGATGGCTGAATGG‑TTTTCGTTG‑CTAAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGGAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAgcg                                                                                                                    <  1:408022/248‑1 (MQ=255)                                                                                                                                  ac‑aacaGATTGCCTGTGTTTCC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGAGGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTGG‑CTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAgcgg                                                                                                                   <  1:36753/251‑1 (MQ=255)                                                                                                                                       caGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTTTTCCCCCGTGAGCGGGACGATAAACAACTGAGGCAAATCGCCAGTGCGCAGAGGGCTGAATGG‑TTTTCGTTG‑CTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCt                                                                                                               <  1:197779/251‑1 (MQ=255)                                                                                                                                       agtgttgtgtgtgTTTAC‑GCT‑GGGATAGCGAAGTTGTCCTCCGGGGAGCTGGGGGATGAAAAAGTGACACAAATCGCAGGGGGGGAGATGGAGGATGGG‑TTTGCTTGG‑CAGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCt                                                                                                               <  2:175726/242‑1 (MQ=255)                                                                                                                                        aGATTGCTTGTGTTTGC‑GCC‑GGGCTACAGTAGTGGTTCTCCGTTGATGTGGCAGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGTTGAATGT‑TTTTGGTTG‑CTGAAAACCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTg                                                                                                              <  1:638000/251‑1 (MQ=255)                                                                                                                                         gATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGCTGACCGACAAAACAGCGGCGCTGc                                                                                                             >  1:797964/1‑251 (MQ=255)                                                                                                                                         gATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGc                                                                                                             >  1:778142/1‑251 (MQ=255)                                                                                                                                           ttcattgtgTTTAC‑GCC‑TGTCTAGCGCGGGTGGTCACCGGTGAGTGGGACGAGGAAAAACTGACGAAAACCGCCAGGGCGAAG‑TGGGTGGATGGGTTTTCGTGG‑TAGAAAAACAAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACTATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGcc                                                                                                            <  1:105336/245‑1 (MQ=255)                                                                                                                                              ccTGTGTTCAC‑CCC‑TGGCTAGCGCAGTTTTTCTCCCGGGAGCGGGACGAGGAAAAACTGACGCAAATCGCCAGTGCGAAGAGGGCTGAATGG‑TTTTCGTTG‑CTGAAAACCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTAt                                                                                                        <  1:805213/251‑1 (MQ=255)                                                                                                                                                    tttcgcccc‑TGGCGAGCGCATTTTTTCCCCCGGGGGCGGGACGATAAAAAACTGACACAAACCGCCAGTGCGCGGAGGGCTGAATGG‑TTTTCGTTG‑CTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCAt                                                                                                    <  2:265336/243‑1 (MQ=255)                                                                                                                                                    tttAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCAt                                                                                                    >  1:265336/1‑249 (MQ=255)                                                                                                                                                    tttAC‑GCC‑GGGCTAGCGCATTTGTTCCCCCGTGGGTTGGACGATGAACAACTGACGAAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATc                                                                                                   <  1:476809/250‑1 (MQ=255)                                                                                                                                                        c‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGGAAACCAGATTTCCCCGCCGACAGTACGTGACGATGCCCGTCTTGAACTCGCCGCGGACTTTTTCGGCCCGGTTCTTATGACCGCCAAACAAGCGGGGCCGCCGTATTCATCGGccc                                                                                              >  2:44276/1‑250 (MQ=255)                                                                                                                                                          gCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGCACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACCGTACGTGGCGATGCCCGTCTGGAACTGGCCGCGGGCCTTTGCGGCCTGGTTTTGATGACCGCCAAACCAGACGCGCTGCCGTATGCCTCGGcccc                                                                                             >  2:334145/1‑249 (MQ=255)                                                                                                                                                              tGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAAACGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACa                                                                                           >  1:340161/1‑250 (MQ=255)                                                                                                                                                               ggCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTAAGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGAACGACAAACAAGCGGCGCTGCCGTATGCATCCGCCTACaa                                                                                          >  1:697246/1‑250 (MQ=255)                                                                                                                                                                            gttCCCCCGTGAGCGGGAGGAAGAACAACTGACGCAAAGCCCCAGTGCGCAGAGGGCTAAAGGG‑TTTTCTTTG‑CTAAAAAACTAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAag                                                                            <  1:591513/251‑1 (MQ=255)                                                                                                                                                                                   cGTGATCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGTTGAATGT‑TTTTCTTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTGTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTaa                                                                     <  1:819882/251‑1 (MQ=255)                                                                                                                                                                                   cGTGAGCTGGACGATGAACAACTGACGCCAAACGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGGTTCTGATGACCGACAAACAAGCGGCGCTGCCGCATGCAACGGCCTAAAACCAGGACGAGCAAGAGAttt                                                                      >  1:735999/1‑249 (MQ=255)                                                                                                                                                                                        gCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGGCTTTTGCGGCCTGTTTCTGATGACCGCCAAACAAGCGGCGCTGCCGTTTGCCTCGGCCTCCAAACAGGGCGAGCAAGAGATTAAACGCt                                                                >  2:194309/1‑251 (MQ=255)                                                                                                                                                                                         cTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACAAGCTGGAAAACCGTATTGCCACGCCGACACGACGAGACGCTGCCCGTCAGGAACAGGCCGCGGGCTTTTGCCGCCTGTTTATGGTGGCCGTCCAACAAAAGGGGCTGCCGGATGCATCCGCCTTCCCACAGGACGAggaaacgaataaccacct                                                               >  2:449298/1‑234 (MQ=255)                                                                                                                                                                                             aCGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTTCCCCGCTGACAGGACGTGACGATGCCCGGCTGGAACTGGCCCCGGACTTTTGCGGGCTGGTTCTGGTGACCGACAAAAAAGCGGCGCTGCCGTATGCAACGGCCTACAAACAGGACGAGCAAGGGGTTAAACggctgttt                                                           >  2:149741/1‑244 (MQ=255)                                                                                                                                                                                             aCGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTACGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTTCCACGCTGACAGTACGTGACGATGCCCGTCTGGGACTGGCCGCGGACTTTTGCGGCCTGTTTTTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGCCCTCCAAACAGGACGAGCAAGGGATTAAACGCCTGtt                                                            >  2:371856/1‑250 (MQ=255)                                                                                                                                                                                                             gcaaatcgctgttGCGCATAGGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATgg                                           <  1:130506/239‑1 (MQ=255)                                                                                                                                                                                                                     ccAGTGCGCAGATGGCTGAGTGT‑TTTTCGTTG‑CAAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAga                                   <  1:44276/251‑2 (MQ=255)                                                                                                                                                                                                                      gggTGCGCGGTTGGGTGAATGG‑TTTTGGTTG‑‑AGAAAAACCAACCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCg                                  <  2:429257/249‑1 (MQ=255)                                                                                                                                                                                                                       agggcgcAGAGGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCgg                                 <  1:82813/248‑1 (MQ=255)                                                                                                                                                                                                                        gTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACA‑GCCGCTCACTGCGGCGGTGAACGGGCTGGAAAAACGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGACGTATGCACCGGCCCACAAACAGGACGAGCAAGGGATTAAAAACTTGTTAGTTGAGGCAGGGATGGCAACCAGAGGc                                >  1:199925/1‑251 (MQ=255)                                                                                                                                                                                                                         gtcgcAGATGGCTGATTGT‑TTTTGTTTG‑CTAAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCa                               <  1:632343/249‑1 (MQ=255)                                                                                                                                                                                                                           ccaggaggggggaagtg‑TTTGCGTTG‑CAAAAAACCGAACC‑GCCGCTCTCTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGAGGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACTAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAAt                             <  2:791850/235‑1 (MQ=255)                                                                                                                                                                                                                               gATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATAGGCCTACAAACAAGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTc                          >  1:725119/1‑250 (MQ=255)                                                                                                                                                                                                                                ggtgacttagtgtttttgtttG‑CTGAAAACAGAACC‑GCCTCTCATTGCGGCGGTGAACGAGCTCGAAAACCGTAGTGCCACTATTACAGTACGTGACGATGCCCGTCTGGAACCGGCCGCGGACTTTTGCTGCCTTGTTCTGATGACCGACAATCAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATCAAACGCTTGTTAGTTGAGGCAGGGATGGAACCCAGCGGCAATTTCa                         <  2:636687/232‑1 (MQ=255)                                                                                                                                                                                                                                        gtctcttccgttcatcaaaaaaACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAAcc                  <  1:817671/233‑1 (MQ=255)                                                                                                                                                                                                                                        aTGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCCAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACTAACCgg                >  1:407181/1‑251 (MQ=255)                                                                                                                                                                                                                                         tGG‑TTTTCGTTG‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGAAGGTGCCCGTGTGGAACTTGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGAGGCGGTGCCGTAAGCATCGGCCTACAAACAGGAAGAGGAAGGGATTAAAAGCTTGTTAGTTTAGGCAGGGATGGAAAACAACGGGAATTTCCACGAACCGGc               >  1:605120/1‑251 (MQ=255)                                                                                                                                                                                                                                                    ag‑CTGAAAAACGTACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCt      <  2:204647/249‑1 (MQ=255)                                                                                                                                                                                                                                                     g‑CTGAAAAACGAACC‑GCCGCTCACTGCGGCGGGGAACGAGATGGAAAAACGGATTTCCACCCTGCAAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGGCTTTTGCGGCCTGTTTCTGATAACCGACAAACAAGCGGCGCTGCCGGATGCATCGGCCTACAAACAGGACGGGCAAAAAATTCAACGCTTGTTAGTTTAGGCAGGGTTGGGAACCAGCGGCAAtttataagaaccaaatcgggccagcg    >  2:817671/1‑228 (MQ=255)                                                                                                                                                                                                                                                        tGAAAAACGAACC‑GCCGCTCACTGCGTCGGTGAACGAGCTGGAAAACCGTATTGCCGCGCTGACAGTACGTGACGTTGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCCTCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGAGGGAAAGCAGCGGCACTTTCAAAGAACCGGCAGACCATCTGGCg  >  1:516449/1‑251 (MQ=255)                                                                                                                                                                                                                                                        aaaaaaCCGACCC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGGCg  <  2:503330/249‑1 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                      GGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGT‑TCCCGC‑GTCGCAGGTGAAA‑TCATGACCAC‑GCTGACAGCAC‑AACAGATTGCCTGTGTTTAC‑GCC‑TGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGG‑TTTTCGTTG‑CTGAAAAGCGAACC‑GCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGGCG  >  NC_000913/1061691‑1062179

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