Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,574,259 A→G pseudogene (341/2937 nt) yjiV → pseudogene, conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,574,2590AG100.0% 269.2 / NA 69pseudogene (341/2937 nt)yjiVpseudogene, conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/32);  total (37/32)

TGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCA  >  NC_000913/4574134‑4574392
                                                                                                                             |                                                                                                                                     
tgtgtatgtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGc                                                                                                                          <  2:85002/139‑1 (MQ=255)
      tgtgAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGaa                                                                                                                    >  1:255482/1‑139 (MQ=255)
          aaCCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                >  1:552982/1‑139 (MQ=255)
          aaCCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                >  1:992748/1‑139 (MQ=255)
          aaCCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                >  1:175569/1‑139 (MQ=255)
            ccAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                              >  1:569997/1‑139 (MQ=255)
            ccAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTctct                                                                                                              >  2:118151/1‑139 (MQ=255)
                 gACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                         <  2:255482/139‑1 (MQ=255)
                 gACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCa                                                                                                         <  2:552982/139‑1 (MQ=255)
                   cAAAATTGTGACTGTGTAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGTCGATGAATTCTCTTTGTAgg                                                                                                       <  1:750871/139‑1 (MQ=255)
                          gtgACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGt                                                                                                <  2:175569/139‑1 (MQ=255)
                          gtgACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGt                                                                                                <  1:118151/139‑1 (MQ=255)
                           tgACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTa                                                                                               <  1:706476/139‑1 (MQ=255)
                                           gTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAtt                                                                                                                  <  2:768905/104‑1 (MQ=255)
                                           gTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAtt                                                                                                                  >  1:768905/1‑104 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAt                                                                                                                   <  2:893735/92‑1 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAAt                                                                                                                   >  1:893735/1‑92 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTAc                                                                                              <  2:1085547/113‑1 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTAc                                                                                              >  1:1085547/1‑113 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAgaagaa                                                                    >  2:275242/1‑139 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAgaagaa                                                                    >  1:605760/1‑139 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAgaagaa                                                                    >  1:482384/1‑139 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAgaagaa                                                                    >  2:478313/1‑139 (MQ=255)
                                                      tGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAgaagaa                                                                    >  2:674121/1‑139 (MQ=255)
                                                        cATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGa                                                                                                      <  2:1104695/103‑1 (MQ=255)
                                                        cATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGa                                                                                                      >  1:1104695/1‑103 (MQ=255)
                                                        cATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACAt                                                                                   <  2:886757/122‑1 (MQ=255)
                                                        cATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACAt                                                                                   >  1:886757/1‑122 (MQ=255)
                                                               tAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGa                                                                                                                     <  2:222387/81‑1 (MQ=255)
                                                               tAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGa                                                                                                                     >  1:222387/1‑81 (MQ=255)
                                                                       gTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGaa                                                   >  2:423886/1‑139 (MQ=255)
                                                                       gTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGaa                                                   >  2:139858/1‑139 (MQ=255)
                                                                         aaTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGaaaa                                                                  >  2:832743/1‑122 (MQ=255)
                                                                         aaTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGaaaa                                                                  <  1:832743/122‑1 (MQ=255)
                                                                         aaTGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAGAAGAGAACa                                                 >  1:707655/1‑139 (MQ=255)
                                                                              gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTAc                                                                                              <  2:534855/89‑1 (MQ=255)
                                                                              gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTAc                                                                                              >  1:534855/1‑89 (MQ=255)
                                                                              gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGa                                                                                 <  1:853987/102‑1 (MQ=255)
                                                                              gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGa                                                                                 >  2:853987/1‑102 (MQ=255)
                                                                              gCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGtttt                                            >  2:481113/1‑139 (MQ=255)
                                                                               ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                <  2:753819/70‑1 (MQ=255)
                                                                               ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATtct                                                                                                                >  1:753819/1‑70 (MQ=255)
                                                                               ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTAc                                                                           <  2:74847/107‑1 (MQ=255)
                                                                               ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTAc                                                                           >  1:74847/1‑107 (MQ=255)
                                                                               ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTc                                           >  1:850323/1‑139 (MQ=255)
                                                                               ccTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTc                                           >  1:228665/1‑139 (MQ=255)
                                                                                  cATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGc                                        <  1:674121/139‑1 (MQ=255)
                                                                                  cATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGc                                        >  2:945975/1‑139 (MQ=255)
                                                                                  cATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGc                                        <  1:139858/139‑1 (MQ=255)
                                                                                              cGGGACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGatta                            >  2:427144/1‑139 (MQ=255)
                                                                                               gggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCtttt                                                                                                           <  2:690767/59‑2 (MQ=255)
                                                                                               gggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTg                                                                                           >  2:928694/1‑75 (MQ=255)
                                                                                               gggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTg                                                                                           <  1:928694/75‑1 (MQ=255)
                                                                                               gggACAAAAAGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGattat                           <  2:707655/139‑1 (MQ=255)
                                                                                               gggACAAAAAGGGCAAAATAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTg                                                                                                           >  1:690767/1‑59 (MQ=255)
                                                                                                 gACAAAAAGGGCAAGCTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTATTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTg                         >  1:671880/1‑139 (MQ=255)
                                                                                                    aaaaaGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCg                      <  2:992748/139‑1 (MQ=255)
                                                                                                    aaaaaGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCg                      <  1:423886/139‑1 (MQ=255)
                                                                                                      aaaGGGCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCa                    <  1:945975/139‑1 (MQ=255)
                                                                                                           gCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCg                      <  2:497329/132‑1 (MQ=255)
                                                                                                           gCAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCg                      >  1:497329/1‑132 (MQ=255)
                                                                                                            cAAACTAGTCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGgagag              <  2:850323/139‑1 (MQ=255)
                                                                                                                    tCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAg                                                                                                        <  2:1040181/41‑1 (MQ=255)
                                                                                                                    tCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAg                                                                                                        >  1:1040181/1‑41 (MQ=255)
                                                                                                                    tCCAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGa      >  2:916995/1‑139 (MQ=255)
                                                                                                                     ccAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAg     <  1:427144/139‑1 (MQ=255)
                                                                                                                      cAGTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAgg    <  1:478313/139‑1 (MQ=255)
                                                                                                                       aGTGGAGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGACAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGc   <  2:794612/139‑1 (MQ=255)
                                                                                                                        gTGGGGAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCa  >  1:744622/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                                     
TGTGTATGTGAACCAGCGACAAAATTGTGACTGTGGAAGCCCTGTATACGAAGTTGCATTCTGTAATGATTGTAATGAGCCTCATCTTCTGGCACGGGACAAAAAGGGCAAACTAGTCCAGTGGGAAAATAAAGGTGGCGATGAATTCTCTTTGCAGGATGAAGTACCTGTTGAACATGACGCTACAGAAGAAAAAGTCGAAAAAGAGAACAGTTTTCAGCCTCCGTTGATTATTGCCGCAGGAGAGACCAGCGAGGCA  >  NC_000913/4574134‑4574392

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