Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,158,773 A→G A49A (GCA→GCG yqhH → lipoprotein YqhH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,158,7730AG100.0% 213.6 / NA 62A49A (GCA→GCGyqhHlipoprotein YqhH
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/33);  total (29/33)

ATTTTCACCGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCACTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATCCCCGGT  >  NC_000913/3158636‑3158900
                                                                                                                                         |                                                                                                                               
aTTTTCACCGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCt                                                                                                                               >  1:84306/1‑140 (MQ=255)
aTTTTCACCGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCt                                                                                                                               >  2:962515/1‑140 (MQ=255)
  tttCACCGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTAc                                                                                                                             <  2:894302/140‑1 (MQ=255)
     cACCGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGcc                                                                                                                          <  2:1419540/140‑1 (MQ=255)
       ccGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCca                                                                                                                        <  2:1194179/140‑1 (MQ=255)
       ccGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCca                                                                                                                        <  1:347611/140‑1 (MQ=255)
          tGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACaa                                                                                                                     >  1:1395984/1‑140 (MQ=255)
                  gttgttCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTAtgc                                                                                                             >  1:1243891/1‑140 (MQ=255)
                  gttgttCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTAtgc                                                                                                             >  1:1421741/1‑140 (MQ=255)
                   ttgttCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTAtgct                                                                                                            >  1:767444/1‑140 (MQ=255)
                          ggCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAACTGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAAt                                                                                                     <  1:1290172/140‑1 (MQ=255)
                          ggCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAAt                                                                                                     <  1:1259986/140‑1 (MQ=255)
                                cTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAg                                                                                               >  2:453743/1‑140 (MQ=255)
                                cTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAg                                                                                               >  2:762155/1‑140 (MQ=255)
                                    ttgctCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTaa                                                                                           <  2:1383284/140‑1 (MQ=255)
                                     tgctCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAAc                                                                                          >  1:812066/1‑140 (MQ=255)
                                      gctCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACa                                                                                         <  1:325247/140‑1 (MQ=255)
                                      gctCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACa                                                                                         <  2:812066/140‑1 (MQ=255)
                                        tCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACaga                                                                                       <  2:1421741/140‑1 (MQ=255)
                                          aGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGc                                                                                     >  1:280571/1‑140 (MQ=255)
                                          aGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGc                                                                                     >  1:1250310/1‑140 (MQ=255)
                                          aGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGc                                                                                     >  1:1066235/1‑140 (MQ=255)
                                            tGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGa                                                                                                 >  2:986618/1‑126 (MQ=255)
                                            tGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGa                                                                                                 <  1:986618/126‑1 (MQ=255)
                                            tGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCAGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGa                                                                                                 >  2:675487/1‑126 (MQ=255)
                                            tGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCAGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGa                                                                                                 <  1:675487/126‑1 (MQ=255)
                                                 gCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATAcg                                                                              >  1:846435/1‑140 (MQ=255)
                                                 gCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATAcg                                                                              <  1:88245/140‑1 (MQ=255)
                                                  cGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAAc                                                                                          <  2:744920/127‑1 (MQ=255)
                                                  cGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAAc                                                                                          >  1:744920/1‑127 (MQ=255)
                                                  cGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATAcgc                                                                             >  2:512138/1‑140 (MQ=255)
                                                      aaTGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCt                                                                         >  1:822759/1‑140 (MQ=255)
                                                                  gTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGa                                                             >  2:465414/1‑140 (MQ=255)
                                                                      aTCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGa                                                             >  1:224748/1‑136 (MQ=255)
                                                                      aTCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGa                                                             <  2:224748/136‑1 (MQ=255)
                                                                          gCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTg                                                     <  2:1250310/140‑1 (MQ=255)
                                                                          gCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTg                                                     <  2:822759/140‑1 (MQ=255)
                                                                          gCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTg                                                     <  2:1424374/140‑1 (MQ=255)
                                                                          gCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTg                                                     <  2:1066235/140‑1 (MQ=255)
                                                                              gCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTg                                                 <  1:962515/140‑1 (MQ=255)
                                                                                 aGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCt                                              >  1:1064742/1‑140 (MQ=255)
                                                                                 aGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCt                                              >  1:1478006/1‑140 (MQ=255)
                                                                                    aaTGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTgcg                                           <  2:481992/140‑1 (MQ=255)
                                                                                    aaTGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTgcg                                           <  1:512138/140‑1 (MQ=255)
                                                                                        tgCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCgctgc                                       >  2:527834/1‑140 (MQ=255)
                                                                                        tgCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCgctgc                                       >  2:831653/1‑140 (MQ=255)
                                                                                        tgCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCgctgc                                       >  2:67585/1‑140 (MQ=255)
                                                                                        tgCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCgctgc                                       >  2:714096/1‑140 (MQ=255)
                                                                                        tgCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCgctgc                                       <  1:527834/140‑1 (MQ=255)
                                                                                         gCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCgctgct                                      <  2:280571/140‑1 (MQ=255)
                                                                                             aCATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCg                                  <  2:846435/140‑1 (MQ=255)
                                                                                              cATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGt                                 >  1:1258069/1‑140 (MQ=255)
                                                                                                ttAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTAt                               <  2:1478006/140‑1 (MQ=255)
                                                                                                       ccAAAATCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCa                        >  2:1324911/1‑140 (MQ=255)
                                                                                                         aaaaTCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTAt                               >  2:838955/1‑131 (MQ=255)
                                                                                                         aaaaTCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTAt                               <  1:838955/131‑1 (MQ=255)
                                                                                                          aaaTCGCCCGGCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGaa                     <  2:1064742/140‑1 (MQ=255)
                                                                                                                    gCTTGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATaaaaaa           <  1:1324911/140‑1 (MQ=255)
                                                                                                                      ttGAGCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAAt         <  1:762155/140‑1 (MQ=255)
                                                                                                                          gCAGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATcccc     <  1:1178648/140‑1 (MQ=255)
                                                                                                                            aGGATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATCCCCgg   <  2:1243891/140‑1 (MQ=255)
                                                                                                                             ggATATGAAAGCGCTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATCCCCGGt  <  2:84306/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                               
ATTTTCACCGTGGGAGCTGTTGTTCTGGCAACCTGCTTGCTCAGTGGCTGCGTCAATGAGCAAAAGGTCAATCAGCTGGCGAGCAATGTGCAAACATTAAATGCCAAAATCGCCCGGCTTGAGCAGGATATGAAAGCACTACGCCCACAAATCTATGCTGCCAAATCCGAAGCTAACAGAGCCAATACGCGTCTTGATGCTCAGGACTATTTTGATTGCCTGCGCTGCTTGCGTATGTACGCAGAATGATAAAAAAATCCCCGGT  >  NC_000913/3158636‑3158900

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: