Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,582,478 A→C 13.3% S191A (TCG→GCG)  hsdM ← DNA methyltransferase M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,582,4780AC13.3% 179.6 / 2.8 78S191A (TCG→GCG) hsdMDNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (52/11);  new base C (2/9);  total (55/23)
Fisher's exact test for biased strand distribution p-value = 5.99e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGA  >  NC_000913/4582250‑4582715
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
ggTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑AACCAGAGTGTTGCCCAGACGGATTGCGCCGCCCTGGGCGGGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGGGCCAGGAGGCGGGTGCCGGGCACCCGTTCGAGGCCGATAAA‑CGCGCGGCGGGGCTGGAAATCC‑TGCGTGTCGCCATCAAGGGCGTCCAGGTCAGTGGTTTG‑GGCCCTAACCTAg                                                                                                                                                                                                                                     >  2:1187335/1‑239 (MQ=255)
           aCAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTtc                                                                                                                                                                                                                          >  1:324157/1‑239 (MQ=255)
                  gcgcCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCaaa                                                                                                                                                                                                                   >  2:233464/1‑239 (MQ=255)
                    gcCTTCGGGA‑GTTTTCCCCGCCGGTGACCCGGGGGTTTGCCCGGGCGGAGTCCACCCCGGGGGCCGAGTTTGCCTTAAATACCGTGCGGCAGGCAGTTCATCAGTGCCAGACGACGGGTCCCGGGCGCCAGTCCGAGGCGGATAAA‑CGCGCGGTGGTTCTGGAAATCC‑TGCGTGTCCCCATCAAGGTCGCCCAGATCATTGGTTTG‑CGCCTTAACATGGCGGTCGGTTTCAATCaaaaa                                                                                                                                                                                                                 <  2:967234/239‑1 (MQ=255)
                      cTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCCATCAAAAAg                                                                                                                                                                                                                >  1:438782/1‑238 (MQ=255)
                       ttCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCTAGGCAGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGAGTGTAGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGc                                                                                                                                                                                                               >  1:435751/1‑238 (MQ=255)
                       ttCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGc                                                                                                                                                                                                               >  2:209439/1‑238 (MQ=255)
                         cGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCG‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGATTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGccc                                                                                                                                                                                                             >  2:940228/1‑238 (MQ=255)
                         cGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGGCCAGATCATTGGTTTG‑CGACTTCACATAGCGGTCGGCTTCACCCAAAAAGccc                                                                                                                                                                                                             >  2:4963/1‑238 (MQ=255)
                         cGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGCTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACCACGGGCGCCGGGCACCAGGTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCACCAAGGTCGTCACGGTCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCGATCAAAAAGccc                                                                                                                                                                                                             >  2:957837/1‑238 (MQ=255)
                         cGGCAGGTTTTCACCGTCCCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGACGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGCCGGCTTCAATCACAAAGccc                                                                                                                                                                                                             >  2:297299/1‑238 (MQ=255)
                                     gccgccgCT‑CCCCAGAGTGTGGCCCGGACGGTTTGCGCGGCGGTGGGCGAGGTTGCCTTCAAGAGCGGGCAGCAGGCGTTTAAGCAGTGCCAGAGGACGGTGGCGGGGCCCCGGTTCGAGGCCGATAAC‑CGCGGGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGCTTTACCAGGGCGGTCGGCTCCAATAAAAAAGCCCGCCGTACCTgccg                                                                                                                                                                                                <  2:1100383/238‑1 (MQ=255)
                                         tCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGCCTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg                                                                                                                                                                                            >  1:222693/1‑239 (MQ=255)
                                         tCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg                                                                                                                                                                                            >  1:207798/1‑239 (MQ=255)
                                               aCCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATGAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCAGTACCTGCCGCCGGGt                                                                                                                                                                                          >  1:562834/1‑236 (MQ=255)
                                                cccAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTACCATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGCCCt                                                                                                                                                                                       >  2:1009603/1‑238 (MQ=255)
                                                  cAGAGGGTAGCCCAGAGGGATGTCCCCACCGTGGTGGAGGTCGCCTCCAATCTGGTGCAGCAGGCAGTTAATCAGGGCCGGAGGGCGGGTGCCGGGCCCCAGTGGGAGGAATAAAAC‑CGGGGGGTGGTTTTGAAACCCT‑TGCGGTTCGCCTTAAAGGTCGTCCAGTTCTTTGGTTGG‑CGCTTTAACAGCGCGGTGGGCTTCAATAAAAAACCCCGCCGTACCTGCCGCCGGGTCCTGca                                                                                                                                                                                    <  2:223953/239‑1 (MQ=255)
                                                     agTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGcacc                                                                                                                                                                                  >  1:925305/1‑238 (MQ=255)
                                                     agTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGca                                                                                                                                                                                    >  1:611577/1‑236 (MQ=255)
                                                          tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGGCTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTc                                                                                                                                                                             >  1:132191/1‑238 (MQ=255)
                                                          tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTc                                                                                                                                                                             >  1:150387/1‑238 (MQ=255)
                                                          tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGGGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGGGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGATGGATCTGGAAATCC‑TGCGGGTCGCCATCAAGGTCGTCCCGTTCATTGGTTTG‑CGACTTAACATATCGGTCGGCTTCAGTCACAAAGGCCGCCGAACCGGCCGCCGGGCCCTGCACCACTTc                                                                                                                                                                             >  2:1033658/1‑238 (MQ=255)
                                                          tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGCTGCCTTCAATATCGTGCAGCAGGCAGATCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCACGGTCGTCCAGATAATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAACCCCGCCGACCCTGCCGACAGGTCCTGCACCACTTc                                                                                                                                                                             >  2:932492/1‑238 (MQ=255)
                                                                  gggagcgcgcCGCCGTGTGCAGGGTGGCCTGCAATAGCGTGCGGCAGGCGGTTCATCAGTGCGAGAGGAGGGTTGCCGGGCCCCAGTTAGAGGCCGACAAA‑CGCGGGGTGGTTCTGGAAACCT‑TGCGTTTCCCAAGCAAGGTCGCCCAAATCTTTGGTTTG‑CGCCTAAACATAGCGGCGGGCTAAAACCAAAAGCCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCt                                                                                                                                                                     <  2:1128621/232‑1 (MQ=255)
                                                                           ccgccgTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGCTTcag                                                                                                                                                           >  2:672441/1‑239 (MQ=255)
                                                                               tgTGGTCGGGGTGGCCTGCAAGAGCGTGCCGCAGGCAGTACATCAGGGCGAGACGAGGGGTGCCGGGCGCCGGTGCGGGGCGGATAAA‑CGCGGGGTGGATCTGAAAATCC‑GGCGTGTCGCAACCAAGGTCATCCAGATCGTTGGTTTG‑CGCTTTAACATGGCGGTCGGCTTCAAACAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTcagcag                                                                                                                                                        <  2:595876/237‑1 (MQ=255)
                                                                               cggggTCGGGGGTGCCTTCAATATGGCGCAGCAGGCGGTAAACAAGTGCCAGCGGAGGGGTGCCGGGCACCAGTTCGAGGCGGAAAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAAATCATTGGTTTG‑CGATTTAACATAGCGGGCGGCTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTcagcag                                                                                                                                                        <  2:925305/235‑1 (MQ=255)
                                                                               cgTGGTCGGGTTTGCCTACAATATCGTGCAGCAGGCAGTTCATCAGGGCAAGCGGCCGGGTGCCGGGCACCAGTTCGAGGCAGATAAA‑CGCGGGGTGGATCGGGAAATCC‑TGCGTGTCGCCATAAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTcagcag                                                                                                                                                        <  1:1046574/238‑1 (MQ=255)
                                                                                           tGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGtt                                                                                                                                           >  2:345422/1‑239 (MQ=255)
                                                                                             ccTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGACTAATGGtttt                                                                                                                                         >  2:768657/1‑239 (MQ=255)
                                                                                                  aaTATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCa                                                                                                                                    >  1:360817/1‑239 (MQ=255)
                                                                                                  aaTATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCa                                                                                                                                    >  1:446189/1‑239 (MQ=255)
                                                                                                         tGCAGCAGGCAGTTCCTCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAACCCGCCGTACCTGCCGCCGGGTCCTGCCCCACTTCACGCGGCGGCGGCTTCAGCAGATGAGTAATGGTTTTAATCAGCggca                                                                                                                              >  1:614125/1‑236 (MQ=255)
                                                                                                         tGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGTGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCTCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTAAGCAGATGAATAATGGGTTTAATCAGCGGAc                                                                                                                              >  2:145302/1‑238 (MQ=255)
                                                                                                         tGCAGCAGGCAGTTCATCAGTGCCAGAAGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGTTCGGCTTCAATCAAAACGCCCCCCGTCCCTGCCGCCGGGTCCTGCACAACTACACGCGGCTGCGGTTTCAGCAGATGACTCCTGGCTTTAATCAGCGGAc                                                                                                                              >  2:1014387/1‑238 (MQ=255)
                                                                                                         cacagcagGCCGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TCCGTGTCGCCATCAGGGGCGTCCAGATCATTGGTTTG‑CGACTTAACTTAGCGCTCGTCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCCGCGGATTCAGAAGTTGAACAATGGTTTTAATGATCGGAc                                                                                                                              >  2:716042/3‑238 (MQ=255)
                                                                                                          gcaggagagagTGCATCGGGGCCAGAGGACGGGGGACGGGCGCGAGGCCGA‑GCCGAAAAAGCGCGCGTTGGTTCTGAAACTCC‑TGGTTGTCGCCATCAGGGTCGTCAAAATCATGGGTTGG‑CCACTAAACATGGCGGTCGGCTACAACAAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGAcg                                                                                                                             <  2:484379/229‑1 (MQ=255)
                                                                                                                 gcgtttCAGCAGGGCCAGACGGGGGGGGCGGGACCCCGGTTGGAGGCCGATAAC‑CGGGGGGTGGATCGGGAAACCC‑TGCGTGCCCCCATAAAGGCCGTCCATATCATTGGTTTG‑CGACTTAACATAGCGGTCGGTTTAAACAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTg                                                                                                                      <  1:1057285/234‑1 (MQ=255)
                                                                                                                 gcagTTCATCAGTGCCAGACGACGTGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATGC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGGTTG‑CGACTTAACATAGCGGTCGGCTTCAGTCAAAAAGCCGGCCGTACCTGCCGCCGGGCCATGCCCCACTTCACGCCGCCGCGGTTTCGGCAGATGAATAATGGTTTTAATGAGCGGACTCGGGGTgt                                                                                                                     >  2:870489/1‑238 (MQ=255)
                                                                                                                     ttCGTCAGTGCCAGGCGGCGGGGGCCGGGCCCCCGCGCGAGGCCGATAAC‑CGGGGGTTGGAT‑TGGAAACCCGCGCGTTTCGCCACCAGGGTCTACCAGATCGTTGGTTGG‑GGTCCAAACATAGCGGTCGTCTTCAACCAAAAACCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTa                                                                                                                 <  2:72715/239‑1 (MQ=255)
                                                                                                                               ccAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGc                                                                                                        >  1:228893/1‑238 (MQ=255)
                                                                                                                               ccAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGAGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGATTCAGCAGATGCATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGc                                                                                                        >  1:789800/1‑238 (MQ=255)
                                                                                                                                     gaggggTGCCGGGCCCCAGTTCGAGGCGTATACA‑GCGGCGGTGGGTCT‑GAAATCC‑TGGGTGGCGCCATCAAGGCCGCCAAGATCATGGGTTTG‑CGTCTTAACATGGCGGTCGGCTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGAc                                                                                                 <  2:1156683/235‑1 (MQ=255)
                                                                                                                                       cgGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACtt                                                                                               >  1:523881/1‑239 (MQ=255)
                                                                                                                                       cgGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGAACTGGCCTGCACCCGACtt                                                                                               >  2:255645/1‑239 (MQ=255)
                                                                                                                                       cgGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGACGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAGAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGGGAAGTACTGGCCTGCACCAGACtt                                                                                               >  2:622679/1‑239 (MQ=255)
                                                                                                                                             ccGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTc                                                                                         >  1:1129178/1‑239 (MQ=255)
                                                                                                                                                    ccAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCTGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTcgcg                                                                                  <  1:304458/239‑1 (MQ=255)
                                                                                                                                                    ccAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTcgcg                                                                                  >  1:76868/1‑239 (MQ=255)
                                                                                                                                                                 gATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACa                                                                     <  1:209439/239‑1 (MQ=255)
                                                                                                                                                                  aTAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACa                                                                     >  1:1037038/1‑238 (MQ=255)
                                                                                                                                                                  aTAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACa                                                                     >  2:708899/1‑238 (MQ=255)
                                                                                                                                                                               tGGATTTGAAAATCC‑GGCGTGTCCCCAACACGGCCGTACAGATTTTTGGTTCG‑CGACTTAACATGGGGGCCGTCTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACa                                                         <  1:1187335/238‑1 (MQ=255)
                                                                                                                                                                                         aaTCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATACTGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAg                                              >  2:63762/1‑239 (MQ=255)
                                                                                                                                                                                            cc‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGCTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATCTCGCCGAAGTCa                                           >  2:874486/1‑239 (MQ=255)
                                                                                                                                                                                             c‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCa                                           >  1:718341/1‑238 (MQ=255)
                                                                                                                                                                                                  gtgtCGCCATCAAGGTTCTCCAGAACTTTGGTTTG‑CGTCTTAACATAGCGGTCGGCTACAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATcgc                                       <  1:552812/238‑1 (MQ=255)
                                                                                                                                                                                                  gtgtCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACatat                                                      >  1:507720/1‑223 (MQ=255)
                                                                                                                                                                                                       gCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACTCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGCCt                                 >  1:642095/1‑239 (MQ=255)
                                                                                                                                                                                                               ggTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgc                         >  1:461502/1‑239 (MQ=255)
                                                                                                                                                                                                                agagcaaggATCATGGGTTTG‑CGACTAAACATGGGGGTCGCTTTAAACAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  2:835038/231‑1 (MQ=255)
                                                                                                                                                                                                                 tggtcCAGATCTTTGGTTTG‑GGCTTAAACCTAGGGGTGGGCTACAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  1:1163259/236‑1 (MQ=255)
                                                                                                                                                                                                                 tcgccCATATCATTGGTTGG‑CACTTTAACATAGCGGTCGGCTACAACCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  1:811559/238‑1 (MQ=255)
                                                                                                                                                                                                                 tcgccCAGTTCA‑TGGTTTGCCGCCTTAACATGGCGTCCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  1:833559/238‑1 (MQ=255)
                                                                                                                                                                                                                  cctcCATATCTTTGGTTCG‑CGATTTAACATAGGGGTCGTCTTAAATCAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcgc                       <  1:230073/236‑1 (MQ=255)
                                                                                                                                                                                                                     ccAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATg                                                                                                                                              >  1:1039108/1‑116 (MQ=255)
                                                                                                                                                                                                                          tCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTCCATATCGCCGAAGTCCTCGCCCGACTTACCGTGCGCGCCGTTga               >  2:310327/1‑237 (MQ=255)
                                                                                                                                                                                                                          tCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGt               >  1:493208/1‑238 (MQ=255)
                                                                                                                                                                                                                           ctttGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAAACAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTAc             <  1:4963/237‑1 (MQ=255)
                                                                                                                                                                                                                               gttttG‑CGCCTAAACGGAGCGGGGGGCTTAAACAAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGt         <  1:957837/237‑1 (MQ=255)
                                                                                                                                                                                                                               ggTTTG‑CGACTTAACATGGCGGTCTGCTTCAAACAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCGTTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGt         <  1:940228/239‑1 (MQ=255)
                                                                                                                                                                                                                                   tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTa              >  1:835038/1‑230 (MQ=255)
                                                                                                                                                                                                                                   tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  1:941535/1‑238 (MQ=255)
                                                                                                                                                                                                                                   tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  1:774718/1‑238 (MQ=255)
                                                                                                                                                                                                                                   tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  1:762383/1‑238 (MQ=255)
                                                                                                                                                                                                                                   tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  2:949551/1‑238 (MQ=255)
                                                                                                                                                                                                                                   tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGCGCGCGCCGTTGTACCAGTCCa      >  2:14154/1‑238 (MQ=255)
                                                                                                                                                                                                                                      cGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCg   >  1:464472/1‑239 (MQ=255)
                                                                                                                                                                                                                                       caCTTACCAGTGCGGTCGGCTCAAAACAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGa  <  1:815903/238‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
GGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGA  >  NC_000913/4582250‑4582715

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