Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 9,754 A→C E150A (GAG→GCG)  mog → molybdopterin adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009139,7540AC98.7% 238.4 / NA 80E150A (GAG→GCG) mogmolybdopterin adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base C (28/49);  minor base T (0/1);  total (28/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATA  >  NC_000913/9507‑9999
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
cATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCggg                                                                                                                                                                                                                                                    >  1:539025/1‑251 (MQ=255)
 atcgggTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCGGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGt                                                                                                                                                                                                                                                   <  2:378128/247‑1 (MQ=255)
 aTCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCCGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGATGCGGGt                                                                                                                                                                                                                                                   >  2:354974/1‑251 (MQ=255)
     ggTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAAc                                                                                                                                                                                                                                                <  1:192175/250‑1 (MQ=255)
     ggTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAAc                                                                                                                                                                                                                                                <  1:520273/250‑1 (MQ=255)
     ggTGCTCACCACGGGCGGAACGGGCCCGGCGCGTCGGGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCGGGAAGGTGTGAAGGACGCTGCGGGTAAc                                                                                                                                                                                                                                                <  2:312779/250‑1 (MQ=255)
      gTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGt                                                                                                                                                                                                                                              >  1:41754/1‑251 (MQ=255)
          tCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTtgtg                                                                                                                                                                                                                                          >  1:351445/1‑251 (MQ=255)
           cacctgggtcGGAACTGGCCTGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGGAGCGGCTCGCGAGATGCCTGGCTG‑TGGTGATCAGATGCGCCAGATAAGCCTGCATTTTGGACCAACTGCGAGCCTGTCGCGTCATGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCATCCGAAGTCTACTAAATAGACGCTGGAAGGTGGGAAGGATGCTGCGGGTAACGTtgt                                                                                                                                                                                                                                           <  2:618164/240‑1 (MQ=255)
             ccacGGGCGGAACTGGCCCGACGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGt                                                                                                                                                                                                                                        >  1:334288/1‑250 (MQ=255)
               ccGGGCGGAACTGGCCCGGCGCGTCGGGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAATCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCGTATCCTTAACTGACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTAGGGGTAACGTTGTGGTaca                                                                                                                                                                                                                                     <  2:56757/250‑1 (MQ=255)
                cGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTacac                                                                                                                                                                                                                                    >  1:366398/1‑251 (MQ=255)
                           ggcgccgtgcgTCGTGAAGTAACGCCCGATGCGACGTTGGCAGTCGCGTACCGCGCTATGCCTTGCTA‑TTGTGAACAGATGCGCCAGATACTCCTGCATTTTGTACCAACTGCGATCCTTTGGCGTCAGGGGGGTGTGATTCGCAAAGAGGCGCTGATCCGTAATGTACCCGGTCAGCTTAAGTCTAGTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTCCACGGTAttttt                                                                                                                                                                                                                           <  2:556853/241‑1 (MQ=255)
                                ggCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcg                                                                                                                                                                                                                    >  1:507486/1‑251 (MQ=255)
                                    cgtcgtGACGTAACGCCCGATGCGAGGCTGGCAGTAGCTGACCGCGAGATGCTGGGCTTGTGG‑GAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGGGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTATAGAGACGCTGGAAGTTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtacc                                                                                                                                                                                                                <  2:579811/251‑1 (MQ=255)
                                     gtcgtGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtaccg                                                                                                                                                                                                               >  1:26212/1‑251 (MQ=255)
                                          gacttAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCATGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATGCGCAAACAGGCGCTGATCCTTAACGTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTGGCCAGCGTACCGTACt                                                                                                                                                                                                           <  2:142216/246‑1 (MQ=255)
                                            cGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGGACCGTACTGCa                                                                                                                                                                                                        <  1:147045/251‑1 (MQ=255)
                                                       aTGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACGGCATTCAGTTGCTg                                                                                                                                                                                             <  2:100363/251‑1 (MQ=255)
                                                           gACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGGCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAg                                                                                                                                                                                         <  1:85740/251‑1 (MQ=255)
                                                                   cAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATAc                                                                                                                                                                                 <  1:257675/251‑1 (MQ=255)
                                                                       agcggacatctagaTGCCTGGCTT‑TGGTGAACAGGTGCGCCAGATAAGCCTGCATTTTGTCCCTACTGCGATCCTTGCGTGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCGTAACTTGCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGGGGTACACGGTATTTTTGCCAGCGTACCGTACTGCAGTCAGTTGCTGGAAGGGCCATACGTTg                                                                                                                                                                             <  2:398220/240‑1 (MQ=255)
                                                                        gCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGa                                                                                                                                                                            <  1:563832/251‑1 (MQ=255)
                                                                         cGGACCGCGAGATGCCAGGCTT‑TGGTGAACAGATGCGCCAGATAACCCCGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGATGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGCATTTTTGCCAGCGTACCGTACTGCATTCCGTTGCTGGAAGGGCCTTACGTTGa                                                                                                                                                                            >  2:436489/1‑250 (MQ=255)
                                                                                    aTGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCAccaga                                                                                                                                                                <  1:503500/251‑4 (MQ=255)
                                                                                     tGCCTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCAccaga                                                                                                                                                                >  2:320665/1‑247 (MQ=255)
                                                                                        cTGGCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCCGAAGTg                                                                                                                                                            >  1:384890/1‑251 (MQ=255)
                                                                                          ggCTT‑TGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGt                                                                                                                                                          >  1:483911/1‑251 (MQ=255)
                                                                                           gCTT‑TGGTGAACCGATGCGCCCGATAAGCCTGCATTTTTTACCCACTGCTATCCTTTCTCGTCAGGTTGGTGTGATTCGCCAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTTAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCTTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGtt                                                                                                                                                         <  2:483297/251‑1 (MQ=255)
                                                                                             tt‑tGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGACGTGGTTgc                                                                                                                                                       >  2:103021/1‑251 (MQ=255)
                                                                                             tt‑tGGTGAACAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTgc                                                                                                                                                       >  1:455944/1‑251 (MQ=255)
                                                                                                       cAGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAg                                                                                                                                              <  1:103021/251‑1 (MQ=255)
                                                                                                       cAGATGCGCCAGATAAGCCTGAATGTTGGGCCAACTGCGGTCCTTTCGCGTCAGGTTGGTGGGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACAGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCGGCATTCAg                                                                                                                                              <  2:170262/251‑1 (MQ=255)
                                                                                                        aGATGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGa                                                                                                                                             <  2:507486/251‑1 (MQ=255)
                                                                                                          aTGCGCCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGAcc                                                                                                                                           >  2:102524/1‑251 (MQ=255)
                                                                                                            gcgcCAGAGAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCg                                                                                                                                          <  2:41754/250‑1 (MQ=255)
                                                                                                              gcCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtacc                                                                                                                                                                                                                <  2:31435/178‑1 (MQ=255)
                                                                                                              gcCAGATAAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGcgtacc                                                                                                                                                                                                                >  1:31435/1‑178 (MQ=255)
                                                                                                                cAGATAAGCCTGCAAGTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGGGTGAGTCGCAAACAGGCGCTGATTCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAATGTTGTGGTACACGGTAGTTTTGCCAGCGTACCGTACTGCATGCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGTTTGCTGCATTCAGACCGAagag                                                                                                                                     <  2:30689/251‑1 (MQ=255)
                                                                                                                   ataaGCCTGCATTTTGTTCCATCTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCTGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACTTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGc                                                                                                                                  <  1:420235/248‑1 (MQ=255)
                                                                                                                     aaGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCa                                                                                                                                 >  1:507196/2‑250 (MQ=255)
                                                                                                                        ccTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAc                                                                                                                             <  2:1314/251‑1 (MQ=255)
                                                                                                                         cTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAcg                                                                                                                            >  1:312972/1‑251 (MQ=255)
                                                                                                                           gCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTTCCCGG‑CAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGa                                                                                                                         <  1:348320/251‑1 (MQ=255)
                                                                                                                           gCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGAcg                                                                                                                            <  2:460749/249‑1 (MQ=255)
                                                                                                                                 tGTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGtt                                                                                                                     <  1:353149/250‑1 (MQ=255)
                                                                                                                                  gTACCAACTGCGATCCTTTCGCGTCAGGTTGGTGTGCTTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGAGGTTGCTGCATTCAGACCGCAGAGTGCAAGACGCGACGtt                                                                                                                     >  2:581030/1‑249 (MQ=255)
                                                                                                                                           gCGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaa                                                                                                           <  1:26291/250‑1 (MQ=255)
                                                                                                                                            cGATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTGACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaaaa                                                                                                         <  1:480779/251‑1 (MQ=255)
                                                                                                                                             gATCCTTTCGCGTCAGGTTTGTGTGATTCGCAAACAGGCGCTTATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCGGCGGGTAACGTTGTGGTACACGGTATTTTTGCCATCGTACCGTTCTTCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGTGTGCTGCAATCAGACCGAAGAGTGCAATACGCGACTGTAGCGAATaaaaaa                                                                                                        <  2:89175/251‑1 (MQ=255)
                                                                                                                                             gATCCTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGCTCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTAAACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATaaaaa                                                                                                         >  2:623371/1‑250 (MQ=255)
                                                                                                                                                 cTTTCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCAATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGACGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcc                                                                                                    >  1:80811/1‑251 (MQ=255)
                                                                                                                                                  tttCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATccc                                                                                                   >  2:496398/1‑251 (MQ=255)
                                                                                                                                                    tCGCGTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATcccc                                                                                                  <  1:542528/250‑1 (MQ=255)
                                                                                                                                                       cgTCAGGTTGGTGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCAGGTCAGCCGAAGTCTATTAAAGAGACGATGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCg                                                                                               >  1:359587/1‑250 (MQ=255)
                                                                                                                                                                   gTGATTCGCAATCAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTAGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGAtc                                                                                   <  2:539025/250‑1 (MQ=255)
                                                                                                                                                                      aTTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGGTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCaa                                                                               <  2:328550/251‑1 (MQ=255)
                                                                                                                                                                          gCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACa                                                                           <  2:175993/251‑1 (MQ=255)
                                                                                                                                                                           cAAACAGGCGCTGATCCTTAACTTTCCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACaa                                                                          <  2:312972/251‑1 (MQ=255)
                                                                                                                                                                              aCAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGATGCATTCAGACCGAAGAGTGCAAGTCGCGACGTTAGAGAATAAAAAAATCCCCCCGAGCGGGGGTATCTCAAATCAATTa                                                                       <  2:38293/251‑1 (MQ=255)
                                                                                                                                                                               cAGGTGCTGATCCGTAACTGACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTTGAAGGTGTGAAGGACGCTGCGGGTAACGGTGTGGTACACGGTATTTTTGCCAGCTGACCGGACTTCATTCAGTTGCTGGAAGGGCCATACGTGGAAACGGCACCAGAAGTTGTTTCTGCAGTCAGAGCTAAGAGTGCAGGACGTGACGGTAGCGAATAAAAAGATCGTCCCGAGCGTGGTGATCTCCGAACAATTAg                                                                      <  2:592982/251‑1 (MQ=255)
                                                                                                                                                                                 ggCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTg                                                                    <  1:354974/251‑1 (MQ=255)
                                                                                                                                                                                    gcTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGa                                                                 <  1:433844/251‑1 (MQ=255)
                                                                                                                                                                                             tAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGGGCAAGACGCGACGTGAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTATTGGGATTCACCa                                                          <  2:384890/249‑1 (MQ=255)
                                                                                                                                                                                              aaCTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAAt                                                        <  1:151950/250‑1 (MQ=255)
                                                                                                                                                                                              aaCTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGACTCACCAAt                                                        <  2:392146/250‑1 (MQ=255)
                                                                                                                                                                                                     ccGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGaa                                                <  1:275935/251‑1 (MQ=255)
                                                                                                                                                                                                       ggTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCAAGCGTACCGTACTGCATGCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCACTCAGACCGAAGAGTGCAAGGCGCGACGTTAGCGAATAAAAAAAACCCCCCGAGCGGGGGGATCGCAAACCAATGAGTGGGATTCACCAAGCGTCAGatcg                                              <  2:56330/251‑4 (MQ=255)
                                                                                                                                                                                                        gTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACgg                                             >  1:523055/1‑251 (MQ=255)
                                                                                                                                                                                                               gAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGAcc                                      >  1:461907/1‑251 (MQ=255)
                                                                                                                                                                                                                        tAAAGAGTCGCTGGAAGGGGTGAAGGACGCTGCGGGTAACGTTGTGGGACACGGTTTTTTTGCCAGCGTACCGTTCTGCATTCAGTTGCTGGAAGGGCCAGACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGAACGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATGAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAGACTGCt                              <  1:496398/250‑1 (MQ=255)
                                                                                                                                                                                                                        tAAAGAGACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGGATTTTTGCCAGCGAACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGCAGTGGTTGCTGCATTCAGACCGAAGAGGGCAAGACGCGACGTTAGCGAATAAAAAAATCACCCCGAGCGGGGGGATCTCACAACAATTAGTGGGATTCACCAATCGGAAGAACGGTGCGACCAAACTGCTc                             >  1:45032/1‑251 (MQ=255)
                                                                                                                                                                                                                           agagACGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATTCGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATGAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGtt                          <  2:281360/251‑1 (MQ=255)
                                                                                                                                                                                                                               aCGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTGTTGCCAGCGTACCGTACGGCATTCAGTTGCGGGAAGGGCCATACGGTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGTGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAgg                      <  2:30822/251‑2 (MQ=255)
                                                                                                                                                                                                                                cGCTGGAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTa                     <  1:433443/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggAAGGTGTGAAGGACGCTGCGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCAGTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTc                 <  2:26212/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gAAGGACGCTGCGGGTAACGTTTTGGTACACTGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGTCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTATTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCGCGTTCAGTACGTCTCCCATCg         <  2:251399/250‑1 (MQ=255)
                                                                                                                                                                                                                                              aaGGACGCTGCGGGTAACGTTGTGGTACACGGTCTTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGcc       >  2:49701/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaGGACGCTGCGGGTAACGATGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGATGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGAGGGGGGATCTCACA‑CAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGCCCAAACTGCTCGTTCAGTACTTCCCCCATCGcc       >  1:245168/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   cGCTGTGGGGAACATTGTGGTTCGCGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCAGAAGTGGTTGCTGCATTCAGACCGAAGAGTTCAAGACGCGACTGTAGCGAATGAAAAAATCCCCCCGAGCGGGGGGATCTCAAAATAATTAGGGGGATTCACTAATCGGCAGTACGGTGCGACCAAACTGCGCGTTCAGTACCTCACCCATCGCCagata  <  2:315155/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
CATCTGGTGCTCACCACGGGCGGAACTGGCCCGGCGCGTCGTGACGTAACGCCCGATGCGACGCTGGCAGTAGCGGACCGCGAGATGCCTGGCTT‑TGGTGAACAGATGCGCCAGATCAGCCTGCATTTTGTACCAACTGCGATCCTTTCGCGTCAGGTGGGCGTGATTCGCAAACAGGCGCTGATCCTTAACTTACCCGGTCAGCCGAAGTCTATTAAAGAGACGCTGGAAGGTGTGAAGGACGCTGAGGGTAACGTTGTGGTACACGGTATTTTTGCCAGCGTACCGTACTGCATTCAGTTGCTGGAAGGGCCATACGTTGAAACGGCACCGGAAGTGGTTGCAGCATTCAGACCGAAGAGTGCAAGACGCGACGTTAGCGAATAAAAAAATCCCCCCGAGCGGGGGGATCTCAAAACAATTAGTGGGATTCACCAATCGGCAGAACGGTGCGACCAAACTGCTCGTTCAGTACTTCACCCATCGCCAGATA  >  NC_000913/9507‑9999

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