Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,520,256 A→G H329R (CAT→CGT)  unknown → Ribonucleotide reductase of class Ia (aerobic), alpha subunit (EC 1.17.4.1)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,520,2560AG100.0% 226.6 / NA 79H329R (CAT→CGT) unknownRibonucleotide reductase of class Ia (aerobic), alpha subunit (EC 1.17.4.1)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (42/37);  total (42/37)

GACGTGTATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCATCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCGTGGCGCCGACCAAATCCAGCAGCTTTATTT  >  1/1520009‑1520501
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gACGTGTATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGtcgtca                                                                                                                                                                                                                                                    >  1:687854/1‑251 (MQ=255)
 aCGTGTATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTcgtcat                                                                                                                                                                                                                                                   <  1:20976/251‑1 (MQ=255)
  cGTGTATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTcgtcat                                                                                                                                                                                                                                                   >  1:261680/1‑250 (MQ=255)
    tgtATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCg                                                                                                                                                                                                                                                 <  1:786399/250‑1 (MQ=255)
    tgtATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGt                                                                                                                                                                                                                                                >  1:577235/1‑251 (MQ=255)
    tgtATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGt                                                                                                                                                                                                                                                >  1:861206/1‑251 (MQ=255)
    tgtATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGt                                                                                                                                                                                                                                                >  1:733004/1‑251 (MQ=255)
    tgtATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGt                                                                                                                                                                                                                                                >  1:338814/1‑251 (MQ=255)
     gtATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTg                                                                                                                                                                                                                                               >  1:892080/1‑251 (MQ=255)
              tCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTggg                                                                                                                                                                                                                                        >  1:739149/1‑249 (MQ=255)
                  aaCATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  <  1:644388/251‑1 (MQ=255)
                  aaCATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  <  1:535909/251‑1 (MQ=255)
                  aaCATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  <  1:285665/251‑1 (MQ=255)
                    cATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGgtg                                                                                                                                                                                                                                <  1:675045/251‑1 (MQ=255)
                     atgatgATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  <  1:559334/248‑1 (MQ=255)
                      tgatgaTTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGtgtg                                                                                                                                                                                                                              >  1:732453/1‑251 (MQ=255)
                            ttCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTgctgggc                                                                                                                                                                                                                        <  1:19800/251‑1 (MQ=255)
                            ttCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTgctgggc                                                                                                                                                                                                                        <  1:129371/251‑1 (MQ=255)
                            ttCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTgctgggc                                                                                                                                                                                                                        <  1:240689/251‑1 (MQ=255)
                                        aTTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTAt                                                                                                                                                                                                            >  1:360705/1‑251 (MQ=255)
                                          ttATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTAttt                                                                                                                                                                                                          <  1:812477/251‑1 (MQ=255)
                                            aTGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGc                                                                                                                                                                                                        <  1:685464/251‑1 (MQ=255)
                                               tAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAAt                                                                                                                                                                                                     >  1:179765/1‑251 (MQ=255)
                                                   gAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGaa                                                                                                                                                                                                 >  1:577656/1‑251 (MQ=255)
                                                                             cGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAg                                                                                                                                                                       <  1:236052/251‑1 (MQ=255)
                                                                                aaGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACTTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCa                                                                                                                                                                    <  1:435413/251‑1 (MQ=255)
                                                                                  gTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTa                                                                                                                                                                                                             >  1:247406/1‑208 (MQ=255)
                                                                                        tgtgtTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAatat                                                                                                                                                            <  1:667620/251‑1 (MQ=255)
                                                                                        tgtgtTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAatat                                                                                                                                                            <  1:20077/251‑1 (MQ=255)
                                                                                            ttgttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTgctgggc                                                                                                                                                                                                                        >  1:543773/1‑187 (MQ=255)
                                                                                            ttgttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAAc                                                                                                                                                         <  1:280262/250‑1 (MQ=255)
                                                                                            ttgttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAAc                                                                                                                                                         <  1:846135/250‑1 (MQ=255)
                                                                                             tgttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACaa                                                                                                                                                       <  1:536380/251‑1 (MQ=255)
                                                                                             tgttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACaa                                                                                                                                                       <  1:412486/251‑1 (MQ=255)
                                                                                             tgttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACaa                                                                                                                                                       <  1:446882/251‑1 (MQ=255)
                                                                                             tgttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACaa                                                                                                                                                       <  1:440713/251‑1 (MQ=255)
                                                                                              gttTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  >  1:175047/1‑175 (MQ=255)
                                                                                                       ctATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACaa                                                                                                                                                       >  1:19176/1‑241 (MQ=255)
                                                                                                       ctATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGCGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGAtttt                                                                                                                                              >  1:416923/1‑250 (MQ=255)
                                                                                                           gAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGc                                                                                                                                                                                                        >  1:489301/1‑188 (MQ=255)
                                                                                                            aaTTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGgtg                                                                                                                                                                                                                                <  1:562819/163‑1 (MQ=255)
                                                                                                                 gCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAAc                                                                                                                                                         <  1:398247/229‑1 (MQ=255)
                                                                                                                  ctctATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  >  1:658846/1‑155 (MQ=255)
                                                                                                                    ctATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAgg                                                                                                                                >  1:798674/1‑251 (MQ=255)
                                                                                                                      aTTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  <  1:484699/151‑1 (MQ=255)
                                                                                                                      aTTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGc                                                                                                                                                                      >  1:338165/1‑211 (MQ=255)
                                                                                                                      aTTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAAc                                                                                                                                                         >  1:363054/1‑224 (MQ=255)
                                                                                                                                aTGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAAc                                                                                                                                                         <  1:327358/214‑1 (MQ=255)
                                                                                                                                  ggAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTaaaa                                                                                                                  >  1:209640/1‑251 (MQ=255)
                                                                                                                                    aaaGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAg                                                                                                                 <  1:772705/250‑1 (MQ=255)
                                                                                                                                           ccGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGAt                                                                                                                                                 >  1:746045/1‑211 (MQ=255)
                                                                                                                                                 cACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGc                                                                                                                                                                                                        <  1:892218/150‑1 (MQ=255)
                                                                                                                                                 cACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAagag                                                                                                      >  1:466429/1‑248 (MQ=255)
                                                                                                                                                                     ttCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAAc                                                                                >  1:45995/1‑250 (MQ=255)
                                                                                                                                                                       ccTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCgg                                                                             >  1:787152/1‑251 (MQ=255)
                                                                                                                                                                                   gatgaGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAg                                                                 >  1:473466/1‑251 (MQ=255)
                                                                                                                                                                                              tCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAAc                                                                                                                                                         <  1:265083/152‑1 (MQ=255)
                                                                                                                                                                                              tCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACaa                                                                                                                                                       >  1:750407/1‑154 (MQ=255)
                                                                                                                                                                                              tCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAaga                                                                                                       >  1:843275/1‑202 (MQ=255)
                                                                                                                                                                                              tCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGt                                                         >  1:574185/1‑248 (MQ=255)
                                                                                                                                                                                              tCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCg                                                       >  1:416090/1‑250 (MQ=255)
                                                                                                                                                                                                   gCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACaa                                                                                                                                                       <  1:152889/149‑1 (MQ=255)
                                                                                                                                                                                                   gCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAgg                                                                <  1:100401/236‑1 (MQ=255)
                                                                                                                                                                                                          gcgcCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATaccac                                           >  1:499939/1‑250 (MQ=255)
                                                                                                                                                                                                             ccAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTa                                       >  1:82249/1‑251 (MQ=255)
                                                                                                                                                                                                                atatTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTgctgggc                                                                                                                                                                                                                        >  1:15456/1‑71 (MQ=255)
                                                                                                                                                                                                                atatTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATg                                    >  1:726669/1‑251 (MQ=255)
                                                                                                                                                                                                                atatTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATg                                    <  1:14278/251‑1 (MQ=255)
                                                                                                                                                                                                                  atTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAgg                                                                                                                                                                                                                                  >  1:517363/1‑59 (MQ=255)
                                                                                                                                                                                                                  atTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCt                                                                                                    <  1:402672/185‑1 (MQ=255)
                                                                                                                                                                                                                  atTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATac                                              <  1:626439/239‑1 (MQ=255)
                                                                                                                                                                                                                  atTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATg                                    >  1:415939/1‑249 (MQ=255)
                                                                                                                                                                                                                  atTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAg                                  >  1:521377/1‑251 (MQ=255)
                                                                                                                                                                                                                    tCAATACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCg                                <  1:787362/251‑1 (MQ=255)
                                                                                                                                                                                                                       aTACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCGCAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAatat                                                                                                                                                            <  1:599725/124‑1 (MQ=255)
                                                                                                                                                                                                                       aTACATGGAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCGTgg                             >  1:111504/1‑251 (MQ=255)
                                                                                                                                                                                                                              gAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGgtg                                                                                                                                                                                                                                >  1:605801/1‑49 (MQ=255)
                                                                                                                                                                                                                              gAAAAAGCCAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCGTGGCGCCGa                       >  1:672191/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ccAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCGTGGCGCCGACCAAATc                <  1:301401/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cAATAAATTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCGTGGCGCCGACCAAATCca              >  1:708395/1‑251 (MQ=255)
                                                                                                                                                                                                                                        taaaaTTCGCCCGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTg                                                                                   >  1:184886/3‑180 (MQ=255)
                                                                                                                                                                                                                                                  ccGTCGTCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCGTGGCGCCGACCAAATCCAGCAGCTTTAttt  <  1:335121/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GACGTGTATAAAGATCGCAACATGATGATTCACGCCTCCAATTTATGTAGTGAAATCGCCTTGCCGTCCAGCGAAGACGAAAGTTTCGTGTGTTGTTTATCTTCTATGAATTTGCTCTATTTTGATGAATGGAAAGACACCGACGCACCGGAAGTGCTGACCTATTTCCTGGATACGGTGATGAGCGAGTTCATTGCACAAAGCGCCAATATTCAATACATGGAAAAAGCCAATAAATTCGCCCGTCATCATCGTGCGCTGGGCTTAGGTGTGCTGGGCTGGCACAGCTATTTGCAATCGAAAAATATCGCTTTCGACAGCTTTGAAGCCATGCAATATAACAATATGATTTTCAAATTGTTACAGGAAAAAACCTTAAAAGCCTCGAAAGAGCTGGCCCATCGTTTCGGTGAACCGGAAATTCTGAAAGGTTATGGTCGCCGTAATACCACCTTAATGAGCGTGGCGCCGACCAAATCCAGCAGCTTTATTT  >  1/1520009‑1520501

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