Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 397,951 C→T R230H (CGC→CAC)  unknown ← Queuine tRNA‑ribosyltransferase (EC 2.4.2.29)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1397,9510CT98.8% 209.3 / ‑5.9 81R230H (CGC→CAC) unknownQueuine tRNA‑ribosyltransferase (EC 2.4.2.29)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (45/35);  minor base G (1/0);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTGTAACAGTCACAGTGCGGATCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGCGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCGTACATT  >  1/397711‑398197
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
gtgtGACAGTCACAGTGCGGATCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCatat                                                                                                                                                                                                                                               >  1:707273/1‑250 (MQ=255)
    aaCAGTCACAGTGCGGATCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATAtctt                                                                                                                                                                                                                                            >  1:471075/1‑249 (MQ=255)
              gTGCGGATCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCg                                                                                                                                                                                                                                 >  1:71524/1‑250 (MQ=255)
              gTGCGGATCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCg                                                                                                                                                                                                                                 >  1:365683/1‑250 (MQ=255)
               tGCGGATCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGcc                                                                                                                                                                                                                               <  1:542073/251‑1 (MQ=255)
                     tCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGc                                                                                                                                                                                                                         <  1:581171/251‑1 (MQ=255)
                          cggGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAAc                                                                                                                                                                                                                    <  1:896851/251‑1 (MQ=255)
                          cggGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAAc                                                                                                                                                                                                                    <  1:856590/251‑1 (MQ=255)
                          cggGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAAc                                                                                                                                                                                                                    <  1:839550/251‑1 (MQ=255)
                          cggGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAAc                                                                                                                                                                                                                    <  1:28407/251‑1 (MQ=255)
                            ggCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCg                                                                                                                                                                                                                  >  1:738344/1‑251 (MQ=255)
                                       tcTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATAtctt                                                                                                                                                                                                                                            >  1:387813/1‑214 (MQ=255)
                                         tttGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGc                                                                                                                                                                                                     <  1:297878/251‑1 (MQ=255)
                                         tttGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGc                                                                                                                                                                                                     <  1:858401/251‑1 (MQ=255)
                                         tttGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGc                                                                                                                                                                                                     <  1:871681/251‑1 (MQ=255)
                                                           gATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGccc                                                                                                                                                                                                                              >  1:594545/1‑208 (MQ=255)
                                                            aTTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCGATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGAt                                                                                                                                                                                  <  1:312059/251‑1 (MQ=255)
                                                            aTTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGAt                                                                                                                                                                                  <  1:21808/251‑1 (MQ=255)
                                                                          gTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTg                                                                                                                                                                                    <  1:324979/235‑1 (MQ=255)
                                                                          gTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTg                                                                                                                                                                    >  1:609483/1‑251 (MQ=255)
                                                                                     aTAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTgcgc                                                                                                                                                         <  1:824632/251‑1 (MQ=255)
                                                                                          tGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAAtt                                                                                                                                                     >  1:610244/1‑250 (MQ=255)
                                                                                            cccATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAAttctt                                                                                                                                                  >  1:55347/1‑251 (MQ=255)
                                                                                            cccATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAAttctt                                                                                                                                                  >  1:427528/1‑251 (MQ=255)
                                                                                              cATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCttt                                                                                                                                                                                                                                        >  1:110125/1‑163 (MQ=255)
                                                                                                       cATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCAc                                                                                                                                                                                                             >  1:363388/1‑181 (MQ=255)
                                                                                                       cATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGcc                                                                                                                                                                            <  1:252457/214‑1 (MQ=255)
                                                                                                       cATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCg                                                                                                                                       >  1:407740/1‑251 (MQ=255)
                                                                                                              ggTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGAtt                                                                                                                                                                                 >  1:482449/1‑202 (MQ=255)
                                                                                                                ttGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGAt                                                                                                                               <  1:564218/250‑1 (MQ=255)
                                                                                                                           cAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAg                                                                                                                   >  1:840122/1‑251 (MQ=255)
                                                                                                                             aTCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGatat                                                                                                                                                                                      <  1:768399/182‑1 (MQ=255)
                                                                                                                                 aaCATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGtt                                                                                                             >  1:667916/1‑251 (MQ=255)
                                                                                                                                    atatCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGc                                                                                                                                                                                                     >  1:861861/1‑160 (MQ=255)
                                                                                                                                     tatCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTg                                                                                                                                                                                    >  1:509546/1‑176 (MQ=255)
                                                                                                                                     tatCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAgg                                                                                                                                            <  1:92749/216‑1 (MQ=255)
                                                                                                                                     tatCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTttg                                                                                                         >  1:756184/1‑251 (MQ=255)
                                                                                                                                       tCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTttggt                                                                                                       >  1:770964/1‑249 (MQ=255)
                                                                                                                                          aTGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATAtcttc                                                                                                                                                                                                                                           <  1:335171/116‑1 (MQ=255)
                                                                                                                                             ccGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTaa                                                                                                 <  1:523908/251‑1 (MQ=255)
                                                                                                                                                cgACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTc                                                                                              <  1:741918/251‑1 (MQ=255)
                                                                                                                                                         ccTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACg                                                                                     <  1:645799/251‑1 (MQ=255)
                                                                                                                                                                 tAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGc                                                                             <  1:667753/251‑1 (MQ=255)
                                                                                                                                                                    aTCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGcc                                                                                                                                                                            <  1:349217/153‑1 (MQ=255)
                                                                                                                                                                               ttCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAg                                                               >  1:531158/1‑251 (MQ=255)
                                                                                                                                                                               ttCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAg                                                               >  1:735908/1‑251 (MQ=255)
                                                                                                                                                                                tCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAg                                                               >  1:66596/1‑250 (MQ=255)
                                                                                                                                                                                       cccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTc                                                                                              >  1:880244/1‑212 (MQ=255)
                                                                                                                                                                                       cccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGa                                                       >  1:654107/1‑251 (MQ=255)
                                                                                                                                                                                       cccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGa                                                       >  1:37814/1‑251 (MQ=255)
                                                                                                                                                                                       cccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGa                                                       >  1:731337/1‑251 (MQ=255)
                                                                                                                                                                                        ccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTTGCCCACGGCTAAACCGCCCAc                                                                                                                                                                                                             >  1:742621/1‑100 (MQ=255)
                                                                                                                                                                                        ccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGatat                                                                                                                                                                                      >  1:479410/1‑123 (MQ=255)
                                                                                                                                                                                        ccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGAt                                                                                                                               <  1:165838/178‑1 (MQ=255)
                                                                                                                                                                                        ccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGa                                                       >  1:806638/1‑250 (MQ=255)
                                                                                                                                                                                        ccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGAc                                                      <  1:96633/251‑1 (MQ=255)
                                                                                                                                                                                        ccATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGGGGTGCATATCTTCTTTTGGGTTGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGag                                                      >  1:151155/1‑250 (MQ=255)
                                                                                                                                                                                             aaaTAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACAttt                                                  >  1:658782/1‑250 (MQ=255)
                                                                                                                                                                                               aTAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAg                                                                                                                                                                      >  1:264228/1‑132 (MQ=255)
                                                                                                                                                                                               aTAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGGTAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCa                                               >  1:445154/1‑251 (MQ=255)
                                                                                                                                                                                                   cgcgGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATgg                                            <  1:590782/250‑1 (MQ=255)
                                                                                                                                                                                                     cgGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTc                                                 <  1:40540/243‑1 (MQ=255)
                                                                                                                                                                                                         ttGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGAt                                          >  1:178061/1‑246 (MQ=255)
                                                                                                                                                                                                           gTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTgg                                   >  1:825037/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaTCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTa                           <  1:107425/250‑1 (MQ=255)
                                                                                                                                                                                                                    aaTCTGCGGGCAGATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTa                           <  1:329468/250‑1 (MQ=255)
                                                                                                                                                                                                                                gATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGcc                                                                                                                                                                            >  1:846741/1‑93 (MQ=255)
                                                                                                                                                                                                                                gATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCTAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGt              >  1:738641/1‑251 (MQ=255)
                                                                                                                                                                                                                                gATGTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAgc                                                              <  1:502453/203‑1 (MQ=255)
                                                                                                                                                                                                                                   gTACTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGGAgg           >  1:115737/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAAc                                                                                                                     >  1:123521/1‑142 (MQ=255)
                                                                                                                                                                                                                                      cTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCAt                                >  1:730110/1‑227 (MQ=255)
                                                                                                                                                                                                                                      cTCCAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCg        >  1:587187/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAAc                                                                                                                                                                >  1:485037/1‑97 (MQ=255)
                                                                                                                                                                                                                                        ccAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCaa                                                                                                                                               <  1:116705/114‑1 (MQ=255)
                                                                                                                                                                                                                                        ccAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGAt                                          >  1:807346/1‑215 (MQ=255)
                                                                                                                                                                                                                                        ccAAAATGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCGTa      >  1:36835/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aaaTGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCGCCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACGTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGGCCAGCGCAGGGACGTTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCGTACAt   <  1:99388/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaTGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCGTACAt   <  1:882535/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaTGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCGTACAt   <  1:99784/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aaTGTGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCGTACAtt  <  1:891350/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
GTGTAACAGTCACAGTGCGGATCCGGCGGGCTGGTGTCGTCTTTGTATTTGGCATTACGGATTTTAACGATGCCGTCGGTGACGAATAAATGCCCATTACGCGCATTGCGGGTTGGCATCACGCAATCGAACATATCAATGCCGCGACGTACGCCTTCCACTAAATCTTCCGGTTTTCCTACACCCATTAAATAGCGCGGTTTGTCGGCGGGAATCTGCGGGCAGATGTACTCCAAAATGCGGTGCATATCTTCTTTTGGTTCGCCCACGGCTAAACCGCCCACCGCATAGCCGTCAAAGCCGATATTGATTAAGCCTTCCAGTGAAACTTTGCGCAATTCTTCAAAGGTTCCGCCTTGGATAATGCCGAACAGCGCGTTTTTGTTGCCTAATTCGTCAAAACGAGCGCGGCTGCGTTGCGCCCAGCGCAGGGACATTTCCATGGATTGCTTGGCATAATTAAACGTGGCGGGGTAGGGCGTACATT  >  1/397711‑398197

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