Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,122,056 A→G A184A (GCT→GCC pyrC ← dihydroorotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,122,0560AG100.0% 216.3 / NA 60A184A (GCT→GCCpyrCdihydroorotase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/29);  total (31/29)

ATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGC  >  NC_000913/1121921‑1122191
                                                                                                                                       |                                                                                                                                       
atTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCAt                                                                                                                                     <  2:231351/140‑1 (MQ=255)
   aCGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCtt                                                                                                                                  <  1:832483/140‑1 (MQ=255)
         gAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTtggtgg                                                                                                                            >  1:468365/1‑140 (MQ=255)
                ggTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGAtgtg                                                                                                                     >  2:306165/1‑140 (MQ=255)
                ggTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGAtgtg                                                                                                                     <  1:673193/140‑1 (MQ=255)
                 gTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGc                                                                                                                    >  2:87509/1‑140 (MQ=255)
                         tACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAAc                                                                                                            >  1:1008594/1‑140 (MQ=255)
                              gTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACtt                                                                                                       >  2:748864/1‑140 (MQ=255)
                              gTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACtt                                                                                                       >  2:494413/1‑140 (MQ=255)
                              gTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACtt                                                                                                       >  1:1089309/1‑140 (MQ=255)
                              gTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACACAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACtt                                                                                                       >  2:329668/1‑140 (MQ=255)
                               tGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACttt                                                                                                      >  2:710703/1‑140 (MQ=255)
                               tGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACttt                                                                                                      >  2:76173/1‑140 (MQ=255)
                                 cGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCa                                                                                                    >  2:764339/1‑140 (MQ=255)
                                 cGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCa                                                                                                    >  1:819138/1‑140 (MQ=255)
                                  ggACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAg                                                                                                   <  1:710703/140‑1 (MQ=255)
                                        acgcCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGt                                                                                             >  2:247541/1‑140 (MQ=255)
                                                 aCCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTg                                                                                    >  2:797429/1‑140 (MQ=255)
                                                      cATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCa                                                                               >  1:133475/1‑140 (MQ=255)
                                                      cATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCa                                                                               >  2:928889/1‑140 (MQ=255)
                                                      cATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATAGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCa                                                                               >  2:130064/1‑140 (MQ=255)
                                                        tatGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCaga                                                                             >  1:1093724/1‑140 (MQ=255)
                                                          tGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAgg                                                                           <  1:1026242/140‑1 (MQ=255)
                                                            gTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGtt                                                                         <  1:748864/140‑1 (MQ=255)
                                                                cGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCAt                                                                     >  2:941621/1‑140 (MQ=255)
                                                                cGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCAt                                                                     >  1:295298/1‑140 (MQ=255)
                                                                 ggTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATc                                                                    >  2:658586/1‑140 (MQ=255)
                                                                  gTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCa                                                                   <  1:87509/140‑1 (MQ=255)
                                                                  gTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCa                                                                   <  2:819138/140‑1 (MQ=255)
                                                                   ttAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCAc                                                                  <  2:295298/140‑1 (MQ=255)
                                                                     aaaCATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGc                                                                <  1:797429/140‑1 (MQ=255)
                                                                     aaaCATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGc                                                                <  2:1093724/140‑1 (MQ=255)
                                                                     aaaCATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGc                                                                <  1:910286/140‑1 (MQ=255)
                                                                             gATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCtata                                                        >  1:954317/1‑140 (MQ=255)
                                                                             gATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCtata                                                        >  1:811940/1‑140 (MQ=255)
                                                                             gATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCtata                                                        >  1:925134/1‑140 (MQ=255)
                                                                             gATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCtata                                                        >  1:82994/1‑140 (MQ=255)
                                                                               tgctgcGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATaaa                                                      <  1:76173/140‑1 (MQ=255)
                                                                                  tgcGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAgcg                                                   <  1:361887/140‑1 (MQ=255)
                                                                                   gcGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAgcgc                                                  >  2:1091608/1‑140 (MQ=255)
                                                                                        gTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTc                                             <  2:1008594/140‑1 (MQ=255)
                                                                                         tGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCa                                            <  1:764339/140‑1 (MQ=255)
                                                                                              gtggCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAg                                                     >  2:602898/1‑126 (MQ=255)
                                                                                              gtggCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAg                                                     <  1:602898/126‑1 (MQ=255)
                                                                                                ggCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCaa                                     <  2:954317/140‑1 (MQ=255)
                                                                                                  cagccagcCGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGCTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCaaaa                                   <  2:1035629/140‑1 (MQ=255)
                                                                                                     ccagccGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATg                                >  2:619594/1‑140 (MQ=255)
                                                                                                         ccGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCg                             <  2:966256/139‑1 (MQ=255)
                                                                                                         ccGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCg                             >  1:966256/1‑139 (MQ=255)
                                                                                                          cGTTCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGat                           <  1:928889/140‑1 (MQ=255)
                                                                                                            ttCATTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGatat                         <  2:925134/140‑1 (MQ=255)
                                                                                                               aTTTCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTg                      <  2:624025/140‑1 (MQ=255)
                                                                                                                 ttCCGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTTACGATCAAAAATGTCGATATCTGCa                    <  2:1089309/140‑1 (MQ=255)
                                                                                                                   ccGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCAtg                  <  1:306165/140‑1 (MQ=255)
                                                                                                                    cGTCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCAtgt                 >  2:726598/1‑140 (MQ=255)
                                                                                                                      tCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGatat                         <  2:647628/130‑1 (MQ=255)
                                                                                                                      tCACGGACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGatat                         >  1:647628/1‑130 (MQ=255)
                                                                                                                          ggACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACt           <  2:133475/140‑1 (MQ=255)
                                                                                                                           gACATAGTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACtt          <  1:1082731/140‑1 (MQ=255)
                                                                                                                                 gTCGGCGGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCAt    <  1:329668/140‑1 (MQ=255)
                                                                                                                                   cggcggCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGc  <  1:1112336/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
ATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGCTGCGGAGTGATGGTGGCAGCCAGCCGTTCATTTCCGTCACGGACATAGTCGGCAGCATCTTTGGTGGTGATGTGCTCAAAAACGACTTTCAGCGCAGTCAGGCGCTGGCGCAGAGGTTCCATCACGCTTTCTATAAAGCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGC  >  NC_000913/1121921‑1122191

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