Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,296,179	A→T	10.1%	I125F (ATT→TTT)	diaA →	DnaA initiator‑associating factor for replication initiation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,296,179	0	A	T	10.1%	227.0 / 4.1	89	I125F (ATT→TTT)	diaA	DnaA initiator‑associating factor for replication initiation
Reads supporting (aligned to +/- strand):  ref base A (60/18);  new base T (2/8);  total (62/27)
Fisher's exact test for biased strand distribution p-value = 6.99e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAATGGAACTTCCGCTGCCAATGCACAGCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACAT  >  NC_000913/3295949‑3296408                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            gTAATGGAACTTCCGCTGCCAATGCACAGCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTaaa                                                                                                                                                                                                                                     >  2:1280249/1‑239 (MQ=255)    aTGGAACTTCCGCTGCCAATGCACAGCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCa                                                                                                                                                                                                                                  >  2:1283393/1‑239 (MQ=255)         aCTTCCGCTGCCAATGCACAGCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACC‑GGTGCGGGCGCTGGGTCATTCGGGAG‑AGGAATTTTTAGCCATTTCCACCCAGGGCAACAACCGCGCAATTGTTAAATAAGttta                                                                                                                                                                                                                             >  1:281224/1‑237 (MQ=255)             ccGCTGCCAATGCACAGCATTTTGCTGCCAGCATGAACAACCGTTTCAAAACGGAGCGGCCCAGTTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTTTATGCAAAACA‑GGGGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAgcc                                                                                                                                                                                                                         <  2:618928/239‑1 (MQ=255)                       tGCACAGCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTAcgc                                                                                                                                                                                                               >  2:687825/1‑239 (MQ=255)                             gCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACGACCGCGATATTGTTAAAGCAGTTGAAGCCGCCGGTACGCGTGata                                                                                                                                                                                                         >  2:1152909/1‑239 (MQ=255)                             gCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATTTTGTTAAAGGAGTTGAAGCCGCCGTTAGCCGTGata                                                                                                                                                                                                         >  2:1176206/1‑239 (MQ=255)                                   tgctgcCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGGGTATGCAAAAAA‑GGTGCGGGCCCTGGGTCTTCCGGGAG‑ATGTATTGTTAGCCCTTTCCCCCCGTGGCAAAAGCGCCGATATTGTTAAAGCAGTTAAAGCCCCCGTTACGCGTGaaaggaaac                                                                                                                                                                                                   >  1:712507/1‑231 (MQ=255)                                        ccAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCCACCGTGGCACACGCCCGGAAAATGTTAAAGCAGTTGAACCCCCCCTTACGCGTGaataagccatattg                                                                                                                                                                                              >  1:1236884/1‑226 (MQ=255)                                        ccAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGACAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCAtgttg                                                                                                                                                                                              >  2:201455/1‑235 (MQ=255)                                         cAGCATGATCAACCGTTTAGAAACGAAACGGCCTAGCTTACCTGCCATTGCACTAACTACTGATACTGTTGTCTTAACGGCGATTGCCAACGATCGCTAACAAGATGAAGTGTATGCAAAAAA‑GGTGCGGGCGCTGGGTCATTCGGGAG‑AAGTATTGCTAGCCATTTCCACCCGTGGAAAAAGACGCGAGTTTGTTAAAGCAGTTGAAGCAGACGTTACGCGTGATATGACCATTGTgg                                                                                                                                                                                             >  2:627557/1‑239 (MQ=255)                                            cATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGGG‑ATGTATTGTTAGCCATTTCCACCCGTGGCCACAGCCGCGATATTGTTAAAGCAGTTGAAGCCCCCGTTACGCGTGCTATGGCCATTGTGGcct                                                                                                                                                                                          >  2:44177/1‑237 (MQ=255)                                                  cAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTTTTGTTAGCCATTTCCACCCGTGGCAAACACCGCGATATTGTGAAACAAGTTTAAACCCCCGTTAGCCCTGATATGACCATTAGTGCATGGACCg                                                                                                                                                                                    >  1:808818/1‑239 (MQ=255)                                                    aCCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGGGATTTGCCCATTGTGGCATTGACCg                                                                                                                                                                                    >  2:499186/1‑237 (MQ=255)                                                      cGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTCTGCAAAACC‑GGTGCGGGCGCAGGGTCATGCGGGCG‑ATTTCTTGTTAGCCATTTCCCCCCGTGGCAACAGCCGCGAGATCGTTAAAGCAGATGAAGCCCCCGTGACGCGTGATATGACCacttaggaattacacggcaa                                                                                                                                                                                >  1:231773/1‑220 (MQ=255)                                                        tttCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGAAATGACCATTGTGGCTTTGACCGGcattg                                                                                                                                                                              >  2:313874/1‑235 (MQ=255)                                                             aaaCGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGAAATGACCATTGTGGCATTGACCGGCTATGAcggc                                                                                                                                                                         >  2:349590/1‑239 (MQ=255)                                                                                cccTGCCATTGGACTAAATCCTGAAAATGTTGTCTTAACGGCGATTGCCAACGATGGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCCCCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGt                                                                                                                                                     <  2:281429/238‑1 (MQ=255)                                                                                  cTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGGGTGAAATCACCCTTGGGGCCTTGACCGGCTATGCATGCGGCGAACTT‑GCCGGTTTGTTa                                                                                                                                                   >  2:602112/1‑239 (MQ=255)                                                                                  cTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCACCCAGTGGCACAAGCCGCGATATTGGTAAAGCATTTGAAGCCGCCGTTACGCGTGATAAGACCATTGTGGCATTTACCGGCATGGCCGGCGCGGAACtttacgggttttta                                                                                                                                                    >  1:225113/1‑227 (MQ=255)                                                                                        ttGAACTAAATCCTAAAAGTGTGGGCTTAACGGCGTTTGCCAACGATCGTTAAAATGATGAGTGGTATG‑AAAACAGGGGGGGGGGGCGGGGTCAGGGGGGAG‑ATTTTTTGTTAGCAATTCCCCCCCGTGGCAACAGCCGCGATATTTTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCca                                                                                                                                             <  2:1290990/239‑1 (MQ=255)                                                                                                aaTACTGATAATGTTGTCTTAACGGCGATTGCCAACGCTCGCTTACATGATGAAGTGTATGCAACACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAATCAGTTGAAGCCGCCGTTACGCGTTATCTGCCCCTTGTGGCATTGACCGGCTATGACGCCGGCAGACAT‑ACAGGTTTGTTAGGGCcaccagaatg                                                                                                                                     >  1:574239/1‑232 (MQ=255)                                                                                                     tGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTTCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGGTATTGTTCAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTTACCAGTTTTGTTAGGCCCAAAggtgtgtggt                                                                                                                                 >  1:626960/1‑231 (MQ=255)                                                                                                      gATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGgaga                                                                                                                                >  2:505352/1‑238 (MQ=255)                                                                                                       aTAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGGGAt                                                                                                                               >  2:868867/1‑238 (MQ=255)                                                                                                       aTAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGACATCCGGGAG‑ATGTATTGTTAGACATTTCCACACGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCAGCCGTTACGCGTGACATGACCATTGATGCCATGACAGGCAAAGACGGGGGCGAAGTT‑GCCGGTTTGTTAAGGCCACAGGAGGGGGAGGTc                                                                                                                              >  1:179269/1‑239 (MQ=255)                                                                                                       aTAATGTTGTCTTAACGGCGATTGCCAACGATCGCATACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGGCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGCGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGCCCCCAGGCTGTGGAGAt                                                                                                                               >  2:112285/1‑238 (MQ=255)                                                                                                        tAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTAAATGATGAAGTGGATGCAAAACA‑GGGGCGGGCGCTGGGGCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCATACGACGAGCGATTGTTAAAGCAGTGTAAGCCGCCGTTACGCGTGATATTCCCATTTTGGCATTGACGGGCTAGGACGGCGGCGACATT‑GGCGGATTCGTGGGGCCACCGGATCTGGGAATcc                                                                                                                             >  1:40688/1‑239 (MQ=255)                                                                                                            gttgtCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTTGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGCGATATGTCACTTGTGGACTTGAACGGCTCTACAGGCCGCGGACTT‑GCGGGGTTTTTAggggccaagggagggggattccgcat                                                                                                                         >  1:418869/1‑216 (MQ=255)                                                                                                            gttgtCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCACGATTGTTAGGGCCACAGGATGTGGAGATCCGCAt                                                                                                                         >  2:77504/1‑239 (MQ=255)                                                                                                               gtCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTAcgc                                                                                                                                                                                                               >  2:386190/1‑151 (MQ=255)                                                                                                               gtCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATtgtg                                                                                                                                                                                              >  2:857790/1‑168 (MQ=255)                                                                                                               gtCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGGGATCCGCAtt                                                                                                                        >  2:845469/1‑237 (MQ=255)                                                                                                                     aCGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCa                                                                                                              >  2:293373/1‑239 (MQ=255)                                                                                                                          gATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTa                                                                                                         >  2:391488/1‑239 (MQ=255)                                                                                                                          gATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTa                                                                                                         >  2:546838/1‑239 (MQ=255)                                                                                                                          gATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTa                                                                                                         >  2:894378/1‑239 (MQ=255)                                                                                                                                    gATCGCTCAAATGATGAATTGTATAAAAACCG‑GGAGCGGGCGTTGGGTCATGGGGGAG‑ATTTATGTTTACCCTTTCCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTAcg                                                                                                                                                                                                                <  1:1235645/129‑1 (MQ=255)                                                                                                                                          ttACATGATGAAGTGTATGCAAAACG‑GGGGGGGGCGCTGGGTCATGGGGGAG‑ATGTATTGGTAGCCATTCCCCCCCGGGGCAACAGCCGGGATATTGTAAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGa                                                                                         <  2:281224/239‑1 (MQ=255)                                                                                                                                           tagttgaggaAGTGTAGGCAAAAGG‑GGGGTGGGTGGGGGGTGATG‑GGGAGTATTTATTGTATGCCATTTCCCCCGGTGGAACCCGCCGCGATTTTGATAAAGCGTTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGa                                                                                         <  1:894378/234‑1 (MQ=255)                                                                                                                                            aCATGATGAAGTGTAGGAAAAAGA‑GGTGGGGGGGGTGGGTCATGCGGGAG‑ATGTTTTGTTAGCCATTCCCACCCGGGGCAACAGCCGCGATTTTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGaaa                                                                                       <  2:712507/239‑1 (MQ=255)                                                                                                                                              atgatgAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGa                                                                                                                                                                                                      >  2:1235645/1‑129 (MQ=255)                                                                                                                                                    aaGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGGGAACTT‑GCAGGTTTGTTAGGGCCCCAGGGTGGTGAGATCCGCCTT‑‑CCTTCCCCACGTGATGGTCCCAATCAGGaaaagaaaatg                                                                               >  1:711660/1‑231 (MQ=255)                                                                                                                                                        gtATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGCGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTg                                                                            >  2:25506/1‑238 (MQ=255)                                                                                                                                                         taggaaaaaCG‑GGTGGGGGCGCTGGGTCATGCGGGAG‑TTGTATTGTTAGCATTTTCCCCCGGGGGCAACAGCCGCGATATTGTAAAGGCAGTTGAGGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGAc                                                                          <  2:7286/234‑1 (MQ=255)                                                                                                                                                          aTGCAAAAGAGGGGGCGGGGGCTGGTTCATGGGGGAG‑TTTTTTTGTTAGCCATTTCCCCCGGTGGCAACAGCCGCAATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGAc                                                                          <  2:1131102/239‑1 (MQ=255)                                                                                                                                                           tGCAAAACG‑GGTGGGGGGGCGGGGTCATGGGGGAG‑ATGTATTGTTAGCCTTTTCCCCCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGga                                                                                                                                  <  2:1265447/183‑1 (MQ=255)                                                                                                                                                            gCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGt                                                                       >  2:1218288/1‑239 (MQ=255)                                                                                                                                                                aaCA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑c                                                                                                                     >  2:1117097/1‑189 (MQ=255)                                                                                                                                                                aaCA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAAt                                                                   >  2:444096/1‑239 (MQ=255)                                                                                                                                                                aaCA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGGGCTCGCATTCAGGAAATGCCTATGCTGACGGTAAAt                                                                   >  2:261978/1‑239 (MQ=255)                                                                                                                                                                 aacgggggggggCGGTGGGTAGTGCGGGAG‑ATGGATTGTAAGCATTTCCCACCGGGGGAACCGGCCGCGATATTGTTAAAGCATTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAAt                                                                   <  2:417396/230‑1 (MQ=255)                                                                                                                                                                 aCA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAAt                                                                   >  2:1213183/1‑238 (MQ=255)                                                                                                                                                                      tggcgggcttggggcatggggggAG‑TTGTATTTTTAGCATTTTCCACCCTTGCAAACGTCTGCGATATTTTTAAAGCAGTTGAAGCCGCCGTTACGCGTTATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGcc                                                               <  2:927015/218‑1 (MQ=255)                                                                                                                                                                         cGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTTTTAGGGACAAAGGTTGGGGAGATCGGCATT‑‑CCTTCGCTCCGTGGGGCTCGCTTTCCGGAAGGGAATATCCGAACGGTAAATTGccgtgg                                                           >  1:902268/1‑235 (MQ=255)                                                                                                                                                                          ggggggtGGGTCAGGGGGGAG‑TTGATTTGTTAGCCTTTTCCCCCCGGGGAAACAGCCGTATTATTTTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCtgtg                                                           <  2:342456/232‑1 (MQ=255)                                                                                                                                                                          ggcgtcTGGGTCATGGGGGAG‑ATGTATTGTGAGCCATTCCCCCCCGTGGCAACAGCCGGGATATTGTTAAGGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCtgtg                                                           <  2:264608/233‑1 (MQ=255)                                                                                                                                                                           ggCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGCGATCCGCATTCCCCTTCGCGTGGTAGTGCGCGCATTCAGGAAAAGCATATGCTGACGGGAAATTGCCtgtg                                                           >  1:116717/1‑239 (MQ=255)                                                                                                                                                                           ggCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCg                                                         >  2:762997/1‑239 (MQ=255)                                                                                                                                                                            cgtgtGGTTCATGGGGGAG‑ATTTATTTTTACCCTTTCCCCCCCGTGGCAACAGCCGGGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCg                                                         <  2:225113/234‑1 (MQ=255)                                                                                                                                                                             cgcTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGGGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAggg                                                                                                                                                >  2:781296/1‑152 (MQ=255)                                                                                                                                                                                  ggTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCAGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACAATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGGTATCCGCAGT‑‑CCTTCGGAGCGGAGCGCTCGCATTACGGAAATGCCTATCCTGACGGTAAATTGCCtggtggatttgat                                                  >  1:264608/1‑228 (MQ=255)                                                                                                                                                                                      aTGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑At                                            >  2:32176/1‑239 (MQ=255)                                                                                                                                                                                       tGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑a                                          >  2:885957/1‑239 (MQ=255)                                                                                                                                                                                            gag‑aTGTCTTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGGACTT‑GCAGGGTTGTTAGGGCCACAGGACGGGGAGACCCGCATT‑‑CCTTTGAATCGTAGTGCTCGCCTTCCGGAAATGCATACTCTGACGGTAAATTAccgggtggctcagttagataaaca                                         >  1:1074603/1‑216 (MQ=255)                                                                                                                                                                                            gag‑aTGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGCTCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCGGGAAAGGCATATGCTGACGGAAAATTGCCCGTGGG‑aagagatcggaaaaa                                         >  1:112871/1‑223 (MQ=255)                                                                                                                                                                                             ag‑aTGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGCGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑Aacac                                      >  2:1209967/1‑237 (MQ=255)                                                                                                                                                                                              g‑aTGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCCCAGGATGTGGGGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGCAAATTGCCTGTGCG‑Atttttccgttacaacgct                                     >  2:235060/1‑222 (MQ=255)                                                                                                                                                                                                    tttGTTAGTCCTCTCCACCCGGGGCACCATCGGGGATTTTGTTAAAGCGTTTAAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTAGGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTcc                              <  1:1218288/238‑1 (MQ=255)                                                                                                                                                                                                        ttAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTTCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTc                           >  2:726657/1‑238 (MQ=255)                                                                                                                                                                                                          aGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGga                                                                                                                                  >  2:352929/1‑138 (MQ=255)                                                                                                                                                                                                          aGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTAAcc                        >  2:207610/1‑239 (MQ=255)                                                                                                                                                                                                          aGCCATTTCCACCCGTGCCAACAGCCGCGCTATTGTTACAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGCCGGCGGCGGACTT‑GCAGGTTTGTTAGGGCCACAGGATTTGGGGATCCGCATT‑‑CCTTCGCTCCGTAGTTCCCGCATTCAGGAAATGCATATGCTGCCGGAAAATTCCCTCTCCG‑ATCTTGTCGAAA‑AACAC‑CTTTTCCCCCAcc                        >  1:884000/1‑239 (MQ=255)                                                                                                                                                                                                           ccccttcccccccGTGGCAACGGGCGATATTTTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCAc                         <  2:116717/229‑1 (MQ=255)                                                                                                                                                                                                            ccATTCCCCGCCGTGGAAACAGCCGGGATATTGTAAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAg                      <  2:751959/239‑1 (MQ=255)                                                                                                                                                                                                             taTTCCCACCCGTGGCAACCGCCGTGAATTTGTTAAAGCAGTTGAAGCCGGCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAgg                     <  2:808818/238‑1 (MQ=255)                                                                                                                                                                                                             cATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGCACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGGTCTGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCCGGAAATGGATATGCTTCCGGAAAATTGCC‑GTGGGAATAGGATCG‑AA‑AACACGCTTTTccccacccggg                     >  1:1021768/1‑231 (MQ=255)                                                                                                                                                                                                             cATTCCCGCCCGTGGCAACAGCCGGGTTTTTGTTAAAGCATTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAg                      <  2:1236884/238‑1 (MQ=255)                                                                                                                                                                                                             cATTCCCCCCCTGGGCAACAGCCGAGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAg                      <  2:1021768/238‑1 (MQ=255)                                                                                                                                                                                                             cATTCCCCCCCGGGGAAACGCCGGCAATATGGTTAAAGCAGTTGAGGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAg                      <  2:556474/238‑1 (MQ=255)                                                                                                                                                                                                               ttcccACCGGGGGCAAAAGCCGGGATATTGTAAAAGGATTGGAGGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGga                    <  2:696471/235‑1 (MQ=255)                                                                                                                                                                                                                ttCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGAATTCAGGAAATGCATATGCGGCCGGTAAATTGCCTGTGCG‑ATCGGATCG‑AT‑AACACGCTTTTCCCTCACCAggcta                  >  2:1128240/1‑236 (MQ=255)                                                                                                                                                                                                                   cccccttgcaaaCCGCCGCATTGTTGTAAAAGCAGTTGAGGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGAtt               <  2:1294352/229‑1 (MQ=255)                                                                                                                                                                                                                     cccGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCCTT‑‑CCTTCGCCTCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAACTTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGAGTGATTaa             >  2:1117195/1‑239 (MQ=255)                                                                                                                                                                                                                      cggTGGAACCAGCCGGGATTTTGTAAAAGCGTTTGAGGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGATTAAg            <  2:626960/237‑1 (MQ=255)                                                                                                                                                                                                                      ccgggacAACAGCCGAGATTTTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGATTAAg            <  2:711660/233‑1 (MQ=255)                                                                                                                                                                                                                         tGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGATTAAGgag         >  2:665929/1‑239 (MQ=255)                                                                                                                                                                                                                         tGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGGTTGTTAGGGCCACAGGATGTGGAGAGCCGCATT‑‑CCTTCGCATCGGAGTGCTCGCATTCAGGAAATGCATATTCTGACGGTAAATTGGCTGGGCG‑AGTTGATCG‑GTAAACAC‑CTTTTCCCTCCCCAGGATGATTAAgggg         >  1:124105/1‑237 (MQ=255)                                                                                                                                                                                                                           gcaaccagcggtatatTGTTAAAGCGGTGGAACCCGCCGTTACGCGTGATATGACCATGGTGGCATTGACCGGCTAGGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGGGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGATTAAGGAGaa       <  1:868867/227‑1 (MQ=255)                                                                                                                                                                                                                                   cgcgTTATTGTTAAAGCTGTTTAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACAt  <  2:207278/236‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                            GTAATGGAACTTCCGCTGCCAATGCACAGCATTTTGCTGCCAGCATGATCAACCGTTTCGAAACGGAGCGGCCCAGCTTACCTGCCATTGCACTAAATACTGATAATGTTGTCTTAACGGCGATTGCCAACGATCGCTTACATGATGAAGTGTATGCAAAACA‑GGTGCGGGCGCTGGGTCATGCGGGAG‑ATGTATTGTTAGCCATTTCCACCCGTGGCAACAGCCGCGATATTGTTAAAGCAGTTGAAGCCGCCGTTACGCGTGATATGACCATTGTGGCATTGACCGGCTATGACGGCGGCGAACTT‑GCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATT‑‑CCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCG‑ATCTGATCG‑AT‑AACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACAT  >  NC_000913/3295949‑3296408

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