Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,941,734 Δ1 bp noncoding (409/2906 nt) rrlC → 23S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,941,7340A.100.0% 152.1 / NA 87noncoding (409/2906 nt)rrlC23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (55/32);  total (55/32)

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACC  >  NC_000913/3941644‑3941771
                                                                                          |                                     
gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                     >  2:5444263/1‑93 (MQ=25)
gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                     <  2:3925993/93‑1 (MQ=18)
gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                     >  2:3943968/1‑93 (MQ=25)
 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggtgga                                    >  2:2080174/1‑92 (MQ=255)
 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                    >  2:1739682/1‑92 (MQ=18)
  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATggggggcc                                   >  2:650014/1‑91 (MQ=16)
  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGAc                                  >  1:1635746/1‑93 (MQ=9)
  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGAc                                  >  1:3642175/1‑93 (MQ=16)
   aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGAcc                                 >  1:2721150/1‑93 (MQ=16)
   aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGAcc                                 <  2:4752217/93‑1 (MQ=9)
   aCAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGAcc                                 <  2:4831469/93‑1 (MQ=9)
    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGAcct                                >  1:4670303/1‑92 (MQ=9)
    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCa                                <  2:1228711/93‑1 (MQ=9)
    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCa                                >  1:3855490/1‑93 (MQ=16)
     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGTGACCAt                               >  1:2281108/1‑93 (MQ=255)
     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAt                               >  2:509263/1‑93 (MQ=16)
     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAt                               >  1:2105097/1‑93 (MQ=9)
      gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtc                              >  1:5796050/1‑93 (MQ=9)
       ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcc                             <  1:2679284/93‑1 (MQ=9)
        gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcct                            >  2:3081780/1‑93 (MQ=12)
        gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcct                            <  2:2197471/93‑1 (MQ=9)
         tGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcctc                           <  2:2118006/93‑1 (MQ=9)
         tGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcctc                           >  2:1217093/1‑93 (MQ=16)
          gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcctcc                          >  1:3813554/1‑93 (MQ=16)
          gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcctcc                          >  2:376510/1‑93 (MQ=16)
          gACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCAtcctcc                          >  1:2642254/1‑93 (MQ=9)
           aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCa                         >  1:1436404/1‑93 (MQ=16)
           aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCa                         <  2:3285936/93‑1 (MQ=9)
           aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCa                         >  2:1951747/1‑93 (MQ=16)
           aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCa                         >  2:5143058/1‑93 (MQ=16)
           aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCa                         >  1:2677405/1‑93 (MQ=16)
             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGGACATCCTCCAAg                       >  1:3687349/1‑93 (MQ=255)
             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAg                       >  1:2859391/1‑93 (MQ=9)
             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAg                       >  2:1952039/1‑93 (MQ=12)
             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAg                       <  1:6454628/93‑1 (MQ=9)
              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAgg                      >  2:3722970/1‑93 (MQ=16)
              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAgg                      >  2:2425379/1‑93 (MQ=16)
              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAgg                      >  2:1147435/1‑93 (MQ=12)
                cccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCt                    <  1:4496214/93‑1 (MQ=9)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGGCATGGGGGATCCTGTCTGAAGATGCGCGGACCATCCTCCAACGCTa                   >  2:3676834/1‑93 (MQ=11)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTa                   >  2:2638971/1‑93 (MQ=16)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTa                   >  1:4516523/1‑93 (MQ=16)
                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTa                   >  1:5438745/1‑93 (MQ=16)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTaa                  <  2:1254097/93‑1 (MQ=9)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTaa                  <  1:2951470/93‑1 (MQ=9)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTaa                  >  1:2228906/1‑93 (MQ=16)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTaa                  <  2:460033/93‑1 (MQ=9)
                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTaa                  <  1:620870/93‑1 (MQ=9)
                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAAt                >  1:4508164/1‑93 (MQ=16)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATa               <  1:3950232/93‑1 (MQ=9)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATa               >  2:550962/1‑93 (MQ=16)
                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATa               >  2:4379204/1‑93 (MQ=16)
                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATAc              <  1:1371089/93‑1 (MQ=9)
                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGGATCCTGTCTGAATATGGGGGGACCATCCTACAAGGCTAAATAc              >  1:243237/1‑93 (MQ=11)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACt             >  2:2251003/1‑93 (MQ=16)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACt             <  2:3914379/93‑1 (MQ=9)
                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACt             >  2:3684428/1‑93 (MQ=16)
                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTc            >  2:6149742/1‑93 (MQ=16)
                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTc            <  1:1692246/93‑1 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTcc           >  2:507600/1‑93 (MQ=16)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTcc           >  1:3913146/1‑93 (MQ=9)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTcc           >  1:3259028/1‑93 (MQ=16)
                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTcc           >  1:465521/1‑93 (MQ=16)
                          aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCct          >  2:165387/1‑93 (MQ=9)
                          aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCct          <  1:6461510/93‑1 (MQ=9)
                           aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctg         >  2:1927351/1‑93 (MQ=12)
                           aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctg         >  1:5668435/1‑93 (MQ=16)
                           aaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctg         >  1:3145304/1‑93 (MQ=16)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctga        <  1:1896434/93‑1 (MQ=9)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctga        <  2:2712284/93‑1 (MQ=9)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctga        >  1:3358629/1‑93 (MQ=16)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctga        <  1:4106435/93‑1 (MQ=9)
                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctga        >  2:1467211/1‑93 (MQ=12)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       >  2:4692908/1‑93 (MQ=16)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       <  1:4890212/93‑1 (MQ=9)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       >  2:5569532/1‑93 (MQ=16)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       >  2:3582865/1‑93 (MQ=9)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       >  2:283009/1‑93 (MQ=16)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       <  1:4191204/93‑1 (MQ=9)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       <  1:1554983/93‑1 (MQ=9)
                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgac       <  1:3241528/93‑1 (MQ=9)
                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgact      >  2:818404/1‑93 (MQ=9)
                               gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgactg     <  2:125972/93‑1 (MQ=9)
                                cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga    >  1:1707803/1‑93 (MQ=9)
                                cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga    <  2:5434741/93‑1 (MQ=9)
                                cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga    <  2:5718154/93‑1 (MQ=9)
                                cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCctgactga    >  1:5891770/1‑93 (MQ=12)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac   <  2:24998/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac   >  2:5482519/1‑93 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac   <  2:1697952/93‑1 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac   >  2:5736602/1‑93 (MQ=9)
                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCtgactgac   <  1:4241133/93‑1 (MQ=9)
                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc  <  1:1585101/93‑1 (MQ=0)
                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc  >  2:5654280/1‑93 (MQ=0)
                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGGGGACCATCCTCCAAGCCTAAATACTCCTGACTGAcc  >  2:1225976/1‑93 (MQ=0)
                                                                                          |                                     
GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACC  >  NC_000913/3941644‑3941771

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: