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→CTAmogmolybdopterin 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

TGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAA  >  NC_000913/9241‑9725
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
tgtctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAACCGTTGTGTGCGCtagtg                                                                                                                                                                                                                                            >  1:660960/1‑251 (MQ=255)
 gtctgtcACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCtagtg                                                                                                                                                                                                                                            <  1:147623/250‑1 (MQ=255)
     gtcACTCCCTTCGCAGTTGATCATTCTGTTTAACGAGATTGTTTAAACGGAAAAATCTTGATTAATACTGTACGTATTGGCTTAGTTTCCATCGCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGa                                                                                                                                                                                                                                       <  2:554604/251‑1 (MQ=255)
                    gtgtATCATTCTTTTGAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACGTTACGTATTGGCTTAGTTGCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCACCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACGGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAGCTTTGTGTTAGCTTGTGGATGCAATGAGTTGCCctct                                                                                                                                                                                                                        <  2:27872/251‑1 (MQ=255)
                      gtATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTg                                                                                                                                                                                                                       <  1:398278/250‑1 (MQ=255)
                       tATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGt                                                                                                                                                                                                                     >  1:572631/1‑251 (MQ=255)
                          cATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTg                                                                                                                                                                                                                    <  1:12684/249‑1 (MQ=255)
                          cATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTg                                                                                                                                                                                                                    >  2:12684/1‑249 (MQ=255)
                          cATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCt                                                                                                                                                                                                                  >  1:535289/1‑251 (MQ=255)
                              cTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGCTGAAATGAGTTGCCATCTGGTGCTcacc                                                                                                                                                                                                              >  2:179175/1‑251 (MQ=255)
                                gTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATGGGCTTAGTTTCCATCGCTGATCGCGCATCCAGCGGCGGTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCTCTTAATCCCCGTTGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTcac                                                                                                                                                                                                               <  1:399990/248‑1 (MQ=255)
                                    aaCGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGc                                                                                                                                                                                                        >  1:554653/1‑251 (MQ=255)
                                          aCTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGGACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACt                                                                                                                                                                                                  >  1:165686/1‑251 (MQ=255)
                                            tGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAACCGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTgg                                                                                                                                                                                                >  2:79070/1‑251 (MQ=255)
                                             gTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAg                                                                                                                                                                                                                                 <  1:140407/217‑1 (MQ=255)
                                             gTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAg                                                                                                                                                                                                                                 >  2:140407/1‑217 (MQ=255)
                                                tAAACGGAAAAAGCTGTATGAATACTGTACGTATGGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTGTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCTGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCCCCACGGGCGGAACTGGcct                                                                                                                                                                                             <  2:660960/250‑2 (MQ=255)
                                                      gAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGCTAAAGGCATCACTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAACCCCGCTTAATCCCCGATGAGCAGGCGATCATAGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGGGTTGCCATATGGTGCTCACCACGGGCGGAACTGGCCGGGcgcgg                                                                                                                                                                                      >  2:595333/1‑250 (MQ=255)
                                                      gAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGcgt                                                                                                                                                                                      >  1:377853/1‑251 (MQ=255)
                                                         caaTCTTGATGAATACTTTACGTATTGGATTCGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGCAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAACCACGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGCGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGGGCTACCCACGGGCGGAACTGGCCCGGCGcgacgt                                                                                                                                                                                   >  2:442312/2‑251 (MQ=255)
                                                                  tgaATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGCCATAGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAATCGTTGTGTGAGCTAGTGGTTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACg                                                                                                                                                                          >  1:563201/1‑251 (MQ=255)
                                                                    aaTACTTTACGAATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCCCGCCGTTTGAACTGGAAACCCGATTAATCCCCGATGAGCAGGCGATCATCGAGCAACAGATGTGTGAGCTAGTGGATGAAATGAGTTGCCACCTAGTGCTCACCACGGGCGGACCTGGCCCGGCGCGGCGGGACGTAACGcc                                                                                                                                                                        >  2:389548/1‑251 (MQ=255)
                                                                     attcattactttttgtcTTAGTTTCCAGCTCTGATCGCGCATACATTGGTGTTAATCAGAATAAAGGAATCCCTGCGCTGGAAGAATGGCTGACATCGGTGCTAACCACGCCGTTTTAACTGGAAACCTGCTTAAGCCCCGAGTATCAGGCTAGCATCGATGATACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCGCGGGTTTAATTGGTCCGGCGCGTGGTGACGTTACGccc                                                                                                                                                                       <  2:376997/235‑1 (MQ=255)
                                                                          ttACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATg                                                                                                                                                                   <  1:232035/250‑1 (MQ=255)
                                                                            aCGTAATGGCTTAGTTTCCGTCTCTGATAGCGCCTACAGCGGCGTTTATCAGGCTAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGTTGAGCAGGCGATCCTCGAGCAAAAGGGGTATGAGCTAGCGGATGAAATGCGTTGCCATCTGGCGCTCACCCCGGGCGGAACTGGCCCGGCGCGGCGTGCAGTCACGCCCGATGCgg                                                                                                                                                                >  2:333589/1‑250 (MQ=255)
                                                                             cGTATTGGCTTAGTTTCCATCTCTGATCGCTCATCCAGCGGCTTTTATCAGTATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAATGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGAc                                                                                                                                                               <  1:173775/251‑1 (MQ=255)
                                                                              gTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCCGGATAAAGGCATCCCTGCGCTGGAAGAATGGATGCCCTCGGCGCTAACCACACCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCCAACGTTGTGTGAGCTAGTGGATGAAATGCGTTGCCATCTGGTGCTCACAACGGGCGGAACTGGCCCGGCGCGTCGTTACGTCACGCCCAATGCGACg                                                                                                                                                              >  2:391535/1‑251 (MQ=255)
                                                                                aTTGGCATAGTTTACATCTCTGATCGTGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGcagagcc                                                                                                                                                            >  2:289267/1‑245 (MQ=255)
                                                                                     cTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGc                                                                                                                                                         >  1:207217/1‑249 (MQ=255)
                                                                                         gTTTCCATCTCTGATCGCGAATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGCAACGCCCGATGCGACGCTggacgtgac                                                                                                                                                   >  2:463704/1‑244 (MQ=255)
                                                                                                 ctctGATCGCGCCTCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGCGCAAACGTTGTGTGAGCTAATGGATGAAATGAGTTGCCATCTGGTGCTAACCACGGGCGGGACTGGCCCGGCGCGTCGTGACGTGACGCCCGATGCGACGCTGGCAGTAGCGGACcgcg                                                                                                                                           >  2:587733/1‑251 (MQ=255)
                                                                                                               tcAGCGGCGTTTATCAGGATAAAGGCACCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGACCTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTcaccac                                                                                                                                                                                                            <  2:231063/171‑1 (MQ=255)
                                                                                                               ccAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTcaccac                                                                                                                                                                                                            >  1:231063/1‑172 (MQ=255)
                                                                                                                    ggcgTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGACACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGCAGTGCCGTAGCGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTGTggggg                                                                                                                        >  2:243531/1‑248 (MQ=255)
                                                                                                                     gcgTTTAGCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACGGGCCCGGCGCGTCGGGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGaa                                                                                                                       <  1:69604/251‑1 (MQ=255)
                                                                                                                       gTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAAc                                                                                                                      >  1:367323/1‑250 (MQ=255)
                                                                                                                          tATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGACGTGACGTAACGCCCGCTGCGACGATGGCAGTAGCGGACCGCGGGAGGCCTGGCTTTGGTGAAccgtt                                                                                                                  >  2:89798/1‑247 (MQ=255)
                                                                                                                                 aTAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAg                                                                                                           <  2:554653/251‑1 (MQ=255)
                                                                                                                                   aaaGGCAACCCTGCGCTGGAAGAATGGCTGCAATCGTCGCCAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCAGCGGGTGCTCACCCCGGGCGGAACTGGCCCGGCGCGGCGTGACGTAACGCCCGATGTGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCTGAt                                                                                                         <  2:535289/251‑1 (MQ=255)
                                                                                                                                       gCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGACGGCGATCATCGAGAAAACGTTGTGGGCGCTAGTGGATGAAATGAGTTGCCATCTGGGGCTCACCAAGGGCGGAACTGAACCGGCGCGTCGTGACGTAACGCCCGATGCGACGATGGCAGTAGCGGTCCGGTAGATGCCTGGCTTTGGGGAACAGATGCGCCGGAta                                                                                                        >  1:572706/1‑247 (MQ=255)
                                                                                                                                        cATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGAGGACCGCGAGATGCGTGGCTTTGGTGCACAGATGCGCCAGATAAGcc                                                                                                    >  2:228040/1‑251 (MQ=255)
                                                                                                                                         aTCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTTGTGAACAGATGCGCCAGATAAGccg                                                                                                   >  2:508426/1‑250 (MQ=255)
                                                                                                                                            ccTGCGCTGGAACAATGGCTGACATAGGCGCTAACCCCGCCGTTTGTACTGGCAACCCGCTTAATCCCCGCTGAGCAGGCGCTCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATATGTTCCCATCTGGTGCTCACCACGGGCGGAACTGGCCAGGCGCTGCGTGACGTAACGCCCGCGGCCACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGCGAACAgctgggcaagagaagcctgc                                                                                                 >  2:476921/1‑231 (MQ=255)
                                                                                                                                                 gcTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCAttttt                                                                                           >  1:342857/1‑250 (MQ=255)
                                                                                                                                                 cccttccgattGGCTCACATCGGCGCTAACCACGCCTTTTGAACTTGAAACCCGCTTAATCCCCGAGGAGCAGTCGATCATCGAGCAAACGGTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCGCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCTCCAGATAAGCCTGCATTTTg                                                                                           <  2:274821/241‑1 (MQ=255)
                                                                                                                                                    ggAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTcacc                                                                                                                                                                                                              <  1:619327/133‑1 (MQ=255)
                                                                                                                                                    ggAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTaacc                                                                                                                                                                                                              >  2:619327/1‑133 (MQ=255)
                                                                                                                                                             gCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGCTGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTg                                                                                           >  1:274821/1‑239 (MQ=255)
                                                                                                                                                               tGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTCATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGCGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCTTCGTGACGTAACGCCCGATGCGACGCTGGCCGTAGCGGCCCGCGGGTGGCCTGGCTTTGGGGAACAGATGCGCCAgttaggcc                                                                                                    >  2:329285/1‑221 (MQ=255)
                                                                                                                                                               tGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGcc                                                                                                    <  1:329285/228‑1 (MQ=255)
                                                                                                                                                                gACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATcc                                                                            <  1:463704/251‑1 (MQ=255)
                                                                                                                                                                gACACCGGCGCTATCCACGTCGTTACAGATGGAAACCCGCTTAATCCCCGATGAGCATGCGATCATCTTGCAAACGTTTTGGGAGTTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCGCGGGCGGAACTGGCCCGGCGCGTTGTGTCGTATCGCGCGATGCGACGCTGGCTTTAGCTGACCGCCAGATTCCTGGCGTAGGTGAACGGATGCACCATATAAGCCTGCATTTTGTTCCAACTGCGAt                                                                              <  2:552968/249‑1 (MQ=255)
                                                                                                                                                                  cATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCACGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGGAGCGGACCGCGGGATGCCTGGCTTTGGTGAACAGATGCGCCCGATAAGACGGCATTTTGTCCCAACTGCGATCCtt                                                                          >  2:207917/1‑251 (MQ=255)
                                                                                                                                                                     gggtgcTAACCTCGCCGGTTGAACGGGAACCCCGCTTAATTTCCGAGGAGCAGGCGATCATCGAGCAAACGTTGGGGGAGCTAGGGGATGAAATGATGTGCAATCGGGCGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACTCCCGATGCGACGCTGGCAGTTGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGAGGCGCCAGATAAGCCTGCATTTTGTTCCAACTGCGTTCCttt                                                                         <  2:245356/245‑1 (MQ=255)
                                                                                                                                                                       gcgcTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGGAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTcgcg                                                                     <  1:333589/251‑1 (MQ=255)
                                                                                                                                                                            aaCCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTTGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCa                                                                  <  1:416798/249‑1 (MQ=255)
                                                                                                                                                                            aaCCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAg                                                                 <  1:441747/250‑1 (MQ=255)
                                                                                                                                                                             aCCACGCCGTTTGAACTGTAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGGTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGAtaag                                                                                                      <  1:380505/212‑4 (MQ=255)
                                                                                                                                                                             aCCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCAAGCTAAGCCTGAATTTTGTACCAACTGCGATCCTTTCGCGTCAggg                                                               >  1:253074/1‑250 (MQ=255)
                                                                                                                                                                               taCGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGAGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGAGACGAAACACCCGATGCGACGCTGGCGGTAACGGAACGCGACATGCCTGGCTTTGGTGAACAGATGCGCAAGATAAGCCTGCATTTTGTACCACATACGTTCCTTTCGCGTAAGGt                                                               >  2:553062/2‑249 (MQ=255)
                                                                                                                                                                                   ccggttGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGTCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCTTCAGGttggtgt                                                         <  1:139313/247‑7 (MQ=255)
                                                                                                                                                                                         gcattggcaccccgtgtAATCCCCGTTGAGCAGGCGCTCATCGAGCAATCGTTGTGTGGGTGAGTGGATGAAATGAGTTGCCATCTGGTGCGTACCAGTGGCGGAACTGGCCCGGCTCGGCGTGACGTTACGACCGATGCGCCGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTGTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTTCCAACTGCGATCCTTTCGCGTCAGTTTGGTGTGATTc                                                    <  2:377853/234‑1 (MQ=255)
                                                                                                                                                                                                    ccGCTTAATCACCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGTTGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGgcgc                                        >  1:284981/1‑251 (MQ=255)
                                                                                                                                                                                                         tAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGAt                                    <  1:389548/250‑1 (MQ=255)
                                                                                                                                                                                                             ccccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCtt                                <  1:290201/250‑1 (MQ=255)
                                                                                                                                                                                                              cccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTaa                              <  1:573022/251‑1 (MQ=255)
                                                                                                                                                                                                              cccGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTa                               >  1:384750/1‑250 (MQ=255)
                                                                                                                                                                                                               ccGATGAGCGGGCGATCACCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCAGCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAAAGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAAc                             <  2:105066/251‑1 (MQ=255)
                                                                                                                                                                                                                 gATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACtt                           >  1:635562/1‑251 (MQ=255)
                                                                                                                                                                                                                       cAGGCGATAATCGAGCAAACGTTATGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGAGAGAGGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGGTGGTGTGATTCGCAAACGGGACATGATACTTACCTTACCCgg                     >  2:47024/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGGTACGTAACGCCCGATGCGACGCTGGCAGTTGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCCTAACTTTCCCGGTCa                  <  2:661725/251‑1 (MQ=255)
                                                                                                                                                                                                                            gATCATCGAGCAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCACCAGGTAAGCCTGCATTTTGTACCAACTTCGGTCCTTTCGCGGCAGGGTGGTGTGATTCGCAAACAGGCGCTGCTCCTTAACTTACCCGGTAAGc                >  2:516798/1‑251 (MQ=255)
                                                                                                                                                                                                                              tcatcgggcacccgatgtgtgAACTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGACGGAACTGGCCCGGCGCGACGTGACGTAACGGCCGATGCGACGCTGGCAGTAGCGGCCCGCGGGATGCCTGGCATTGGTGAACAGGTGCGCCAGAAAAGGCTGCGATTTGGACCAACTGAGAGCCTTTTACGTCAGGATGGTGTGGTTAGCCAAAAGGCCCTGATCCTTCACTTAACCGGTCAGCCg              >  2:98341/16‑251 (MQ=255)
                                                                                                                                                                                                                                   gAGCAAACGTTGTGGGAGCTAGGGGATGAAATGAGTTGCCATCTGGTGCTCACCGCGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTc         <  2:165686/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cAAACGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAt      >  1:231068/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aaaCGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAtt     <  1:508426/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaCGTTGTGTGAGCTAGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAtt     >  1:17926/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aaaCGTTGTGGGAGCTAGTGGATGAAATGAGTGGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTAtt     <  1:566057/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTTGTGTGAGCTAGTGGATGAAATGAGTAGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTaaa  <  1:587733/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
TGTCTGTCACTCCCTTCGCAGTGTATCATTCTGTTTAACGAGACTGTTTAAACGGAAAAATCTTGATGAATACTTTACGTATTGGCTTAGTTTCCATCTCTGATCGCGCATCCAGCGGCGTTTATCAGGATAAAGGCATCCCTGCGCTGGAAGAATGGCTGACATCGGCGCTAACCACGCCGTTTGAACTGGAAACCCGCTTAATCCCCGATGAGCAGGCGATCATCGAGCAAACGTTGTGTGAGCTGGTGGATGAAATGAGTTGCCATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTTTGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAA  >  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: