Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,600,215 +C intergenic (‑26/+23) yjjA ← / ← dnaC DUF2501 domain‑containing protein YjjA/DNA replication protein DnaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,600,2141.C97.6% 284.7 / ‑6.4 84intergenic (‑25/+24)yjjA/dnaCDUF2501 domain‑containing protein YjjA/DNA replication protein DnaC
Reads supporting (aligned to +/- strand):  ref base . (2/0);  major base C (44/36);  minor base A (2/0);  total (48/36)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGC  >  NC_000913/4599964‑4600457
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
gatgggcgtgcgtgtggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCg                                                                                                                                                                                                                                                       >  2:218521‑M2/245‑251 (MQ=255)
gatgggcgtgcgtgtggtcaccggtatcagtttcggatttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCg                                                                                                                                                                                                                                                       >  2:687876‑M2/245‑251 (MQ=255)
  ggggcgtgcgtgtggccaccggtatcagtttcggttttgtcttctacgtcctggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCgcc                                                                                                                                                                                                                                                    <  2:809761‑M2/9‑2 (MQ=255)
             gtggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGcc                                                                                                                                                                                                                                         <  1:441468‑M2/20‑1 (MQ=255)
               ggtccccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcaccccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCa                                                                                                                                                                                                                                        <  1:578853‑M2/21‑1 (MQ=255)
               ggtcaccggtatcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAg                                                                                                                                                                                                                                       <  1:202637‑M2/22‑1 (MQ=255)
                         atcagtttgggttttttcttctacgtacgggaccagttcgtcggcccgcttacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAAc                                                                                                                                                                                                                                   <  2:209043‑M2/26‑1 (MQ=255)
                         atcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAAc                                                                                                                                                                                                                                   >  1:209043‑M2/220‑245 (MQ=255)
                         atcagtttcggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGa                                                                                                                                                                                                                              <  1:687876‑M2/31‑1 (MQ=255)
                           cagtttgggttttgctttctacgtacgggaccagaccgccggcccgctgacgttggtttatggcaccccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGAtt                                                                                                                                                                                                                            <  2:534795‑M2/33‑1 (MQ=255)
                             gtttcggttttgtcttctacgtactggaccagatcttcggcccgctgaagatggtttatggcatccagccgatcataggcgcgctgttgccaagcgccagcttcttcttactcagcctgtggctgttcatgagcaaatcgtaaccgaatgaacgcataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGACAGTATACCGCCAGAAACGGGGATTAAt                                                                                                                                                                                                                         >  1:876982‑M2/216‑251 (MQ=255)
                                  ggttttgtcttctacgtactggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATActct                                                                                                                                                                                                                    >  1:202350‑M2/211‑251 (MQ=255)
                                         tcttctacgtactggaccagatcttcggcccgctgcagttggtttatggcatcccgccgatcatcggcgcactgttgccacgcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTg                                                                                                                                                                                                             >  2:473494‑M2/204‑251 (MQ=255)
                                               ccgtacgggaccagatcttcggccggctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgctgcgttcattcctgaTGTTTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTAcc                                                                                                                                                                                                        <  2:899207‑M2/53‑1 (MQ=255)
                                                cgtacgggaccagatctccggcccgctgacgtgggtttagggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCg                                                                                                                                                                                                      <  2:523172‑M2/55‑1 (MQ=255)
                                                      ggaccagatcttcggcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACg                                                                                                                                                                                                >  1:160422‑M2/191‑251 (MQ=255)
                                                        accagatcttcggccctctgccgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctgtattattgttagcaaagtgtgcttcgttcattcccgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGt                                                                                                                                                                                              <  1:17481‑M2/63‑1 (MQ=255)
                                                                    gcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGt                                                                                                                                                                                  >  1:622133‑M2/177‑251 (MQ=255)
                                                                    gcccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcagtcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGt                                                                                                                                                                                  <  1:473494‑M2/75‑1 (MQ=255)
                                                                     cccgctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGtt                                                                                                                                                                                 >  1:609782‑M2/176‑251 (MQ=255)
                                                                       cgctgacgttggtttatggcatcccgccgatcatcggcgaactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgcataaagggagcgaaaactccctttattattgttagcaacgtgtgcttcgttcattcctgtTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTg                                                                                                                                                                                >  1:552786‑M2/174‑250 (MQ=255)
                                                                         ctgacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGa                                                                                                                                                                               >  1:357160‑M2/172‑249 (MQ=255)
                                                                           gacgttggtttatggcatcccgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGAt                                                                                                                                                                           <  1:730522‑M2/82‑1 (MQ=255)
                                                                                 ggtttatggcatcccgccgataatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAccc                                                                                                                                                                      >  1:303569‑M2/164‑250 (MQ=255)
                                                                                     tatgtcatcccgccgatcatcggcgcacgtttgccaagctccagcttcttattaatcagcctgtggctgttaatgataaaatcgtaaccgaatgaacgaatacagggagcgaaagctccctttattattgttagcaaagtgttcttcgttcattcctgaTGTGTCGCCAGTTCACCGTCAGAAACGGGGATTAATACTCTTTACCTGTTTCCCGGCTACGGTAGTTAACCCTGATGAAGATCACCCACaa                                                                                                                                                                  <  2:164928‑M2/91‑1 (MQ=255)
                                                                                             ccagccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTAGACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAcc                                                                                                                                                         <  1:416779‑M2/100‑1 (MQ=255)
                                                                                               cgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCt                                                                                                                                                                                          <  1:320024‑M2/67‑1 (MQ=255)
                                                                                               cgccgatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCt                                                                                                                                                                                          >  2:320024‑M2/150‑216 (MQ=255)
                                                                                                   gatcatcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGgcg                                                                                                                                                   >  1:234483‑M2/146‑251 (MQ=255)
                                                                                                      catcggcgcactgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCAt                                                                                                                                                <  2:830872‑M2/109‑1 (MQ=255)
                                                                                                                ctgttgccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATc                                                                                                                                      >  2:549745‑M2/133‑251 (MQ=255)
                                                                                                                     gccaagcgccagcttcgtcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttatgattggtagcaaagtgtgcttcgttcattcctgaTGTGTCgcc                                                                                                                                                                                                                                                    <  1:855751‑M2/9‑2 (MQ=255)
                                                                                                                      ccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCATTATACCGCCAGAAACGTGGATTAATACTCTTTACCTGTGACCCTGCTACGGTAGCTATCCCAg                                                                                                                                                                                   <  1:360451‑M2/74‑1 (MQ=255)
                                                                                                                      ccaagcgccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAg                                                                                                                                                                                   >  2:360451‑M2/127‑200 (MQ=255)
                                                                                                                             ccagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGccc                                                                                                                          >  1:269013‑M2/120‑250 (MQ=255)
                                                                                                                              cagcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGAGTAATACTCTTTACCTGTTACCCGGCTACGGTAGGTATCCCAGTTGAAGATCACCCACAAACTGTTTCCCAGGCGCATACGGGCCATCTCGCGTTCGCCCAg                                                                                                                        <  2:479895‑M2/133‑1 (MQ=255)
                                                                                                                                gcttcttcttaatcagcctgtggctgttaatgagaaaatcgtacccgaatgaacgaattaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATCCGGTCCATCACGCGTTCGCCCAGc                                                                                                                       >  2:410783‑M2/117‑250 (MQ=255)
                                                                                                                                gcttcttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCa                                                                                                                      <  1:567463‑M2/135‑1 (MQ=255)
                                                                                                                                    cttcttaatcagcctgtggatgctaatgccaaaatcgtaaccgaatgaacgaatacagggagcgaaagctccctttattattgttagaaaagagtgcttcgttgattcctgaTGTGTAGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGCAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCTTCCAGCGttcgcccacaacctt                                                                                                                  >  2:55556‑M2/113‑244 (MQ=255)
                                                                                                                                     ttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATc                                                                                                                                                                          >  1:309137‑M2/112‑194 (MQ=255)
                                                                                                                                     ttcttaatcagcctgtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGGTACCCGGCTACGGTAGCTATCCCAGTTGAAGATc                                                                                                                                                                          <  2:309137‑M2/83‑1 (MQ=255)
                                                                                                                                           atcagccggtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgtgcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTTTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAtc                                                                                                           <  2:297668‑M2/146‑1 (MQ=255)
                                                                                                                                              agcctgtggctgttaatgagaaaatcgtaaccgaattaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgctgcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCCCCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCatctca                                                                                                        <  2:288139‑M2/149‑2 (MQ=255)
                                                                                                                                                   gtggctgttaatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCat                                                                                                   <  2:160422‑M2/154‑1 (MQ=255)
                                                                                                                                                   gtggctgttaatgagaaaagcgtaaccgaaggaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTGGAAGATCACCCACAAACTGTTTCCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTGGGTCATCTCTTCCat                                                                                                   <  2:609782‑M2/154‑1 (MQ=255)
                                                                                                                                                        tgttaatgagaaactcgtaaccgaatgaacgaataaagggagcgacagatccctttattattgttcgcaaagtgtgcttcgttcattcctgaTGTGTCGCTCGTATACCGCCAGAACCGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCCGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATCCGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTAc                                                                                              >  1:565555‑M2/93‑251 (MQ=255)
                                                                                                                                                           taatgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTg                                                                                            <  1:410783‑M2/161‑1 (MQ=255)
                                                                                                                                                              tgagaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTgg                                                                                        >  1:917035‑M2/87‑251 (MQ=255)
                                                                                                                                                                agaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCAAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTc                                                                                      >  1:815775‑M2/85‑251 (MQ=255)
                                                                                                                                                                 gaaaatcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCa                                                                                     <  1:77217‑M2/168‑1 (MQ=255)
                                                                                                                                                                    aatcgtaaccgaatgaacgaataaagggcgcgcaagctccctttattcttgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGCTTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACACACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCACCTTGGTCATCTCTTCCATATTACTGTTGGTCAGc                                                                                   >  2:366493‑M2/81‑250 (MQ=255)
                                                                                                                                                                     atcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCAt                                                                                 >  1:33421‑M2/80‑251 (MQ=255)
                                                                                                                                                                     atcgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCAt                                                                                 >  2:68626‑M2/80‑251 (MQ=255)
                                                                                                                                                                       cgtaaccgaatgaacgaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATcc                                                                               >  1:811732‑M2/78‑251 (MQ=255)
                                                                                                                                                                          aaccgaatgaacgaatacagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCgg                                                                            >  1:735830‑M2/75‑251 (MQ=255)
                                                                                                                                                                                atgaacgaataaagggagctaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGc                                                                      <  2:269013‑M2/183‑1 (MQ=255)
                                                                                                                                                                                      gaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCCGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGc                                                                      >  2:143432‑M2/63‑245 (MQ=255)
                                                                                                                                                                                      gaataaagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGc                                                                      <  1:143432‑M2/183‑1 (MQ=255)
                                                                                                                                                                                           aagggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaaga                                                            >  2:518500‑M2/58‑250 (MQ=255)
                                                                                                                                                                                            agggagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGCTCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGgaaga                                                            >  2:368829‑M2/57‑249 (MQ=255)
                                                                                                                                                                                               gagcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgc                                                       >  1:730271‑M2/54‑251 (MQ=255)
                                                                                                                                                                                                agcgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcg                                                      >  2:701026‑M2/53‑251 (MQ=255)
                                                                                                                                                                                                  cgaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggc                                                    <  1:218521‑M2/201‑1 (MQ=255)
                                                                                                                                                                                                   gaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                                   >  2:272454‑M2/50‑251 (MQ=255)
                                                                                                                                                                                                   gaaagctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcggcg                                                   >  1:355635‑M2/50‑251 (MQ=255)
                                                                                                                                                                                                      agctccctttattattgttagcaaagtgtgcttcggtcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCGCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATc                                                <  2:702367‑M2/205‑1 (MQ=255)
                                                                                                                                                                                                       gctccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGTATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcg                                               <  1:366493‑M2/206‑1 (MQ=255)
                                                                                                                                                                                                         tccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATcgac                                             >  1:338070‑M2/44‑251 (MQ=255)
                                                                                                                                                                                                          ccctttattattttgagcaaagtgtgcctctttcatttcggaTGTGTCGCCAGTATACCGCCAAAAACTGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACTCAGGCGCATACGGTCCATCACGCTTTCGCCGAGCAACTTGTTCATCTCTTCCATATTCCTGTTGGTCATCATCCCGGTTGGGCGTTTGTAAGAAGAGCGGCGATcgacg                                            <  1:414914‑M2/209‑1 (MQ=255)
                                                                                                                                                                                                          ccctttattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAccc                                                                                                                                                        >  1:32708‑M2/43‑143 (MQ=255)
                                                                                                                                                                                                          ccctgtattattgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACTGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTTccc                                                                                                                                                        <  2:32708‑M2/101‑1 (MQ=255)
                                                                                                                                                                                                                   attgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGGTCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAAATTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtg                                   >  2:568869‑M2/34‑251 (MQ=255)
                                                                                                                                                                                                                    ttgttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCAAAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtg                                   >  2:617113‑M2/33‑250 (MQ=255)
                                                                                                                                                                                                                     tgttagcacagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTCATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtgat                                 >  2:116216‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                                     tgttagcaaagtgtggttcgttcattcatgaTGTGTCGACAGTATAACGCCAGAAACGGGGCTTAATACTCTTTACCTGTTACCAGGCTACGGTAGCTCTCCCAGTTGAAGATCACCCACAAACTGTTCCCCAGGAGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCACGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGTCGGTCTGGTtgat                                 >  2:744768‑M2/32‑251 (MQ=255)
                                                                                                                                                                                                                      gttagcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTtgatg                                >  1:18283‑M2/31‑251 (MQ=255)
                                                                                                                                                                                                                        tagcaaagtgtgattggttccttcctgaTGTGTCGCCAGTATACCGCCAGAACCGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCCTATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTgatgat                              >  2:111804‑M2/29‑251 (MQ=255)
                                                                                                                                                                                                                         agcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATc                             <  2:735830‑M2/224‑1 (MQ=255)
                                                                                                                                                                                                                          gcaaagtgtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGAtcc                            >  1:430038‑M2/27‑250 (MQ=255)
                                                                                                                                                                                                                            aaattgtgcttcgtgcattccggaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTGTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGGCCAGCTCGCGTTCGCCCAGCATCTTGGTCATCTCTTCCATATTGCTGTTTGTCAGCATCTCGGTGGGGCGTTTGGGAGAAGAGCGGCCATCGACGATCTGGTTTATGATCACt                          <  2:565555‑M2/227‑1 (MQ=255)
                                                                                                                                                                                                                                 gtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGAt                                                 <  2:519074‑M2/204‑1 (MQ=255)
                                                                                                                                                                                                                                 gtgcttcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGAt                                                 >  1:519074‑M2/20‑223 (MQ=255)
                                                                                                                                                                                                                                      tcgttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTGGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTAtt                >  2:705254‑M2/15‑251 (MQ=255)
                                                                                                                                                                                                                                         ttcattcctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCg             >  2:316871‑M2/12‑251 (MQ=255)
                                                                                                                                                                                                                                               cctgaTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCgg       <  1:116216‑M2/246‑1 (MQ=255)
                                                                                                                                                                                                                                                   aTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGCCGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTg   >  1:162668‑M2/2‑251 (MQ=255)
                                                                                                                                                                                                                                                   aTGTGTCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTg   >  1:56054‑M2/2‑251 (MQ=255)
                                                                                                                                                                                                                                                     gtgtCGCCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCCCTTTTTCGTATTTCGATTCGGTCTGc  >  2:695023/1‑250 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    
CGGCATTGTTCATGTTATCTGCGCTTAAGGCCTGGTTTCCGCTGCTAAGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCGGCGATCGACGATCTGGTTGATGATCACTTTTTCGTATTTCGATTCGGTCTGC  >  NC_000913/4599964‑4600457

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