Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 326,812 T→C 15.5% S146G (AGC→GGC)  betA ← choline dehydrogenase, a flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913326,8120TC15.5% 185.6 / 2.1 85S146G (AGC→GGC) betAcholine dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base T (51/18);  new base C (4/12);  total (55/30)
Fisher's exact test for biased strand distribution p-value = 4.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT  >  NC_000913/326575‑327048
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
tgtgATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGCTTGACGCCGGGCTTGGAGGCAGT‑GACGCt                                                                                                                                                                                                                                                >  1:267479/1‑238 (MQ=255)
tgtgATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGGGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGCGGTAGT‑GACGCt                                                                                                                                                                                                                                                >  1:410816/1‑238 (MQ=255)
  tgATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTCCGCGGGTAGCCCGCCTGCACGACCGCTTCAATCATCGCTTCAACCAGCGGCTTGACGCCCGCTTTGGAGTGTGT‑CACGCT‑Ca                                                                                                                                                                                                                                             >  2:120323/1‑238 (MQ=255)
    aTCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCCTCGGGCCAAAACCTTCCTGCTGATAACC‑GCTGAGATCGTCCGTGCGCGGGGAGCCCGCCTGCACGCCCGCCTCAATCATCGGTTCAAACAGCGGATTGACGCCGGGTTTTGAGGTAGA‑GAAGCC‑CACCg                                                                                                                                                                                                                                          >  2:753094/1‑239 (MQ=255)
    aTCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGCGGTAGT‑GACGCT‑CACCg                                                                                                                                                                                                                                          >  1:1418045/1‑239 (MQ=255)
         tCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCCCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCa                                                                                                                                                                                                                                     >  2:877363/1‑239 (MQ=255)
          cATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGCAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGTGGTAGT‑GACGCT‑CACCGGGCCa                                                                                                                                                                                                                                     >  2:1260928/1‑238 (MQ=255)
           aTAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGCACCAAAACCTTCCTGCTGACAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCATCAATCATCGCTTCAAACAGCGGGTTTACGCCGCGTTCGGAGGTAGT‑GCCGAC‑CACCGCGCCA‑t                                                                                                                                                                                                                                   >  2:1001109/1‑238 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGt                                                                                                                                                                                                                          >  1:77986/1‑239 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGt                                                                                                                                                                                                                          >  1:1115695/1‑239 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGt                                                                                                                                                                                                                          >  1:442371/1‑239 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGt                                                                                                                                                                                                                          >  1:1068610/1‑239 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGt                                                                                                                                                                                                                          >  1:625382/1‑239 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGt                                                                                                                                                                                                                          >  1:1211698/1‑239 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCC‑CACCGGGCCA‑TCACCGCCGt                                                                                                                                                                                                                          >  2:1388970/1‑239 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGAAGT‑GACGCT‑CACAGGGCCA‑TCAccgccgg                                                                                                                                                                                                                          >  1:268570/1‑238 (MQ=255)
                   aGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGGGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCCCCGCCGt                                                                                                                                                                                                                          >  2:584004/1‑239 (MQ=255)
                          aTGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTc                                                                                                                                                                                                                   >  2:1157737/1‑239 (MQ=255)
                                gggTTGGGACGCGTTTTGGCCGGAAGGGGATAGCC‑CGCCGGGTGGTTGGCGCGCGGGCCCTGGGGCGTGACGTGGCG‑ATCAAGCGACCCAAACCCTTCCTGTGGAAACCC‑TTGGAGACGGCCGGGGCGGGGTGACCCCCCCGCCCCGCCCTCTTAAATCACCGCTTCAACCAGCGGATTGAGGCGGGTTTGGGAGGTGGT‑GCCGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGtttt                                                                                                                                                                                                              <  2:1406107/237‑1 (MQ=255)
                                  gTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCa                                                                                                                                                                                                            >  1:411953/1‑238 (MQ=255)
                                       gACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCAccgcc                                                                                                                                                                                                                            >  1:1039460/1‑217 (MQ=255)
                                       gACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCata                                                                                                                                                                                                      >  1:1015160/1‑239 (MQ=255)
                                                 ggCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtc                                                                                                                                                                                             >  1:246852/1‑238 (MQ=255)
                                                 ggCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑AACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGtc                                                                                                                                                                                             >  2:576289/1‑238 (MQ=255)
                                                        ctgagtggcccccgcgcgGT‑GCGGGCGCGCCGCCCCGGGGGGGTGGGGTTGCG‑ACCCACCGGACAAAACCCTCCCTTCGAATACCC‑TTGGAACGCGCCCGGGCGGGGGAGCCCCGCTGGCACGCCCGCTTCAATCACCGCTTAAAACAGCGGATGGACGCGGGGTTGGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCt                                                                                                                                                                                      <  2:1444216/226‑1 (MQ=255)
                                                            gaTAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCgg                                                                                                                                                                                 >  1:250248/1‑239 (MQ=255)
                                                                 ccACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGgtagta                                                                                                                                                                            >  1:1444216/1‑239 (MQ=255)
                                                                      cgcgGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggca                                                                                                                                                                       >  1:94451/1‑239 (MQ=255)
                                                                      cgcgGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggca                                                                                                                                                                       >  1:1368909/1‑239 (MQ=255)
                                                                      cgcgGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCC‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGCGCCTCGGCCTTGCGGCGGTAggccg                                                                                                                                                                       >  2:1372466/1‑236 (MQ=255)
                                                                           t‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGgcaggcag                                                                                                                                                                  >  1:110093/1‑239 (MQ=255)
                                                                           ggtggggCGCGGGGCCCCTGGGGGGTGGCGGGGCG‑CTCCAGGGGCCCAAACCCTCCTGGCTGATAACG‑GTGGAGATCGCCCGTGCGGGGGTGGCCCGCCGGCCGCCCCGTTTAAATCACCGCTTAAAACAGGGGATTGACGCGGGGTTGGGGGGTAGT‑GACGCT‑CCCGGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGggcaggc                                                                                                                                                                    <  2:583140/233‑1 (MQ=255)
                                                                                   cgcgcgcccgggcggcgTGGCGG‑GCGCATCCAGCGGCCCAACCCCCTCGCGCTG‑TAACCTGTTGAGTGCCCCCGAGCGCGGCTGGCCCGCCTGCACGCCCCCTTCAATCACCGCTCCAAGCCGCGGATTGACGCGGGGTTGGGAGGTAGT‑GACGCT‑CACCGGCCCC‑TCGCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGt                                                                                                                                                           <  2:548588/228‑1 (MQ=255)
                                                                                      aCGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACCGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACCCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGGTTTCACCCATACCGCGCGTCTCGGCCTTCCGGTCGAAGGGCAGGCAGTCgggggac                                                                                                                                                         >  2:778144/1‑235 (MQ=255)
                                                                                      aCGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGCCGCCGGGGTTGGTGGCAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATCGTCGTTTTCACCCATATCGCGAGCCTCGGCCCTGCGGTCGTAGGGCAGGCGGTCgggtgt                                                                                                                                                          >  2:870525/1‑232 (MQ=255)
                                                                                      aCGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGtt                                                                                                                                                                                                                >  1:510737/1‑183 (MQ=255)
                                                                                      aCGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCCACCATATCGCGAGTCTCGGCCTTGCGGTAGCAGGGCAGGCAGTCGAgggcg                                                                                                                                                         >  2:204993/1‑235 (MQ=255)
                                                                                      aCGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTa                                                                                                                                                          >  1:255686/1‑237 (MQ=255)
                                                                                      aCGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTa                                                                                                                                                          >  1:126299/1‑237 (MQ=255)
                                                                                      aCGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGCAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCCTGAGTCTCGGCCTTGCGGCAGTAGGGAAGGCAGTCGAGGGAg                                                                                                                                                         >  2:1368973/1‑238 (MQ=255)
                                                                                                    tGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTCGTAGGGCAGGCAGTCGAGGTAGCCCCAGGTCTCCAg                                                                                                                                           >  2:834888/1‑238 (MQ=255)
                                                                                                    tGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGCAGTAGGGCAGGCAGTCGCGGTAGCTCCAGTTCTcccg                                                                                                                                           >  2:231107/1‑236 (MQ=255)
                                                                                                                    ccGGCCCAAACCCTTCCGGCTGATACCC‑GTTGAGACCTCCCGTGCGGGGGTAGCCCGCCTGCACGCCCGCTTAAATCATCGCTTAAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCC‑CCCCGGCCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGcc                                                                                                                           <  1:231107/238‑1 (MQ=255)
                                                                                                                    acgacccAAACCCTTCCTGCTGAGACCC‑GTTGAGTTCGCCGGTACGCGGGTAGGCCGCGTGCGCGCCCGCTTCAACCAGCGCTTAAACCGGCGGTTTGACTCGGGGTTGGGGGGTGGA‑CACGCT‑CCCCGGCCCA‑CCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTACTTGCGcc                                                                                                                           <  2:639827/234‑1 (MQ=255)
                                                                                                                       gACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGGTTGGCGGCAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGCGTCTCGGCCTTGCGGCAGTAGGGCAGGCAGACGAGGTAGCTCCAGGTATCCAGACCGGGCGCTTGAGCTCAg                                                                                                                        >  2:267040/1‑239 (MQ=255)
                                                                                                                                  tCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGAt                                                                                                              >  1:1428594/1‑238 (MQ=255)
                                                                                                                                  tCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGAt                                                                                                              >  2:590950/1‑238 (MQ=255)
                                                                                                                                    ctgctgATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCa                                                                                                           >  2:192142/1‑239 (MQ=255)
                                                                                                                                    ctgctgATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGATTGGAGGTAGT‑GACGCT‑CACCGCGCCA‑TCACCGCCGTGATAGTAGTTTTCACCCATATCGCGCGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTGGCTCCAGTTCTCCAGACCGGGCTCTTGCGCCCAGTTATCGAGTTCCa                                                                                                           >  2:420777/1‑239 (MQ=255)
                                                                                                                                          aTAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCAt                                                                                                     >  1:857996/1‑239 (MQ=255)
                                                                                                                                              ccgtttGAGATCGCCCGTGCGGGGGTGGCCCGCCTGCCCCCCCGCTTCAATCCTCGCTTCAAACAGCGGATTGGCCCGGGGTTTGGAGGTAGT‑GACGCT‑CCCCGGGCCA‑TCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGc                                                                                                  <  1:540565/235‑1 (MQ=255)
                                                                                                                                              cc‑GTGGGGAGCGCCCGTGCGCGGGTACCCCGCCGGCCCGCCCTCTTCAATCTCCGCTCCAAACAGCGGATGGACGCGGGGTTTGGAGGTAGT‑GACGCC‑CCCCGGGCCA‑CCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGcc                                                                                                 <  2:1115695/239‑1 (MQ=255)
                                                                                                                                               ccggtgagcggcccggggcgcGGGTGGCCCGCCGGCCCGCCCGCTTCAACCTCCGCTCCAACCGGCGGATGTGCGCCGGGTTTGGGGGTAGT‑GCCGCC‑CCCCGGCCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCa                                                                                                      <  2:1277622/217‑1 (MQ=255)
                                                                                                                                                 gTTGAGACCGCCCGTGCGGGGGTACCCCGCCTGCACGCCCGCTTCAATCATCGCTCCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCAc                                                                                               <  1:90315/239‑1 (MQ=255)
                                                                                                                                                  ttGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATCGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACg                                                                                              >  1:449395/1‑239 (MQ=255)
                                                                                                                                                  ttGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACg                                                                                              >  2:448215/1‑239 (MQ=255)
                                                                                                                                                  ttGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACg                                                                                              >  1:172109/1‑239 (MQ=255)
                                                                                                                                                          gTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGc                                                                                      >  1:472084/1‑239 (MQ=255)
                                                                                                                                                          gTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGGAGCCCCAGGTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGACACATCGCATTGCCACGGCTGTAg                                                                                       >  2:72257/1‑238 (MQ=255)
                                                                                                                                                          gTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGCTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCAACCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGTGATCAAGCGCATTGCCACGGATGTAGc                                                                                      >  1:1881/1‑239 (MQ=255)
                                                                                                                                                                   ggggggccccGCCGGCACGCCCGCTTAAATCATCCCTTCAAACAGCGGTTTGCCGCGGGTTTTGGAGGTAGT‑GACGCC‑CCCCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCg                                                                              <  1:1157737/231‑1 (MQ=255)
                                                                                                                                                                     ggTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGCAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGca                                                                                     >  1:390885/1‑229 (MQ=255)
                                                                                                                                                                        cgCCCGCCTCTACGCCCGCTGCAACCACCGCTTAAACCAGCGGATGTCCGCCGGGTTTGGAGGTAGT‑GCCGCC‑CACCGGGCCA‑TCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGAt                                                                         <  1:167452/237‑1 (MQ=255)
                                                                                                                                                                            ggcgggCACGCCGGCTTCAATCACCGCTACAAACAGCGGAGTGA‑GCGGGGGTTGGGGGGGGTGGACGCT‑CCCGGGGCAC‑CCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGcg                                                                    <  2:671407/234‑1 (MQ=255)
                                                                                                                                                                                      cccTCTTCAATCAGCGCTGCAAACAGCGGATTGGCGGGGGGGTTGGTGGTAGG‑GCCGCT‑CCCCGGGCCC‑GCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCAcc                                                          <  2:625382/239‑1 (MQ=255)
                                                                                                                                                                                      cccGTTTCAATCAGCCCTTCAAACGGGGGATGGACGCGGGGTTGGGGGTTGGT‑GACGCT‑CCCCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCAcc                                                          <  1:304021/239‑1 (MQ=255)
                                                                                                                                                                                      cccGCTTCAATCACCGCTTAAAACAGCGGTTTGACGCGGGGTTTGGGGTGAGT‑GACGCC‑CCCCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCAcc                                                          <  1:778144/239‑1 (MQ=255)
                                                                                                                                                                                       ccGCTTAAAGCAGCGCTTAAAACAGCGGATTGACGCGGGGTTGGGAGGTAGT‑GACGCC‑CACGGGGCCA‑GCGCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCAcc                                                          <  1:1096758/238‑1 (MQ=255)
                                                                                                                                                                                        ttcttaaaccctcccgaaaaaCGGCGGTGTGACGCGGGGTTTGAAGTTAG‑‑GCAGCTCCCCCGGGCCA‑CCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCa                                                        <  2:449395/221‑1 (MQ=255)
                                                                                                                                                                                        ggctcaattccgcGCTTCAAACGGCGGATGGCCGCGGGGTTGGGTGGTAGC‑GCCGCC‑CCCCGGCCCT‑TCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCa                                                        <  2:94451/227‑1 (MQ=255)
                                                                                                                                                                                        cgttccaatacacccTTAAAACGGCGGTTTGCCGCCGGTTTTGGAGGTAGT‑GCCGCT‑CACCGGCCCA‑CCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCa                                                        <  1:1195197/225‑1 (MQ=255)
                                                                                                                                                                                        cctttCAACCACCGCTTAAACGAGGGGTTTGGGGCGGGGTTGGGGGGTGGT‑CCCGCT‑C‑CCGGGCCACTCTCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCa                                                        <  2:255686/236‑1 (MQ=255)
                                                                                                                                                                                        cGCTTAAATCACCACTCAAACCAGCGATTTGGCGCGGGGTTTGGGGTGAGT‑GACGCC‑CCCCGGGCCA‑TCACCGCCGTGATAGTCATTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCa                                                        <  1:337249/239‑1 (MQ=255)
                                                                                                                                                                                         gattCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAg                                                       >  2:577196/3‑239 (MQ=255)
                                                                                                                                                                                                    cGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACcgc                                            >  1:1013665/1‑239 (MQ=255)
                                                                                                                                                                                                          aaaCAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GCCGCT‑CCCCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGc                                      <  1:1388970/239‑1 (MQ=255)
                                                                                                                                                                                                                    gttaCGCGGGGTGTGAGGCTAGT‑GCCGCC‑CTCCGGCCCC‑GCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATgc                             <  2:442371/235‑1 (MQ=255)
                                                                                                                                                                                                                          cgggttttggggagaGT‑GACCCC‑CCCCGGCCCA‑CCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGt                       <  2:782814/224‑1 (MQ=255)
                                                                                                                                                                                                                            gggTTTGGGTGTAGG‑GGCCCC‑CCCGGGCCCA‑CCCCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTa                     <  1:591853/238‑1 (MQ=255)
                                                                                                                                                                                                                            cggggtgggaggcgggaaCGCT‑CAGCGGGCCT‑ATGCCGCCGTGATAGTCGTATTCACCGATATCGCAAGTCGCGGCCTTTCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGATTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACAGATGTAGCACATGCCGTTGATCAGCGGCGATCCACCCAGGCCTTTACCGCGTCCGCACTCCATGCGGGGGTTa                     <  2:222552/222‑1 (MQ=255)
                                                                                                                                                                                                                                 tGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATa               >  1:281969/1‑239 (MQ=255)
                                                                                                                                                                                                                                     ggTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAAc            <  1:192142/238‑1 (MQ=255)
                                                                                                                                                                                                                                        aGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGtt        >  2:627365/1‑239 (MQ=255)
                                                                                                                                                                                                                                        aGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGCGATTCATAAACggct        >  2:1359521/1‑237 (MQ=255)
                                                                                                                                                                                                                                             aCGCT‑CACCGGGCGA‑TTGCCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAgg    <  1:1001109/239‑1 (MQ=255)
                                                                                                                                                                                                                                              cGCT‑CCCCGGCCCT‑CCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAgg    <  1:877363/238‑1 (MQ=255)
                                                                                                                                                                                                                                                cT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCCCGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGCCCTTTACCGCGTCAGCACTCCCTGCGGCGGTTCTTCACTACCGGTTACGGtt  >  2:433546/1‑238 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
TGTGATCGGTCATAGCGTGAGTACGAATGGTCAGGTTAGGACGCGATTTGGCCTGATCGAGATAGCCACGCGCGGT‑GCTGGCGCGACGGCCCTGCGGCGTGACGGTGCG‑ATCCATCGGACCAAAACCTTCCTGCTGATAACC‑GTTGAGATCGTCCGTGCGCGGGTAGCCCGCCTGCACGCCCGCTTCAATCATCGCTTCAAACAGCGGATTGACGCCGGGTTTGGAGGTAGT‑GACGCT‑CACCGGGCCA‑TCACCGCCGTGATAGTCGTTTTCACCCATATCGCGAGTCTCGGCCTTGCGGTAGTAGGGCAGGCAGTCGAGGTAGCTCCAGTTCTCCAGACCGGGTTCTTGCGCCCAGTTATCGAGATCCAGCGCATTGCCACGGATGTAGCACATGCCGTTGATCAGCGACGATCCACCCAGACCTTTACCGCGTCCGCACTCCATGCGGCGGTTATTCATAAACGGTTCAGGTT  >  NC_000913/326575‑327048

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