Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,215,670 G→T D65Y (GAT→TAT)  aceB → malate synthase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,215,6700GT100.0% 126.5 / NA 36D65Y (GAT→TAT) aceBmalate synthase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (21/15);  total (21/15)

GTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTT  >  NC_000913/4215440‑4215913
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
gTGGGGAAGTTTTCGGATACATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATCTTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCt                                                                                                                                                                                                                                 >  2:79538/1‑251 (MQ=255)
        gTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTcg                                                                                                                                                                                                                         >  1:134721/1‑251 (MQ=255)
             cGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATg                                                                                                                                                                                                                    <  1:249397/251‑1 (MQ=255)
                                       tGACTGAACAGGCAACAACCACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAACTTCTTACTGCCGAAGCGGTAGAATTTCTGAATGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTATGGCAGAGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCt                                                                                                                                                                                          >  2:215043/1‑251 (MQ=255)
                                         aCTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCTGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGc                                                                                                                                                                                        <  1:525578/251‑1 (MQ=255)
                                              aCAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACt                                                                                                                                                                                   >  1:7533/1‑251 (MQ=255)
                                                      caacaacCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATTAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGAc                                                                                                                                                                           >  1:550116/1‑251 (MQ=255)
                                                                      ggCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGAt                                                                                                                                                             <  1:364716/249‑1 (MQ=255)
                                                                                 ggCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGt                                                                                                                                                 >  2:141201/1‑250 (MQ=255)
                                                                                                   agaAGCAAATTCTTACTGCCGAAGCGTTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACATCGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCCGAATGGAAAATTCGCGGGATTCCTGCGGACTGAGAAGACCGCCGCGTAGAGATAACTGTCCCGGTAGAGCGCAAGATGTTGATc                                                                                                                              <  2:550116/251‑1 (MQ=255)
                                                                                                                      cGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAAtg                                                                                                           >  2:427961/1‑251 (MQ=255)
                                                                                                                       gAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCACtgt                                                                                                          >  2:136105/1‑251 (MQ=255)
                                                                                                                            ggTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAg                                                                                                     <  2:7533/251‑1 (MQ=255)
                                                                                                                              tAGAGTTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTACATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCCCAACGCCAGTGTGAAAGTc                                                                                                   >  2:216375/1‑251 (MQ=255)
                                                                                                                                   tttCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTAt                                                                                              >  1:194006/1‑251 (MQ=255)
                                                                                                                                    ttCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGGCTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTAt                                                                                              <  1:90631/250‑1 (MQ=255)
                                                                                                                                        gactgaGCTGGAGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTCTTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCACCGCCAATGTGCAAGTCTTTATGGCCg                                                                                         >  2:345045/1‑251 (MQ=255)
                                                                                                                                                   tGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGa                                                                                  >  1:388379/1‑247 (MQ=255)
                                                                                                                                                        cATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACt                                                                         <  2:277789/251‑1 (MQ=255)
                                                                                                                                                              aCGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCAc                                                                    <  1:427961/250‑1 (MQ=255)
                                                                                                                                                               cGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCAcc                                                                   >  1:428167/1‑250 (MQ=255)
                                                                                                                                                                  cacgaCGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGAc                                                               >  1:196336/5‑251 (MQ=255)
                                                                                                                                                                                  tCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCGCCAGACTGGAACAAAGTGATc                                                <  2:376085/250‑1 (MQ=255)
                                                                                                                                                                                  tCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGAt                                                 <  2:272976/249‑1 (MQ=255)
                                                                                                                                                                                     ggCAGCGCGCATTCATCAGCAGCAAGATATTGATCACGGAACGTTGCCTTATGTTATTTCGGAAACAGCTGCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATTTGAAAGTCTTTATGGCCGATTTCGAATATTCACTGGCTCCAGACTGGAACAAAGTGATCGACg                                            <  2:14867/251‑1 (MQ=255)
                                                                                                                                                                                         gcgcgcATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCa                                        <  1:99197/251‑1 (MQ=255)
                                                                                                                                                                                          cgcgcATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGATGGGCaa                                       <  2:388379/251‑1 (MQ=255)
                                                                                                                                                                                                   agcagcagcAAGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTg                              >  2:172831/1‑251 (MQ=255)
                                                                                                                                                                                                             aGATATTGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGACGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGAGAAAGTCTTTATGTCCGATTTCGAAGATTCACTGGCACCa                                                                  >  1:181721/1‑205 (MQ=255)
                                                                                                                                                                                                             aGATATTGATAACGGAACGTTGCCTTATTTTAGTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCATCGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCa                                                                  <  2:181721/205‑1 (MQ=255)
                                                                                                                                                                                                                  ttGATAGCGGAACGGTGCCTTATTTTATTGCGGAAACAGCTTCCATTAGCGATGCTGATTGGAAAATTCGCGGAATTCCTGCGGACTTAGAAGACAGCCGCGTAGAGATAACTGGCACGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCCAGGTGGAAGTCTTTATGGCCGAGTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAAc               >  2:34248/1‑251 (MQ=255)
                                                                                                                                                                                                                  ttGATAACGGAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGTGATTCCTGCGGACTTAGAGGACCGCCGCGTAGAGATAACTGGCCCGTTAGCGCGCAAGAGGGTGATCAACGCGCTCAACGCCAATGGGAAAGCCTTTATGGCGGATTTGGAAGTTTCCCTGGCACCAGACGGGACAAAAGAGATCAACGGGCAAATTAACTTACGGGATGCGGTTaa                >  2:411758/1‑250 (MQ=255)
                                                                                                                                                                                                                          ggaccGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCACGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACACCAATGTGACAGTCTTTATGGACGATTTCGAAGATTCACTGGCACCAGACTGGAACAAATTGATCGACGGGCCAATTAACCTGCGAGATGCGGTTAACGGCACCAt       >  2:523770/5‑251 (MQ=255)
                                                                                                                                                                                                                          ggAACGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCAt       >  1:170564/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCGTTGCCTTATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAg    <  1:74496/251‑1 (MQ=255)
                                                                                                                                                                                                                                ttGCCTTATTTTATTTCGGACACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGCCCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCACCGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGtt  >  1:450060/1‑250 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
GTGGGGAAGTTTTCGGATCCATGACGAGGAGCTGCACGATGACTGAACAGGCAACAACAACCGATGAACTGGCTTTCACAAGGCCGTATGGCGAGCAGGAGAAGCAAATTCTTACTGCCGAAGCGGTAGAATTTCTGACTGAGCTGGTGACGCATTTTACGCCACAACGCAATAAACTTCTGGCAGCGCGCATTCAGCAGCAGCAAGATATTGATAACGGAACGTTGCCTGATTTTATTTCGGAAACAGCTTCCATTCGCGATGCTGATTGGAAAATTCGCGGGATTCCTGCGGACTTAGAAGACCGCCGCGTAGAGATAACTGGCCCGGTAGAGCGCAAGATGGTGATCAACGCGCTCAACGCCAATGTGAAAGTCTTTATGGCCGATTTCGAAGATTCACTGGCACCAGACTGGAACAAAGTGATCGACGGGCAAATTAACCTGCGTGATGCGGTTAACGGCACCATCAGTT  >  NC_000913/4215440‑4215913

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