Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 243,735 G→T E65* (GAG→TAG)  lpcA → D‑sedoheptulose 7‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913243,7350GT100.0% 184.8 / NA 52E65* (GAG→TAG) lpcAD‑sedoheptulose 7‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/35);  total (17/35)

TTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAG  >  NC_000913/243602‑243867
                                                                                                                                      |                                                                                                                                     
tttAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGa                                                                                                                                 <  2:1458869/140‑1 (MQ=255)
  tAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGAcc                                                                                                                               <  1:1029478/140‑1 (MQ=255)
     aaGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGt                                                                                                                            <  1:668113/140‑1 (MQ=255)
           aCGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTAc                                                                                                                      <  1:1050905/140‑1 (MQ=255)
                  tatTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTTCCCTTGCGACGCTATGCACTTTTCCGAATAGTTGACCGGTCGCTACCCTGgaa                                                                                                              >  1:268107/1‑140 (MQ=255)
                  tatTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCCTTTCGACGCTATTCACTTTTCCCGATAGTTGACCGGGCGCTACCCTGaaa                                                                                                              >  2:763648/1‑140 (MQ=255)
                       aCGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGt                                                                                                         <  1:626212/140‑1 (MQ=255)
                       aCGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGt                                                                                                         <  1:349477/140‑1 (MQ=255)
                            aTTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCggg                                                                                                    <  2:621822/140‑1 (MQ=255)
                            aTTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGc                                                                                                   <  1:429628/140‑1 (MQ=255)
                             ttCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTTCCCTTTCGAAGCTATGCACTTTTCCCAATAATTTACCGGGCGCTACCCTGAAAAACCTCCcggg                                                                                                   >  2:1049935/1‑139 (MQ=255)
                              tCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTTCCCTTTCGGCGCCATGCACTTTTCCCGAAAAGTGACCCGGCGCTACCCTGGAAAACCGCCCGGGc                                                                                                   >  2:257226/1‑140 (MQ=255)
                                  cgcgcgGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTAccc                                                                                              <  2:1463756/140‑1 (MQ=255)
                                      cggcggTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCg                                                                                          <  1:1487582/140‑1 (MQ=255)
                                      cggcggTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCg                                                                                          <  2:492693/140‑1 (MQ=255)
                                      cggcggTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCg                                                                                          <  2:1425190/140‑1 (MQ=255)
                                       ggcggTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGa                                                                                         <  2:1356048/140‑1 (MQ=255)
                                           gTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGc                                                                                     <  1:1464497/140‑1 (MQ=255)
                                           gTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGc                                                                                     <  2:1052076/140‑1 (MQ=255)
                                           gTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGc                                                                                     <  2:1470700/140‑1 (MQ=255)
                                                     cAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGAc                                                                           <  1:6050/140‑1 (MQ=255)
                                                      aGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCCCTTCGACGCCATACACTTTTCCGGATAGTTGACCGGTCGCTACCGCGAAAAACCTCCGGGCGACCCCCCGATTGCTAttttttccg                                                                          >  2:919915/1‑134 (MQ=255)
                                                       gACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGt                                                                         <  1:1099718/140‑1 (MQ=255)
                                                       gACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGt                                                                         <  1:782298/140‑1 (MQ=255)
                                                         cAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTa                                                                       <  1:1132420/140‑1 (MQ=255)
                                                              ttAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCat                                                                  <  1:1075846/140‑1 (MQ=255)
                                                                      ggTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTTCCCTTTCGGCGCTCTGCACCTTTCCGAATAATTTTCCGGGCGCCCCCCTGCAAAACCCCCCGGCCACCCCGCGGTTGCTATTTCTTACGCTTGTCATAttttccc                                                          >  2:951427/1‑136 (MQ=255)
                                                                      ggTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTTCCCTTGCGACGCCATGCACTTTTCCCCATAGTTGACCGGGCGCTACCCCGAAAACCCTCCGGGCCACCCCGCGATTGCTATTTTTTACGTTAGTCATATTTTCtt                                                          >  2:279019/1‑139 (MQ=255)
                                                                      ggTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCCTTGCGACGCTATGCACTTTGCCGAATAGATGACCGGGCGCTACCGTGAAAAACGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTTCtt                                                          >  2:813542/1‑139 (MQ=255)
                                                                      ggTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTTCCGAATAGTTGACCGGGCGCTACCGTGGAAACCGTCCGGGCTACCCCGCGGTTGCTATTTTTTACGTTAGTCATATTTCCTg                                                          >  2:46107/1‑140 (MQ=255)
                                                                         ggCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGt                                                       <  2:430692/140‑1 (MQ=255)
                                                                         ggCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGt                                                       <  1:951427/140‑1 (MQ=255)
                                                                                     tCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGAttt                                           <  2:1108926/140‑1 (MQ=255)
                                                                                        tGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGt                                       <  1:456584/140‑1 (MQ=255)
                                                                                         gCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGt                                       <  1:1049935/140‑1 (MQ=255)
                                                                                           ggCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGttt                                     <  1:871693/140‑1 (MQ=255)
                                                                                                cggcggTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATg                                <  1:887069/140‑1 (MQ=255)
                                                                                                   cggTTCCCATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGata                             <  1:1332705/140‑1 (MQ=255)
                                                                                                        cccATTGCGACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTtc                        <  1:208659/140‑1 (MQ=255)
                                                                                                                gACGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTa                <  2:488133/140‑1 (MQ=255)
                                                                                                                 aCGCTATGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTAc               <  2:1259131/140‑1 (MQ=255)
                                                                                                                       tGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGaa         >  2:1294957/1‑140 (MQ=255)
                                                                                                                       tGCACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGaa         >  2:897002/1‑140 (MQ=255)
                                                                                                                         cACTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGc       <  1:279260/140‑1 (MQ=255)
                                                                                                                           cTTTGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCgg     >  2:1180572/1‑140 (MQ=255)
                                                                                                                              ttGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTa   >  1:916574/1‑140 (MQ=255)
                                                                                                                              ttGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTa   >  1:1091834/1‑140 (MQ=255)
                                                                                                                              ttGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTa   >  1:332831/1‑140 (MQ=255)
                                                                                                                              ttGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTa   >  2:1336243/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCg                   >  1:722753/1‑123 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCg                   <  2:722753/123‑1 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTa   <  1:373765/139‑1 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTa   >  2:373765/1‑139 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  2:399336/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  2:39314/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  1:459380/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  1:1369205/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  2:785738/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  1:739425/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  1:848790/1‑140 (MQ=255)
                                                                                                                               tGCCGAATAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAg  >  2:1210410/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TTTAAAAGATGACGCCAATATTCACGCCATTCAGCGCGCGGCGGTCCTGTTAGCAGACAGCTTTAAAGCCGGTGGCAAAGTGCTTTCCTGCGGCAACGGCGGTTCCCATTGCGACGCTATGCACTTTGCCGAAGAGTTGACCGGTCGCTACCGTGAAAACCGTCCGGGCTACCCGGCGATTGCTATTTCTGACGTTAGTCATATTTCCTGCGTCGGTAATGATTTCGGTTTCAATGATATTTTCTCCCGCTACGTTGAAGCGGTAG  >  NC_000913/243602‑243867

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