Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 9,488 G→A L61L (CTG→CTA)  mog → molybdopterin adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,4880GA97.4% 221.4 / NA 77L61L (CTG→CTA) mogmolybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (42/32);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTG‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAA  >  NC_000913/9241‑9725
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
tgtctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAACCGTTGTGTGCGCta‑gtg                                                                                                                                                                                                                                            >  1:660960/1‑251 (MQ=255)
 gtctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCta‑gtg                                                                                                                                                                                                                                            <  1:147623/250‑1 (MQ=255)
     gtcACTCCCTTCGCAGTTGATCATTCTGTTTAACGAGATTGTTTAAACGGAAAAATCTTGATTAATACTGTACGTATTGGCTTAGTTTCCATCGCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGa                                                                                                                                                                                                                                       <  2:554604/251‑1 (MQ=255)
                    gtgtATCATTCTTTTGAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACGTTACGTATTGGCTTAGTTGCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCACCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACGGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAGCTTTGTGTTAGCTT‑GTGGATGCAATGAGTTGCCctct                                                                                                                                                                                                                        <  2:27872/251‑1 (MQ=255)
                      gtATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTg                                                                                                                                                                                                                       <  1:398278/250‑1 (MQ=255)
                       tATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGt                                                                                                                                                                                                                     >  1:572631/1‑251 (MQ=255)
                          cATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTg                                                                                                                                                                                                                    <  1:12684/249‑1 (MQ=255)
                          cATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTg                                                                                                                                                                                                                    >  2:12684/1‑249 (MQ=255)
                          cATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCt                                                                                                                                                                                                                  >  1:535289/1‑251 (MQ=255)
                              cTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGCTGAAATGAGTTGCCATCTGGTGCTcacc                                                                                                                                                                                                              >  2:179175/1‑251 (MQ=255)
                                gTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATGGGCTTAGTTTCCATCGCTGATCGCGCATCCAGCGGCGGTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCTCTTAATCCCCGTTGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTcac                                                                                                                                                                                                               <  1:399990/248‑1 (MQ=255)
                                    aaCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGc                                                                                                                                                                                                        >  1:554653/1‑251 (MQ=255)
                                          aCTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGGACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACt                                                                                                                                                                                                  >  1:165686/1‑251 (MQ=255)
                                            tGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAACCGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTgg                                                                                                                                                                                                >  2:79070/1‑251 (MQ=255)
                                             gTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAg                                                                                                                                                                                                                                 <  1:140407/217‑1 (MQ=255)
                                             gTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAg                                                                                                                                                                                                                                 >  2:140407/1‑217 (MQ=255)
                                                tAAACGGAAAAAGCTGTATGAATACTGTACGTATGGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTGTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCTGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCCCCACGGGCGGAACTGGcct                                                                                                                                                                                             <  2:660960/250‑2 (MQ=255)
                                                      gAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGCTAAAGGCATCACTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAACCCCGCTTAATCCCCGATGAGCAGGCGATCATAGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGGGTTGCCATATGGTGCTCACCACGGGCGGAACTGGCCGGGcgcgg                                                                                                                                                                                      >  2:595333/1‑250 (MQ=255)
                                                      gAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGcgt                                                                                                                                                                                      >  1:377853/1‑251 (MQ=255)
                                                         caaTCTTGATGAATACTTTACGTATTGGATTCGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGCAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAACCACGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGCGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGGGCTACCCACGGGCGGAACTGGCCCGGCGcgacgt                                                                                                                                                                                   >  2:442312/2‑251 (MQ=255)
                                                                  tgaATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGCCATAGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAATCGTTGTGTGAGCTA‑GTGGTTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACg                                                                                                                                                                          >  1:563201/1‑251 (MQ=255)
                                                                    aaTACTTTACGAATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCCCGCCGTTTGAACTGGAAACCCGATTAATCCCCGATGAGCAGGCGATCATCGAGCAACAGATGTGTGAGCTA‑GTGGATGAAATGAGTTGCCACCTAGTGCTCACCACGGGCGGACCTGGCCCGGCGCGGCGGGACGTAACGcc                                                                                                                                                                        >  2:389548/1‑251 (MQ=255)
                                                                     attcattactttttgtcTTAGTTTCCAGCTCTGATCGCGCATACATTGGTGTTAATCAGAATAAAGGAATCCCTGCGCTGGAAGAATGGCTGACATCGGTGCTAACCACGCCGTTTTAACTGGAAACCTGCTTAAGCCCCGAGTATCAGGCTAGCATCGATGATACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCGCGGGTTTAATTGGTCCGGCGCGTGGTGACGTTACGccc                                                                                                                                                                       <  2:376997/235‑1 (MQ=255)
                                                                          ttACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATg                                                                                                                                                                   <  1:232035/250‑1 (MQ=255)
                                                                            aCGTAATGGCTTAGTTTCCGTCTCTGATAGCGCCTACAGCGGCGTTTATCAGGCTAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGTTGAGCAGGCGATCCTCGAGCAAAAGGGGTATGAGCTA‑GCGGATGAAATGCGTTGCCATCTGGCGCTCACCCCGGGCGGAACTGGCCCGGCGCGGCGTGCAGTCACGCCCGATGCgg                                                                                                                                                                >  2:333589/1‑250 (MQ=255)
                                                                             cGTATTGGCTTAGTTTCCATCTCTGATCGCTCATCCAGCGGCTTTTATCAGTATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAATGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGAc                                                                                                                                                               <  1:173775/251‑1 (MQ=255)
                                                                              gTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCCGGATAAAGGCATCCCTGCGCTGGAAGAATGGATGCCCTCGGCGCTAACCACACCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCCAACGTTGTGTGAGCTA‑GTGGATGAAATGCGTTGCCATCTGGTGCTCACAACGGGCGGAACTGGCCCGGCGCGTCGTTACGTCACGCCCAATGCGACg                                                                                                                                                              >  2:391535/1‑251 (MQ=255)
                                                                                aTTGGCATAGTTTACATCTCTGATCGTGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGcagagcc                                                                                                                                                            >  2:289267/1‑245 (MQ=255)
                                                                                     cTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGc                                                                                                                                                         >  1:207217/1‑249 (MQ=255)
                                                                                         gTTTCCATCTCTGATCGCGAATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGCAACGCCCGATGCGACGCTggacgtgac                                                                                                                                                   >  2:463704/1‑244 (MQ=255)
                                                                                                 ctctGATCGCGCCTCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGCGCAAACGTTGTGTGAGCTA‑ATGGATGAAATGAGTTGCCATCTGGTGCTAACCACGGGCGGGACTGGCCCGGCGCGTCGTGACGTGACGCCCGATGCGACGCTGGCAGTAGCGGACcgcg                                                                                                                                           >  2:587733/1‑251 (MQ=255)
                                                                                                               tcAGCGGCGTTTATCAGGATAAAGGCACCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGACCTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTcaccac                                                                                                                                                                                                            <  2:231063/171‑1 (MQ=255)
                                                                                                               ccAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTcaccac                                                                                                                                                                                                            >  1:231063/1‑172 (MQ=255)
                                                                                                                    ggcgTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGACACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGCAGTGCCGTAGCGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTGTggggg                                                                                                                        >  2:243531/1‑248 (MQ=255)
                                                                                                                     gcgTTTAGCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACGGGCCCGGCGCGTCGGGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGaa                                                                                                                       <  1:69604/251‑1 (MQ=255)
                                                                                                                       gTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAAc                                                                                                                      >  1:367323/1‑250 (MQ=255)
                                                                                                                          tATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGACGTGACGTAACGCCCGCTGCGACGATGGCAGTAGCGGACCGCGGGAGGCCTGGCTTTGGTGAAccgtt                                                                                                                  >  2:89798/1‑247 (MQ=255)
                                                                                                                                 aTAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAg                                                                                                           <  2:554653/251‑1 (MQ=255)
                                                                                                                                   aaaGGCAACCCTGCGCTGGAAGAATGGCTGCAATCGTCGCCAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCAGCGGGTGCTCACCCCGGGCGGAACTGGCCCGGCGCGGCGTGACGTAACGCCCGATGTGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCTGAt                                                                                                         <  2:535289/251‑1 (MQ=255)
                                                                                                                                       gCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGACGGCGATCATCGAGAAAACGTTGTGGGCGCTA‑GTGGATGAAATGAGTTGCCATCTGGGGCTCACCAAGGGCGGAACTGAACCGGCGCGTCGTGACGTAACGCCCGATGCGACGATGGCAGTAGCGGTCCGGTAGATGCCTGGCTTTGGGGAACAGATGCGCCGGAta                                                                                                        >  1:572706/1‑247 (MQ=255)
                                                                                                                                        cATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGAGGACCGCGAGATGCGTGGCTTTGGTGCACAGATGCGCCAGATAAGcc                                                                                                    >  2:228040/1‑251 (MQ=255)
                                                                                                                                         aTCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTTGTGAACAGATGCGCCAGATAAGccg                                                                                                   >  2:508426/1‑250 (MQ=255)
                                                                                                                                            ccTGCGCTGGAACAATGGCTGACATAGGCGCTAACCCCGCCGTTTGTACTGGCAACCCGCTTAATCCCCGCTGAGCAGGCGCTCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATATGTTCCCATCTGGTGCTCACCACGGGCGGAACTGGCCAGGCGCTGCGTGACGTAACGCCCGCGGCCACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGCGAACAgctgggcaagagaagcctgc                                                                                                 >  2:476921/1‑231 (MQ=255)
                                                                                                                                                 gcTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCAttttt                                                                                           >  1:342857/1‑250 (MQ=255)
                                                                                                                                                 cccttccgattGGCTCACATCGGCGCTAACCACGCCTTTTGAACTTGAAACCCGCTTAATCCCCGAGGAGCAGTCGATCATCGAGCAAACGGTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCGCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCTCCAGATAAGCCTGCATTTTg                                                                                           <  2:274821/241‑1 (MQ=255)
                                                                                                                                                    ggAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTcacc                                                                                                                                                                                                              <  1:619327/133‑1 (MQ=255)
                                                                                                                                                    ggAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTaacc                                                                                                                                                                                                              >  2:619327/1‑133 (MQ=255)
                                                                                                                                                             gCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTg                                                                                           >  1:274821/1‑239 (MQ=255)
                                                                                                                                                               tGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTCATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGCGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCTTCGTGACGTAACGCCCGATGCGACGCTGGCCGTAGCGGCCCGCGGGTGGCCTGGCTTTGGGGAACAGATGCGCCAgttaggcc                                                                                                    >  2:329285/1‑221 (MQ=255)
                                                                                                                                                               tGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGcc                                                                                                    <  1:329285/228‑1 (MQ=255)
                                                                                                                                                                gACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATcc                                                                            <  1:463704/251‑1 (MQ=255)
                                                                                                                                                                gACACCGGCGCTATCCACGTCGTTACAGATGGAAACCCGCTTAATCCCCGATGAGCATGCGATCATCTTGCAAACGTTTTGGGAGTTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCGCGGGCGGAACTGGCCCGGCGCGTTGTGTCGTATCGCGCGATGCGACGCTGGCTTTAGCTGACCGCCAGATTCCTGGCGTAGGTGAACGGATGCACCATATAAGCCTGCATTTTGTTCCAACTGCGAt                                                                              <  2:552968/249‑1 (MQ=255)
                                                                                                                                                                  cATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCACGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGGAGCGGACCGCGGGATGCCTGGCTTTGGTGAACAGATGCGCCCGATAAGACGGCATTTTGTCCCAACTGCGATCCtt                                                                          >  2:207917/1‑251 (MQ=255)
                                                                                                                                                                     gggtgcTAACCTCGCCGGTTGAACGGGAACCCCGCTTAATTTCCGAGGAGCAGGCGATCATCGAGCAAACGTTGGGGGAGCTA‑GGGGATGAAATGATGTGCAATCGGGCGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACTCCCGATGCGACGCTGGCAGTTGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGAGGCGCCAGATAAGCCTGCATTTTGTTCCAACTGCGTTCCttt                                                                         <  2:245356/245‑1 (MQ=255)
                                                                                                                                                                       gcgcTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGGAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTcgcg                                                                     <  1:333589/251‑1 (MQ=255)
                                                                                                                                                                            aaCCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTTGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCa                                                                  <  1:416798/249‑1 (MQ=255)
                                                                                                                                                                            aaCCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAg                                                                 <  1:441747/250‑1 (MQ=255)
                                                                                                                                                                             aCCACGCCGTTTGAACTGTAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGGTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGAtaag                                                                                                      <  1:380505/212‑4 (MQ=255)
                                                                                                                                                                             aCCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCAAGCTAAGCCTGAATTTTGTACCAACTGCGATCCTTTCGCGTCAggg                                                               >  1:253074/1‑250 (MQ=255)
                                                                                                                                                                               taCGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGAGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGAGACGAAACACCCGATGCGACGCTGGCGGTAACGGAACGCGACATGCCTGGCTTTGGTGAACAGATGCGCAAGATAAGCCTGCATTTTGTACCACATACGTTCCTTTCGCGTAAGGt                                                               >  2:553062/2‑249 (MQ=255)
                                                                                                                                                                                   ccggttGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGTCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCTTCAGGttggtgt                                                         <  1:139313/247‑7 (MQ=255)
                                                                                                                                                                                         gcattggcaccccgtgtAATCCCCGTTGAGCAGGCGCTCATCGAGCAATCGTTGTGTG‑GGTGAGTGGATGAAATGAGTTGCCATCTGGTGCGTACCAGTGGCGGAACTGGCCCGGCTCGGCGTGACGTTACGACCGATGCGCCGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTGTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTTCCAACTGCGATCCTTTCGCGTCAGTTTGGTGTGATTc                                                    <  2:377853/234‑1 (MQ=255)
                                                                                                                                                                                                    ccGCTTAATCACCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGTTGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGgcgc                                        >  1:284981/1‑251 (MQ=255)
                                                                                                                                                                                                         tAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGAt                                    <  1:389548/250‑1 (MQ=255)
                                                                                                                                                                                                             ccccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCtt                                <  1:290201/250‑1 (MQ=255)
                                                                                                                                                                                                              cccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTaa                              <  1:573022/251‑1 (MQ=255)
                                                                                                                                                                                                              cccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTa                               >  1:384750/1‑250 (MQ=255)
                                                                                                                                                                                                               ccGATGAGCGGGCGATCACCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCAGCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAAAGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAAc                             <  2:105066/251‑1 (MQ=255)
                                                                                                                                                                                                                 gATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACtt                           >  1:635562/1‑251 (MQ=255)
                                                                                                                                                                                                                       cAGGCGATAATCGAGCAAACGTTATGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGAGAGAGGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGGTGGTGTGATTCGCAAACGGGACATGATACTTACCTTACCCgg                     >  2:47024/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCGATCATCGAGCAAACGTTGTGTGAGCTG‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGGTACGTAACGCCCGATGCGACGCTGGCAGTTGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCCTAACTTTCCCGGTCa                  <  2:661725/251‑1 (MQ=255)
                                                                                                                                                                                                                            gATCATCGAGCAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCACCAGGTAAGCCTGCATTTTGTACCAACTTCGGTCCTTTCGCGGCAGGGTGGTGTGATTCGCAAACAGGCGCTGCTCCTTAACTTACCCGGTAAGc                >  2:516798/1‑251 (MQ=255)
                                                                                                                                                                                                                              tcatcgggcacccgatgtgtgAACTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGACGGAACTGGCCCGGCGCGACGTGACGTAACGGCCGATGCGACGCTGGCAGTAGCGGCCCGCGGGATGCCTGGCATTGGTGAACAGGTGCGCCAGAAAAGGCTGCGATTTGGACCAACTGAGAGCCTTTTACGTCAGGATGGTGTGGTTAGCCAAAAGGCCCTGATCCTTCACTTAACCGGTCAGCCg              >  2:98341/16‑251 (MQ=255)
                                                                                                                                                                                                                                   gAGCAAACGTTGTGGGAGCTA‑GGGGATGAAATGAGTTGCCATCTGGTGCTCACCGCGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTc         <  2:165686/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAAACGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAt      >  1:231068/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aaaCGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAtt     <  1:508426/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaCGTTGTGTGAGCTA‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAtt     >  1:17926/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aaaCGTTGTGGGAGCTA‑GTGGATGAAATGAGTGGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAtt     <  1:566057/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTTGTGTGAGCTA‑GTGGATGAAATGAGTAGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTaaa  <  1:587733/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
TGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTG‑GTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAA  >  NC_000913/9241‑9725

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