Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 501,526 A→G intergenic (+97/+36) gsk → / ← ybaL inosine/guanosine kinase/putative transporter YbaL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5260AG100.0% 154.4 / NA 41intergenic (+97/+36)gsk/ybaLinosine/guanosine kinase/putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/14);  total (27/14)

GACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTT  >  NC_000913/501397‑501662
                                                                                                                                 |                                                                                                                                        
gACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTc                                                                                                                                <  2:682893/140‑1 (MQ=255)
 aCAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCa                                                                                                                               <  1:1016880/140‑1 (MQ=255)
 aCAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCa                                                                                                                               <  1:1734591/140‑1 (MQ=255)
 aCAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCa                                                                                                                               <  2:1743400/140‑1 (MQ=255)
    gCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCatat                                                                                                                            >  2:1475065/1‑140 (MQ=255)
    gCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCatat                                                                                                                            >  1:917470/1‑140 (MQ=255)
     ccTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATAtt                                                                                                                           >  1:447094/1‑140 (MQ=255)
         gAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTAc                                                                                                                       <  1:603637/140‑1 (MQ=255)
         gAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTAc                                                                                                                       <  2:1376273/140‑1 (MQ=255)
               tCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCgg                                                                                                                 >  2:294914/1‑140 (MQ=255)
                cTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGc                                                                                                                <  1:441115/140‑1 (MQ=255)
                    cTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAacc                                                                                                            <  1:294914/140‑1 (MQ=255)
                     tGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAacca                                                                                                           <  1:318756/140‑1 (MQ=255)
                         aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATcc                                                                                                                   >  2:174739/1‑128 (MQ=255)
                         aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATcc                                                                                                                   <  1:174739/128‑1 (MQ=255)
                         aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGt                                                                                                       >  1:94422/1‑140 (MQ=255)
                         aTCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGt                                                                                                       >  2:1018486/1‑140 (MQ=255)
                                      cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGt                                                                                              >  1:391707/1‑136 (MQ=255)
                                      cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGt                                                                                              <  2:391707/136‑1 (MQ=255)
                                      cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacc                                                                                          >  1:719554/1‑140 (MQ=255)
                                            ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                    >  1:1053844/1‑140 (MQ=255)
                                            ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                    >  1:897804/1‑140 (MQ=255)
                                                 aGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgc                                                                               >  2:1094725/1‑140 (MQ=255)
                                                        aTAGGGCTTTCCCGCCGCATCTGGCTTTTGGCTCTCGATGCCTGAGGCGACGCGGGCGCGTCTTATCATGCCTGCATTTTTTTAATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCACCACCTCACCCGCCGGTGGc                                                                        <  2:1146220/140‑1 (MQ=255)
                                                        aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                        >  2:1088953/1‑140 (MQ=255)
                                                        aTAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGc                                                                        >  2:548671/1‑140 (MQ=255)
                                                             gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                   >  1:1217292/1‑140 (MQ=255)
                                                             gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                   >  1:771388/1‑140 (MQ=255)
                                                             gCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                   >  1:359526/1‑140 (MQ=255)
                                                                ttCCCGCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGTGTTTTATCATGCCTGAATTTTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCcag                                                                <  2:359526/140‑1 (MQ=255)
                                                                tcccccccgcTTCTGGCTTTGGCCTCCCGATGCCTGATGCGCCGCTGGCGCTTCTTATCATCCCTGCATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCcag                                                                <  2:447094/134‑1 (MQ=255)
                                                                 tcccgccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCGGGCGCGTCTTATCATGCCTGCATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                               <  2:94422/137‑1 (MQ=37)
                                                                 tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                               >  1:722468/1‑140 (MQ=255)
                                                                 tCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                               >  1:810330/1‑140 (MQ=255)
                                                                 ccACCCCGCATCGGGTATTTGGCTCTCGTTGCCTGATGCGCCGCGGGCGCGTTTTATCATGCCTGCATTTTTTTCATTTTTTCCATCCGCCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCcagc                                                               <  2:964514/139‑1 (MQ=255)
                                                                  cACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                              >  1:1253843/1‑140 (MQ=255)
                                                                      ccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc                                                          >  1:1166799/1‑140 (MQ=255)
                                                                       cgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCGGGCGCTTTTTATCATGCCTGCATTTTTTTCATTTTTTACATCCGGCAACCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTcc                                                         <  2:922992/140‑1 (MQ=255)
                                                                       cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                         >  2:1468792/1‑140 (MQ=255)
                                                                       cgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTcc                                                         >  1:1315940/1‑140 (MQ=255)
                                                                                  tttGCCTCTCGATGCCGGATGCCCCGCTGGCGCTTCTTATCATGCCTGCATATTTTTCATATTTTCCACCCGGCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACg                                              <  2:719554/140‑1 (MQ=255)
                                                                                    tGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTGCATTTTTTTCATATTTTACACCCGGCAACCACCTTTTACCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACggg                                            <  1:1468792/140‑1 (MQ=39)
                                                                                             ttGCCTGAGGCGCCGCGGGCGCGTCTTATCATGCCTGCATTTTTTTCATTTTTTACATCCGGCACCCACCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCAc                                   <  2:917470/139‑1 (MQ=255)
                                                                                               cccTGATGCGACGCGGGCGCGTCTTATCATGCCTGCTTTTTTTTAATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                 <  2:1694711/139‑1 (MQ=255)
                                                                                                ccggaggCGACGCGGGCGCGTTTTATCATGCCTGAATTTTTTTCATATTTTACATCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                <  2:988011/134‑1 (MQ=255)
                                                                                                  tgttGCGACGCTGGCGCTTCTTATCATGCCTCCATTTTTTTCATATTTTACACCCGCCAACCCCCTTTCCCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                              <  2:1053844/137‑1 (MQ=255)
                                                                                                  tGATGCGCCGCTGGCGCGTTTTATCATGCCTGCATTTTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                              <  2:897804/140‑1 (MQ=255)
                                                                                                  tGATGCGACGCTGGCGCGTCTTATCATGCCTGCATTTTTTTCATATTTTACACCCGGCACCCACCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCg                              <  1:1094725/140‑1 (MQ=255)
                                                                                                   gATGCGCCGCTGGCGCGTCTTATCATGCCTGCATTTTTTTCATTTTTAACCTCCGCCAACCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                             <  2:322392/140‑1 (MQ=255)
                                                                                                   gATGCGCCGCTGGCGCGTCTTATCATGCCTGCATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCGCGCTCGc                             <  1:922579/140‑1 (MQ=255)
                                                                                                    aTGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                            >  1:670421/1‑140 (MQ=255)
                                                                                                    aTGCGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                            >  1:821409/1‑140 (MQ=255)
                                                                                                     ggcggcgCGGGCGCGTCTTATCATGCCTGCATTTTTTTCATATTTTACACCCGACACCCCCCGTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGccc                           <  2:930998/139‑1 (MQ=255)
                                                                                                       cGACGCTGGCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCAt                         >  1:988011/1‑140 (MQ=255)
                                                                                                         gcgcgGGCGCGTTTTATCATCCCTGCATTTTTTTCATTTTTTACACCCGCCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                       <  2:1315940/139‑1 (MQ=255)
                                                                                                             gggCCCGTCTTATCTTGCCTGCTTTTTTTTCATATTTTACACCCGGAACCCACCTTTACCCCCGCCCCCGCCCCGCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                   <  2:1253843/139‑1 (MQ=255)
                                                                                                              ggCGCGTCTTATCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAcc                  >  1:679733/1‑140 (MQ=255)
                                                                                                                         tCATGCCTGCATATTTTTCATATTTTACATCCGGCAACCACCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGcacc       <  1:1544509/140‑1 (MQ=255)
                                                                                                                            tcccTGCATTTTTTTCATTTTTTACATCCGGCACCCCCCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg    <  2:821409/138‑1 (MQ=255)
                                                                                                                             gCCTGCATTTTTTTCATTTTTTACATCCGGCACCCACCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt   <  1:1670457/140‑1 (MQ=255)
                                                                                                                              ccTGCATTTTTTTCATTTTTTACATCCGGCACCCACTGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt  <  1:1088953/140‑1 (MQ=255)
                                                                                                                              ccTGCATATTTTTCATATTTTACACCCGGCACCCACCGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt  <  1:548671/140‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                        
GACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTT  >  NC_000913/501397‑501662

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