Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,790,310 C→T A110V (GCG→GTG)  lhgD → L‑2‑hydroxyglutarate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,790,3100CT98.7% 256.4 / NA 75A110V (GCG→GTG) lhgDL‑2‑hydroxyglutarate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (31/43);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCGTTGCTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGCCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGATACGTACCCGTCAGGGCGGCG  >  NC_000913/2790069‑2790549
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
gctttGCTGGAAAAAGGGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACgg                                                                                                                                                                                                                                        <  2:126040/248‑1 (MQ=255)
 ggttgcgtgAAAAAGAGTCCGCCCCGGCCTGTCCCCAGGCGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGt                                                                                                                                                                                                                                       <  2:489344/243‑1 (MQ=255)
  gTTGCTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGt                                                                                                                                                                                                                                       >  1:314043/1‑250 (MQ=255)
         gAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAgcg                                                                                                                                                                                                                               >  1:342698/1‑251 (MQ=255)
            caaGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGcgcg                                                                                                                                                                                                                             >  2:204266/2‑250 (MQ=255)
              agcgTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGGATCGAGCGCGAAt                                                                                                                                                                                                                          >  2:321146/4‑251 (MQ=255)
              agagTCCCCCCCGGCCTGTCCCCAGACGGGCGACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAAt                                                                                                                                                                                                                          <  2:392159/251‑1 (MQ=255)
                   ccgcccgtgCCTGTCACCGGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACCCGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCtt                                                                                                                                                                                                                     <  2:236107/243‑1 (MQ=255)
                    ccccccGGCGTGTCATCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTGTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGAGGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTa                                                                                                                                                                                                                    <  2:82230/249‑1 (MQ=255)
                     gCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGGCTATTACACGTCCGACAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTACGATCTCGAAATGGACCGGATGCGCGCCTTCTGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGCCTTa                                                                                                                                                                                                                    >  2:516997/1‑250 (MQ=255)
                            ccTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCcgac                                                                                                                                                                                                            >  2:309929/1‑251 (MQ=255)
                            ccTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCcgac                                                                                                                                                                                                            <  1:254099/251‑1 (MQ=255)
                              tGTCTCCAGACGGGCTACACCAGCGTCGTGATCCAGGCCGGGGTCTATTACACGTCCTGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTGTGCGATCAAAACGGCATTCGCTACGACAACTGCGGTAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGAGGCGCGCCTTATGGGAACGCACAGTGGCGAACGGCATCGAGCGCGAATGGCTTAACGCcgac                                                                                                                                                                                                            <  2:221003/249‑1 (MQ=255)
                                tCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAAc                                                                                                                                                                                                        >  1:75627/1‑251 (MQ=255)
                                         gggCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGACCGCACCGTGGCGAACGGTATCGAGCGCGAAAGGCTTAACGCCGACGAACTgcgcgagc                                                                                                                                                                                               >  2:81464/1‑251 (MQ=255)
                                             cacaACAGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTgcgcgagcgcg                                                                                                                                                                                            >  1:145400/1‑250 (MQ=255)
                                              acaacaGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGaa                                                                                                                                                                                          <  1:384256/251‑1 (MQ=255)
                                               caacaGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc                                                                                                                                                                                         <  2:482599/251‑1 (MQ=255)
                                               caacaGCGGCGTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAAc                                                                                                                                                                                         <  2:132255/251‑1 (MQ=255)
                                                      ggcgTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAatat                                                                                                                                                                                  >  1:37602/1‑251 (MQ=255)
                                                        cgTGATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGCACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATc                                                                                                                                                                                 >  2:98694/1‑250 (MQ=255)
                                                           gATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCAcc                                                                                                                                                                              <  2:2817/250‑1 (MQ=255)
                                                           gATCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCAcc                                                                                                                                                                              <  1:13417/250‑1 (MQ=255)
                                                            aTCCATGCCGGGGTCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCgg                                                                                                                                                                            <  1:198445/251‑1 (MQ=255)
                                                                         tCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCAtt                                                                                                                                                               >  1:284796/1‑251 (MQ=255)
                                                                         tCTATTACACGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCAtt                                                                                                                                                               >  1:178249/1‑251 (MQ=255)
                                                                                cacGTCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACAGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGc                                                                                                                                                        >  1:523940/1‑251 (MQ=255)
                                                                                 acgtccGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGcc                                                                                                                                                       >  2:344994/5‑251 (MQ=255)
                                                                                    cccGGCAGCCTGAAGGCGCAGGTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTc                                                                                                                                                    <  2:300817/251‑1 (MQ=255)
                                                                                        gCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCCCCGGGCTCGGCGGCGTTTTTGTGCCGTCCAGc                                                                                                                                                <  1:107093/251‑1 (MQ=255)
                                                                                        gCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGc                                                                                                                                                >  1:306922/1‑251 (MQ=255)
                                                                                          aGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCCCCACAGCTTTTTGCGATCAAAACGGCAATCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCATTATGGAAACTCACAGCGGCGAACGGCATCGAGCGCGAATGGCTTAACGCCGCCGAACTGCGGCAGCGCGACCCGCATATCACCGGGCTCGGTGGCATTCTTGAGCAGTcccgccg                                                                                                                                              >  2:191169/1‑246 (MQ=255)
                                                                                                aaGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGt                                                                                                                                        >  1:66306/1‑251 (MQ=255)
                                                                                                               cTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACgg                                                                                                                         >  2:228582/1‑250 (MQ=255)
                                                                                                                     ggAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGAcggc                                                                                                                   <  1:508223/251‑1 (MQ=255)
                                                                                                                     ggAAACCGCGCCACCAAAGCCTTCTGCGAGCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTAGCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGTTACCGCGACGTGAcggc                                                                                                                   <  2:350908/251‑1 (MQ=255)
                                                                                                                       aaatcgcgCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGAcggcgg                                                                                                                 >  2:222154/5‑251 (MQ=255)
                                                                                                                        aaCCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTGTTGTGCCGTCCAGCGGCATTGTCAGCTGCCGCGACGTGACggcggc                                                                                                                <  2:517887/251‑1 (MQ=255)
                                                                                                                            gcgcCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCCTTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGagg                                                                                                            >  2:285096/1‑249 (MQ=255)
                                                                                                                             cgcCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATgg                                                                                                           <  1:105461/251‑1 (MQ=255)
                                                                                                                             cgcCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACGGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATgg                                                                                                           <  1:370622/251‑1 (MQ=255)
                                                                                                                                caccaAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCa                                                                                                         >  1:182570/1‑250 (MQ=255)
                                                                                                                                    aaaGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAAt                                                                                                    <  2:145400/251‑1 (MQ=255)
                                                                                                                                        ccTTTTGCGATCAAAACGGCCTTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTtc                                                                                                <  2:486008/251‑1 (MQ=255)
                                                                                                                                        ccTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGCTGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGCTGGCAAAAATCTtc                                                                                                >  1:43770/1‑251 (MQ=255)
                                                                                                                                         attttGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGCTGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAACTCTtcc                                                                                               >  2:72217/2‑251 (MQ=255)
                                                                                                                                           tttGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTtcc                                                                                               >  2:508019/1‑249 (MQ=255)
                                                                                                                                            ttGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTtccagt                                                                                            >  2:10757/1‑251 (MQ=255)
                                                                                                                                            ttGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTtccag                                                                                             <  1:321146/250‑1 (MQ=255)
                                                                                                                                             tGCGATCAAAACGGCATTCGCTAAGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTtccagt                                                                                            <  2:229958/250‑1 (MQ=255)
                                                                                                                                              gCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTtccagtcc                                                                                          <  2:411468/251‑1 (MQ=255)
                                                                                                                                                   cAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCagag                                                                                      <  2:192440/250‑1 (MQ=255)
                                                                                                                                                        cGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGCCGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCgggg                                                                                >  1:142756/1‑249 (MQ=255)
                                                                                                                                                           cATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTCTGGGGACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCgggtgt                                                                              >  2:375405/1‑246 (MQ=255)
                                                                                                                                                            aTTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCggggcaa                                                                             >  2:321428/1‑245 (MQ=255)
                                                                                                                                                            aTTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGGAAt                                                                            >  1:455481/1‑251 (MQ=255)
                                                                                                                                                                     gACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATaa                                                                   <  2:326054/251‑1 (MQ=255)
                                                                                                                                                                            gCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGa                                                             <  1:309929/250‑1 (MQ=255)
                                                                                                                                                                            gCGGCAAGATGCTGGTTGCCACGGCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGaa                                                            <  2:284796/251‑1 (MQ=255)
                                                                                                                                                                            gCGGCAAGATGCGGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGTGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGa                                                             <  2:139811/250‑1 (MQ=255)
                                                                                                                                                                             cGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCgggggaaaaatctttaaaaccca                                                             >  2:33262/1‑228 (MQ=255)
                                                                                                                                                                                 aaGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGGGAAAATATCtattaagacggagac                                                         >  2:95972/1‑237 (MQ=255)
                                                                                                                                                                                  aGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGAGGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGc                                                      <  2:280042/251‑1 (MQ=255)
                                                                                                                                                                                    aTGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCgg                                                    <  1:370910/251‑1 (MQ=255)
                                                                                                                                                                                      gCTGGTTGCCACGTCCGAGCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTGAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGGCAGCGGGc                                                  <  2:314043/251‑1 (MQ=255)
                                                                                                                                                                                                  gTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATa                                      <  1:222154/251‑1 (MQ=255)
                                                                                                                                                                                                           cGAAATGGAACGGATGCGCGCCTTATGGGAACGCGCATTGGCGAACGGTATCTAGCGGGAATGGCATAACGCCGACGAACTGCGCGAGCGCTAACCCAATATCACCCGGCTCGCCGGCAGTTTTGTTCCGTCGAGCGGCATTGGCAGCTACCGCGACGTGACGGCGGCGATGGCACAAATCTACCAGTCCACATGCGGTGAAATTACCTATAACGCCTAAGTCAGCGGGCTTAATGAGCTTaaaaag                                 <  2:220108/247‑2 (MQ=255)
                                                                                                                                                                                                           cGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCTAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTGCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATaaaaa                                  <  1:285096/246‑1 (MQ=255)
                                                                                                                                                                                                           cGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCAt                                       <  1:367991/241‑1 (MQ=255)
                                                                                                                                                                                                               aTGGAACGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCgtggt                         <  2:509114/251‑1 (MQ=255)
                                                                                                                                                                                                                    aCGGATGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGATAc                    <  1:392334/251‑1 (MQ=255)
                                                                                                                                                                                                                        aTGCGCGCCTTATGGGAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGTATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGATACGTAc                <  1:452345/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGATACGTACCCGTCAGGgcggcg  <  1:516997/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCAGTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGATACGTACCCGTCAGGgcggcg  <  2:390386/251‑1 (MQ=255)
                                                                                                                                                                                                                                      ggAACGCACAGTGGCGAACGGTATCGAGCGCGAATGGCGTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCAGTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGATACGTACCCGTCAGGgcggcg  <  2:499153/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GCGTTGCTGGAAAAAGAGTCCGCCCCGGCCTGTCACCAGACGGGCCACAACAGCGGCGTGATCCATGCCGGGGTCTATTACACGCCCGGCAGCCTGAAGGCGCAGTTTTGCCTGGCGGGAAACCGCGCCACCAAAGCCTTTTGCGATCAAAACGGCATTCGCTACGACAACTGCGGCAAGATGCTGGTTGCCACGTCCGATCTCGAAATGGAACGGATGCGCGCCTTATGGGAACGCACAGCGGCGAACGGTATCGAGCGCGAATGGCTTAACGCCGACGAACTGCGCGAGCGCGAACCGAATATCACCGGGCTCGGCGGCATTTTTGTGCCGTCCAGCGGCATTGTCAGCTACCGCGACGTGACGGCGGCGATGGCAAAAATCTTCCAGTCCAGAGGCGGTGAAATTATCTATAACGCCGAAGTCAGCGGGCTTAATGAGCATAAAAACGGCGTGGTGATACGTACCCGTCAGGGCGGCG  >  NC_000913/2790069‑2790549

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