Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,529,304	A→C	9.3%	pseudogene (1403/2037 nt)	rhsE →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,529,304	0	A	C	9.3%	257.5 / 3.8	100	pseudogene (1403/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base A (48/40);  new base C (2/10);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAA  >  NC_000913/1529085‑1529522                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           gTGAAACCGGGGCGGGGCTTGCGCGGTGGGGGCTGAGGGGGGAGCAACGGGCCAGACA‑GGTGGAGCCGGAATCCCCCCCGGCGCAAAAAGTTCATTTTTCCCCCT‑GCACCCCCGGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGCACCAGCCg                                                                                                                                                                                                                     <  2:952443/227‑1 (MQ=14)            gCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAACGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGGG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCCTCACCTGGCCCAGCCGTACCGTCTGCCa                                                                                                                                                                                                         >  1:375712/1‑228 (MQ=255)            gCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGCCTGCCa                                                                                                                                                                                                         >  1:478641/1‑228 (MQ=255)                    ggcgaggTGCGGTCTGGCGGGGGAGCAACGGCCAAAACG‑GGGGGAGCCGGAATCCCCCCCGGCGCGAAAAGTTCATTTTTACCACG‑GCGCCCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGCACCAGCCGTACCGTCTGCCAGGgcagca                                                                                                                                                                                                 <  2:999438/222‑1 (MQ=17)                           tGCGGGCTGACGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑GATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatgat                                                                                                                                                                                          >  1:824706/1‑228 (MQ=255)                             cGGGCTGACGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGGGAACCCGCATCACCTGCCCCAGCCGAACCGTCTGCCCAGGCAGCAGCaggattg                                                                                                                                                                                        >  2:152188/1‑226 (MQ=255)                                    aCGGTGGAGAACGTGGCCAGCAA‑GGTGGGCCCGAAATCCACACGGGCGAAAAAAGTTATTTTTTCCCCCG‑GCGCCCCCGGGGCCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCg                                                                                                                                                                                 <  2:913090/228‑1 (MQ=17)                                     cGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgca                                                                                                                                                                                                    >  1:1902807/1‑208 (MQ=255)                                     cGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGACGTAACGTCTGCCAGGGCAGCAGCATGATGTGGCGCCgg                                                                                                                                                                                >  2:1039531/1‑228 (MQ=255)                                     cGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCCTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGCCTGCCCGAGCAGCAGCctgat                                                                                                                                                                                          >  2:1073666/1‑218 (MQ=255)                                     cGGTGGACCAACTGCCAAGCCA‑GGTGGCGCCGAAATCCCCCCCGGGGCGAAAATTTCATTTTTACCCCG‑GCGCCCACGGGGGCCGGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCg                                                                                                                                                                                 <  2:1145931/227‑1 (MQ=17)                                      ggTGGAGAAACTGCCCAGACA‑GGTGGACCCGAAATACCCCCCGGCGAGAAAAGTTCATTTTTCCCCCG‑GCGCCCCCGGGGCCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAAGGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcag                                                                                                                                                                                                   <  2:1899440/208‑1 (MQ=17)                                       gTGGACAAACTGGCCAGACG‑GGTGGCGCGGAAATCCACCCCGGCGCGAAAAGTTCATTTTTCCCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAGTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgg                                                                                                                                                                                <  1:1161053/226‑1 (MQ=17)                                         ggcacAACT‑GCCCGACAGGGTGGACCCGAAAGCCACACCGGCGAAAAAAGTTCATTTTGCCCACG‑GCGCCCCCCGGGGCCGGCCGGTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGAGGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                                                             <  2:375094/223‑1 (MQ=17)                                         ggAGAAACTGCCCAGACA‑GGGGGAGCGGAAATACCCCCCGGCGCGAAAAGTACATTTTTACCACG‑GCGACCACGGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAAGGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                                                             <  1:1075589/227‑1 (MQ=17)                                         gagacatgtGGCGAGACG‑GGGGGGCCGGGAATACCCACGGGCGGGAAATGTTCATTTTTCCCCCG‑GCGCCCACGGGGGCCTGCCGCTGGCGCTCAGCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAAGGGGGCAACCAGCTTAAGGAGGAGACCCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                                                             <  2:1776807/219‑1 (MQ=17)                                         cgagcactgggCCAGACA‑GGTGGGCCCGAAATCCACCCCGGCGCGAAAAGTGCATTTTTACCACT‑GCGCCCCCCGGGGCCTGCCGCTGGCGCTCAACAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                                                             <  2:398785/218‑1 (MQ=17)                                              aaCTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGGACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTAc                                                                                                                                                                        >  1:765860/1‑227 (MQ=255)                                                  ggCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGGACCGTCTGCCCGGGCAGCAGCATGATGAGGAGTCGGGGCTGtgc                                                                                                                                                                        >  2:306109/1‑221 (MQ=255)                                                       gACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTc                                                                                                                                                              >  1:1642293/1‑228 (MQ=255)                                                              tGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCATCATGATGAGGAGTCGGGGGCGGACCATAACCGTCAc                                                                                                                                                            >  2:1410173/1‑224 (MQ=255)                                                              tGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGGACTATAACCGTCaa                                                                                                                                                            >  1:1016951/1‑223 (MQ=255)                                                                   ccGGACTACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGCGGCTGTACTCTAACCGGCACCGGCACTAc                                                                                                                                                   >  1:2023007/1‑228 (MQ=255)                                                                         accccgccGGCGGGAAAAGTTCATTTTGACCACG‑GCGACCCCGGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTa                                                                                                                                                    <  1:1062986/215‑1 (MQ=255)                                                                                gggcggcaaaaTTTCTTTTTTACCCCTGGCGCCCCCGGGGCCCTGCCGCTGGCGCTCATCGGCGAAGCCGGC‑AATGC‑GGCGGGGCGCGG‑GGAG‑TATGATGAAGGGGGCAACCAGCTTAATGGGGAGAACCCGCATCCCCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAg                                                                                                                                       <  2:1046201/221‑1 (MQ=17)                                                                                cGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgca                                                                                                                                                                                                    >  1:523215/1‑166 (MQ=255)                                                                                 ggCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGTGGCTGCACTATAACCGTTACCGGCACTACGATCCGTTGCAgg                                                                                                                                      >  2:696267/1‑227 (MQ=255)                                                                                 ggCGCAAAAAGTTCATTTTTCCCACT‑GCCACCCCGGGGCCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAAGGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcag                                                                                                                                                                                                   <  2:522354/166‑1 (MQ=17)                                                                                   cgcgAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAATACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGCATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTc                                                                                                                                                                                                              >  1:942380/1‑153 (MQ=255)                                                                                   cgcgAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCggg                                                                                                                                                                               >  1:1174655/1‑184 (MQ=255)                                                                                   cgagAAAAGTTCTTTTATCCCACG‑GCGACCACGGGGCCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATCC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCACCCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                                                                                    <  2:764567/224‑1 (MQ=17)                                                                                    gcgAAAAGTTCATTTTTCCCCCT‑GCGCCCCCCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGCACCAGCCGTACCGTCt                                                                                                                                                                                                             <  2:940709/153‑1 (MQ=255)                                                                                    caaaacgtttCATTTTCCCCGCG‑CCGACCACGGGGACCTGCCGCTGGCACTCATCAGCGAATACGGC‑AATAC‑GGGGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCg                                                                                                                                  <  2:1767454/220‑1 (MQ=17)                                                                                     cgAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTc                                                                                                                                                                                  >  2:1788599/1‑179 (MQ=255)                                                                                      gAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCg                                                                                                                                                                                 <  1:1791768/179‑1 (MQ=255)                                                                                            ttCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCa                                                                                                                                        >  1:1048088/1‑214 (MQ=255)                                                                                             tCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCa                                                                                                                                                             >  2:364069/1‑192 (MQ=255)                                                                                              cATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCAc                                                                                                                                                            <  1:364676/192‑1 (MQ=255)                                                                                                      ccACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGGG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                                                              >  2:1193045/1‑166 (MQ=255)                                                                                                       cACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGGG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                                                             <  1:1195186/166‑1 (MQ=255)                                                                                                          t‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCa                                                                                                                                                             >  1:81117/1‑179 (MQ=255)                                                                                                            gCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGCATGGGGCAACCCGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGCTGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGAt                                                                                                                                                >  1:1728304/1‑191 (MQ=255)                                                                                                            gCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCg                                                                                                                                                                >  1:1767293/1‑175 (MQ=255)                                                                                                             cgcccccgggggCCTGCCGCGGGCGCTCAGCAGCGAAGACGGC‑AATAC‑GGCGGGGCGCGG‑GGAG‑TATGAGGAAGGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCGGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATc                                                                                                                                               <  2:1725246/183‑1 (MQ=17)                                                                                                                     ggggCCTGCCGCTGGCGCTCATCAGCGAAGACGGA‑AATCC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGAAACCAGCTTAATGGGGAGACCCCGCATCCCCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGagag                                                                                                  <  2:1121531/228‑1 (MQ=255)                                                                                                                      gggCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCa                                                                                                                                                             >  1:1927621/1‑168 (MQ=255)                                                                                                                       ggCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcag                                                                                                                                                                                                >  1:1494246/1‑132 (MQ=255)                                                                                                                        gCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagcagc                                                                                                                                                                                               <  2:1491577/132‑1 (MQ=255)                                                                                                                           ggCCGCTGGCGCTCATCAGCAAAGACGGA‑AATAC‑GGCCTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAAGGAGGAGAACCCGCATCACCGGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtgga                                                                                             <  2:2019395/226‑1 (MQ=255)                                                                                                                              ggCTGGCGTTAACCAGCGAAG‑CGGC‑AAAGCGGGCGGGGCGCGG‑GGAT‑TATGAGGAGTGGGGAAACCAGCTTAAGGGGGAAACCCCGCATCCCCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggatgga                                                                                         <  2:158513/227‑1 (MQ=17)                                                                                                                              aggtGGCGCTCATCCGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GAAG‑TATGATGAATGGGGCACCCGGCTTAAGGAGGAGAACCCGCATCACCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggatgga                                                                                         <  1:1468433/225‑1 (MQ=255)                                                                                                                               ggtGGCGCCCACGAGCGAAGGCGGCAAATGC‑GG‑GGGGCGCGGCGGAT‑AATGATGAAGGGGCAAACCGGTTTAATGAGGAGACCCCGCATCCCCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggatgga                                                                                         <  2:1475944/226‑1 (MQ=17)                                                                                                                                 tGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCCGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAACCCGGACCCGATTGGGTTGACAGGtgga                                                                                             >  1:963635/1‑221 (MQ=255)                                                                                                                                          tcaGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGt                                                                                                                                                                                   <  2:1168521/126‑1 (MQ=255)                                                                                                                                          acaGCGAAGACGCA‑CATCC‑GGCGTGGCGCGG‑GAGT‑TAGGAAGAAGGGGCCAACCAGCTTAAGGAGGAGACCCCGCATCCCCGGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGt                                                                             <  2:1015095/227‑1 (MQ=17)                                                                                                                                             gcaaaGGCGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAAGGAGGAGACCCCGCATCCCCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTAt                                                                           <  2:823284/224‑1 (MQ=255)                                                                                                                                             aggAAGACGAA‑CATCC‑GGCGTGGCGCGG‑GGAT‑GATGATGAATGGGGCAACCGGCTTAAGGAGGAGAACCCGCATCACCTGCCCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTAt                                                                           <  2:477866/225‑1 (MQ=17)                                                                                                                                                 gaccGGC‑AATGC‑GGCGGGGCGGGG‑GGAT‑TATGATGAATGGGGCACCCGGTTGAAGGAGGAGAACCTGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGtt                                                                      <  2:1641000/225‑1 (MQ=18)                                                                                                                                                   aCGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                                                                                                    >  2:1858977/1‑116 (MQ=255)                                                                                                                                                   aCGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATAACCCCGGAACCGATTGGGTTGTGAGGTGGAGGGAATATGGACTAGGATCCGTTgg                                                                    >  2:537121/1‑227 (MQ=17)                                                                                                                                                    cGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGt                                                                                                                                                                                   <  1:1862267/116‑1 (MQ=255)                                                                                                                                                                 cGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCCCCTGACCCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACaa                                                         <  2:1195496/227‑1 (MQ=255)                                                                                                                                                                 cGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAg                                                        >  2:771878/1‑228 (MQ=255)                                                                                                                                                                     ggcacg‑ggAGTTATGATGAATGGGGCCACCAGCTTAATGGGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCAGACaa                                                         <  1:1041413/219‑1 (MQ=255)                                                                                                                                                                     gctcgG‑GGAG‑TAGGATGAATGGGCCAACCAGTTTAATGGGGAGAACCCGGACCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACa                                                          <  1:1412707/219‑1 (MQ=255)                                                                                                                                                                                 tATGATGAATGGGGCAACCAGTTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAAt                                           <  1:1056767/227‑1 (MQ=255)                                                                                                                                                                                 tATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGc                                                                                                                                   >  2:1132871/1‑139 (MQ=255)                                                                                                                                                                                  atgatgAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCg                                                                                                                                  <  1:1134905/139‑1 (MQ=255)                                                                                                                                                                                  atgatgAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGGGAGGTGGATGGAATATGTATCAGTATCCGTTGAATTCCCTACAAGTGATAGCCCCAATgg                                         >  1:2005186/1‑228 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                                                              >  2:217756/1‑88 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                                                              >  1:499537/1‑88 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCg                                                                                                                                                                >  2:700620/1‑102 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAggg                                                                                                                                     >  1:2040242/1‑129 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCg                                                                                                               >  2:71239/1‑151 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGa                                                                    >  2:200365/1‑194 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGa                                  >  2:274622/1‑228 (MQ=255)                                                                                                                                                                                         aTGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGa                                  >  2:757615/1‑228 (MQ=255)                                                                                                                                                                                          tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                                                             <  1:218128/88‑1 (MQ=255)                                                                                                                                                                                          tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                                                             <  2:498727/88‑1 (MQ=255)                                                                                                                                                                                          tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCt                                                                                                                                                                            >  1:1521185/1‑89 (MQ=255)                                                                                                                                                                                          tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGt                                                                                                                                                               <  1:701795/102‑1 (MQ=255)                                                                                                                                                                                          tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                                                                                    <  2:2036597/129‑1 (MQ=255)                                                                                                                                                                                          tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGa                                                                                                              <  1:71362/151‑1 (MQ=255)                                                                                                                                                                                          tGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGaa                                                                   <  1:200707/194‑1 (MQ=255)                                                                                                                                                                                           ggggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCt                                                                                                                                                                            >  2:639381/1‑88 (MQ=255)                                                                                                                                                                                           ggggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTg                                                                                                                                                                           <  2:1518468/89‑1 (MQ=255)                                                                                                                                                                                           ggggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAggg                                                                                                                                     >  1:808581/1‑127 (MQ=255)                                                                                                                                                                                            gggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTg                                                                                                                                                                           <  1:640440/88‑1 (MQ=255)                                                                                                                                                                                            gggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                                                                                    <  2:807197/127‑1 (MQ=255)                                                                                                                                                                                             ggCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCg                              >  1:234826/1‑228 (MQ=255)                                                                                                                                                                                               cAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCTTACAAGTGATAGACCCAATGGGGTTAGATGCGa                             >  1:237195/1‑227 (MQ=255)                                                                                                                                                                                                      cTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGGTGCGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACACGTGATAGACCCAATGGGGTTAGATGCGCTTGAGAat                     >  1:2026856/1‑228 (MQ=255)                                                                                                                                                                                                        tAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAatat                   >  1:673432/1‑228 (MQ=255)                                                                                                                                                                                                            ggggagACCCCTCATCACCGGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACa               <  2:961924/226‑1 (MQ=255)                                                                                                                                                                                                                agaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcag                                                                                                                                                                                                   >  2:1936188/1‑44 (MQ=255)                                                                                                                                                                                                                agaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGc                                                                                                                                                        >  2:1492785/1‑87 (MQ=255)                                                                                                                                                                                                                 gaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagc                                                                                                                                                                                                  <  1:1939627/44‑1 (MQ=255)                                                                                                                                                                                                                 gaACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCa                                                                                                                                                       <  1:1495456/87‑1 (MQ=255)                                                                                                                                                                                                                         aTCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGACACTACGATCCGTTGCAGGGGCGGGATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGGGAATATGACATCAGGTGGACTaa  >  1:961408/1‑228 (MQ=255)                                                                                                                                                                                                                              cTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAg           <  1:758894/214‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                           GTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACA‑GGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACT‑GCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGC‑AATAC‑GGCGTGGCGCGG‑GGAG‑TATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATATGACATCAGGTGGACTAA  >  NC_000913/1529085‑1529522

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