Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 746,411 T→C T231A (ACC→GCC)  unknown ← Molybdopterin molybdenumtransferase (EC 2.10.1.1)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1746,4110TC100.0% 133.8 / NA 51T231A (ACC→GCC) unknownMolybdopterin molybdenumtransferase (EC 2.10.1.1)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/25);  total (26/25)

CACTAATTGATAAAATGTCACCAACGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGTGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCACGTTTTCAATGCCAAGGGAGGCGAGTAACGGCAGGGAAATCGCATTTAAT  >  1/746181‑746652
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
cACTAATTGATAAAATGTCACCAACGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGt                                                                                                                                                                                                                               <  1:205137/251‑1 (MQ=255)
    aaTTGATAAAATGTCACCAACGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTAt                                                                                                                                                                                                                            >  1:419600/1‑250 (MQ=255)
          tAAAATGTCACCAACGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCaa                                                                                                                                                                                                                     <  1:526025/251‑1 (MQ=255)
              aTGTCACCAACGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATa                                                                                                                                                                                                                 <  1:16860/251‑1 (MQ=255)
                        cGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCa                                                                                                                                                                                                       >  1:460151/1‑251 (MQ=255)
                        cGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCa                                                                                                                                                                                                       >  1:600823/1‑251 (MQ=255)
                          ccGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAac                                                                                                                                                                                                     <  1:882002/251‑1 (MQ=255)
                             acacCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCaa                                                                                                                                                                                                                     <  1:640469/232‑1 (MQ=255)
                             acacCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAacac                                                                                                                                                                                                   <  1:635282/250‑1 (MQ=255)
                             acacCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAaca                                                                                                                                                                                                    >  1:503837/1‑249 (MQ=255)
                               acCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTc                                                                                                                                                                                                <  1:554822/251‑1 (MQ=255)
                                 cGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGc                                                                                                                                                                                              <  1:217216/251‑1 (MQ=255)
                                          cGGCAAACCGCAGAGCCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCGCTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAAtt                                                                                                                                                                                     <  1:790677/251‑1 (MQ=255)
                                          cGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAAtt                                                                                                                                                                                     >  1:208315/1‑251 (MQ=255)
                                          cGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAAtt                                                                                                                                                                                     <  1:661757/251‑1 (MQ=255)
                                                 ccGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATAc                                                                                                                                                                                                                <  1:786761/217‑1 (MQ=255)
                                                 ccGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCaa                                                                                                                                                                              >  1:240027/1‑251 (MQ=255)
                                                                                       aGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAAttt                                                                                                                                        >  1:222423/1‑251 (MQ=255)
                                                                                        ggTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTg                                                                                                                                       >  1:805153/1‑251 (MQ=255)
                                                                                                  gTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAAtt                                                                                                                                         >  1:183223/1‑239 (MQ=255)
                                                                                                     tAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGc                                                                                                                          <  1:849157/251‑1 (MQ=255)
                                                                                                       aTGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCgg                                                                                                                        <  1:267210/251‑1 (MQ=255)
                                                                                                       aTGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCgg                                                                                                                        <  1:756318/251‑1 (MQ=255)
                                                                                                                       aaTTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATTCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCAc                                                                                                         >  1:711755/1‑250 (MQ=255)
                                                                                                                            aTTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCGCGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATtc                                                                                                   <  1:549184/251‑1 (MQ=255)
                                                                                                                               ttGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATtcatc                                                                                                <  1:865111/251‑1 (MQ=255)
                                                                                                                                tGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATtcatca                                                                                               >  1:140229/1‑251 (MQ=255)
                                                                                                                                 gCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAac                                                                                                                                                                                                     >  1:56738/1‑148 (MQ=255)
                                                                                                                                  cccACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCAc                                                                                              <  1:461742/250‑1 (MQ=255)
                                                                                                                                    cACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCgg                                                                                           >  1:708911/1‑251 (MQ=255)
                                                                                                                                         tctcTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGaa                                                                                      <  1:149653/251‑1 (MQ=255)
                                                                                                                                          ctctAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAacac                                                                                                                                                                                                   <  1:246567/141‑1 (MQ=255)
                                                                                                                                          ctctAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGaaa                                                                                     >  1:367079/1‑251 (MQ=255)
                                                                                                                                            ctAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTc                                                                                                                                                <  1:114534/190‑1 (MQ=255)
                                                                                                                                            ctAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAAt                                                                                                                                          >  1:307258/1‑196 (MQ=255)
                                                                                                                                              aaaaTATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCAc                                                                                                                                                                  <  1:472308/170‑1 (MQ=255)
                                                                                                                                               aaaTATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAacac                                                                                                                                                                                                   >  1:748673/1‑136 (MQ=255)
                                                                                                                                               aaaTATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCAc                                                                                 >  1:565370/1‑250 (MQ=255)
                                                                                                                                                         ggTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTa                                                                      >  1:251788/1‑251 (MQ=255)
                                                                                                                                                            aaaacccGCTTCCCCCACGAAACCGCCGCCGCTGGTAATCCCCAAACCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTAAAATTGCGCCTCGTTATCTGGCAAAATCCCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACAAATTCATCCCCGGTAGAAAGCACCGCCACTTTTAAAc                                                                   <  1:718375/246‑1 (MQ=255)
                                                                                                                                                                             cGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCACg                                                  >  1:677110/1‑251 (MQ=255)
                                                                                                                                                                              ggAACCGCCGCCGCTGGTAACCACCAAACCCGCCTGTTGTGGCGCTTAAATAAAGGCGTTTTCAAATGGCGCCTCGTTATCTGGCAAAATACCGAAATCCACCACTTCGCAGTTTAATTTTGCAAATAATAATTTCACGGCGAAGCGGTTGGTGTCATAAATTGGCCCGGTTTGTAGCGGTTTGCCCACCGGCCCCAATTCATCACCGGAAAAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCAc                                                   <  1:144914/249‑1 (MQ=255)
                                                                                                                                                                                             ggTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCACGTTTTCAATGCCAAggg                                  >  1:366043/1‑251 (MQ=255)
                                                                                                                                                                                              gTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCAc                                                                                              <  1:321786/190‑1 (MQ=255)
                                                                                                                                                                                              gTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCACGTTTTCAATGCCAAGGGa                                 >  1:291698/1‑251 (MQ=255)
                                                                                                                                                                                                   cACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCACGTTTTCAATGCCAAGGGAGGCGa                            >  1:144571/1‑251 (MQ=255)
                                                                                                                                                                                                               cTGTTGTTGCGCTTCAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGc                                                                              >  1:613800/1‑189 (MQ=255)
                                                                                                                                                                                                                             cAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGgtg                                                                                                                                                  >  1:211287/1‑107 (MQ=255)
                                                                                                                                                                                                                             cAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCgg                                                                                                                        <  1:681408/133‑1 (MQ=255)
                                                                                                                                                                                                                             cAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGtt                                                                                                                      >  1:561962/1‑135 (MQ=255)
                                                                                                                                                                                                                             cAATAAAGGCGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGGTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCACGTTTTCAATGCCAAGGGAGGCGAGTAACGGCAGGGAAATCGCATTtaat  >  1:38249/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
CACTAATTGATAAAATGTCACCAACGCCGACACCGGGTTGCCCGGCAAACCGCAGAACCAGGCGTTTTTTAATTTGCCGAAAGCAAAAGGTTTACCCGGTTTAATGGCGATTTTCCAGAAATTGATTTTGCCCACTTTCTCTAAAATATGTTTGGTGAAATCCGCTTCGCCCACGGAAACGCCGCCGCTGGTAATCACCAAATCCGCCTGTTGTTGCGCTTCAATAAAGGTGTTTTCAAATTGCGCCTCGTTATCTGGCAAAATACCGAAATCCAACACTTCGCAGTTTAATTTTTCCAATAATAATTTCACGGCGAAACGGTTGGTGTCATAAATTTGCCCGGTTTGTAGCGGTTTGCCCACCGGCACCAATTCATCACCGGTAGAAAGCACCGCCACTTTTAAACGCGGGAAAACCTTCACGTTTTCAATGCCAAGGGAGGCGAGTAACGGCAGGGAAATCGCATTTAAT  >  1/746181‑746652

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