Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,022,355 A→G L510S (TTA→TCA)  yccS ← putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,022,3550AG100.0% 215.8 / NA 59L510S (TTA→TCA) yccSputative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/32);  total (27/32)

CCAGTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTAACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCA  >  NC_000913/1022209‑1022499
                                                                                                                                                  |                                                                                                                                                
ccAGTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGAc                                                                                                                                                <  2:914097/149‑1 (MQ=255)
ccAGTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGc                                                                                                                                              <  1:211294/151‑1 (MQ=255)
 cAGTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCt                                                                                                                                             <  1:176093/151‑1 (MQ=255)
   gTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCGCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTAc                                                                                                                                           <  2:749424/151‑1 (MQ=255)
   gTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTAc                                                                                                                                           >  2:214759/1‑151 (MQ=255)
   gTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTAc                                                                                                                                           >  1:857085/1‑151 (MQ=255)
     aTGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACtt                                                                                                                                         >  2:608384/1‑151 (MQ=255)
     aTGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACtt                                                                                                                                         >  1:297426/1‑151 (MQ=255)
       ggCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCa                                                                                                                                       <  1:884459/151‑1 (MQ=255)
         cATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCaa                                                                                                                                      >  2:520521/1‑150 (MQ=255)
             gagaTACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCt                                                                                                                                 >  1:322402/1‑151 (MQ=255)
                aTACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTc                                                                                                                               >  1:305285/1‑150 (MQ=255)
                aTACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCa                                                                                                                              <  1:88942/151‑1 (MQ=255)
                 tACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCAc                                                                                                                             >  2:735316/1‑151 (MQ=255)
                   ccGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCAccc                                                                                                                           >  1:675745/1‑151 (MQ=255)
                        aGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTa                                                                                                                       >  2:577297/1‑150 (MQ=255)
                             gCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTaa                                                                                                                 <  2:310091/151‑1 (MQ=255)
                              ccTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTaaa                                                                                                                <  2:150651/151‑1 (MQ=255)
                               ctctGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTaaa                                                                                                                <  2:378144/150‑1 (MQ=255)
                                       cgcgTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAg                                                                                                       <  2:279010/151‑1 (MQ=255)
                                           ttCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCAcc                                                                                                   >  2:321031/1‑151 (MQ=255)
                                              gAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCaa                                                                                                 >  2:508117/1‑150 (MQ=255)
                                               aGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCaaa                                                                                                >  1:157739/1‑150 (MQ=255)
                                               aGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCaaa                                                                                                >  1:535811/1‑150 (MQ=255)
                                                        ggCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGtgatga                                                                                      <  1:508019/151‑1 (MQ=255)
                                                          cAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAAc                                                                                    >  1:793718/1‑151 (MQ=255)
                                                           aGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACa                                                                                   <  1:657673/151‑1 (MQ=255)
                                                            gATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACAt                                                                                  >  1:532049/1‑151 (MQ=255)
                                                                gCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGtt                                                                              <  1:94783/151‑1 (MQ=255)
                                                                 cGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTg                                                                             <  2:581116/151‑1 (MQ=255)
                                                                  gAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGc                                                                            >  2:660525/1‑151 (MQ=255)
                                                                        gCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGc                                                                      <  1:303093/151‑1 (MQ=255)
                                                                          cAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGt                                                                    <  1:1058996/151‑1 (MQ=255)
                                                                           aGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGt                                                                    >  2:693456/1‑150 (MQ=255)
                                                                              cAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTg                                                                <  2:185783/151‑1 (MQ=255)
                                                                               aGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGc                                                               >  1:402152/1‑151 (MQ=255)
                                                                                 gCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCAc                                                             <  2:260481/151‑1 (MQ=255)
                                                                                           gCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAg                                                   <  2:808062/151‑1 (MQ=255)
                                                                                              cACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGc                                                 <  1:993385/150‑1 (MQ=255)
                                                                                               aCTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCa                                                >  1:911905/1‑150 (MQ=255)
                                                                                                  gccgccCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCaaaaa                                            <  1:14417/151‑1 (MQ=255)
                                                                                                  gccgccCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCaaaa                                             >  1:1017061/1‑150 (MQ=255)
                                                                                                  gccgccCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCaaaa                                             >  2:267883/1‑150 (MQ=255)
                                                                                                       ccACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGc                                       >  2:654498/1‑151 (MQ=255)
                                                                                                           gCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGc                                   <  2:467411/151‑1 (MQ=255)
                                                                                                            cAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGcc                                  <  2:50535/151‑1 (MQ=255)
                                                                                                                 gCACAACCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGtaa                             <  2:254394/151‑1 (MQ=255)
                                                                                                                      aCCAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAAcc                        <  1:1043588/151‑1 (MQ=255)
                                                                                                                       ccAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATCAACAt                                                                                  >  2:364820/1‑92 (MQ=255)
                                                                                                                       ccAATCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATCAACAt                                                                                  <  1:364831/92‑1 (MQ=255)
                                                                                                                         aaTCAGCGTATCGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACcagc                     <  2:157734/151‑1 (MQ=255)
                                                                                                                          aTCAGCGTATCGATTACGCGAGGGGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACAAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACcagc                     >  2:670123/1‑150 (MQ=255)
                                                                                                                                    cGATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTg          <  2:342667/151‑1 (MQ=255)
                                                                                                                                     gATTACGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGc         >  1:209174/1‑151 (MQ=255)
                                                                                                                                         aCGCGAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCt      <  2:472574/150‑1 (MQ=255)
                                                                                                                                          cgcgAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTc    <  1:882381/151‑1 (MQ=255)
                                                                                                                                           gcgAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTcc   <  2:234238/151‑1 (MQ=255)
                                                                                                                                            cgAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCa  <  1:906166/151‑1 (MQ=255)
                                                                                                                                             gAGGTGACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCACCTGCCCTTcc   >  2:662219/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCAGTATGGCATCGAGATACCGACAGTTGGCCTCTGTGGCGCGTTCGAGCATGCGCGGCAGATTGCGAAACTGCCAGTCAGGCCAGATGTAGCTCACTGCCGCCCACGCAATGGCACAACCAATCAGCGTATCGATTACGCGAGGTAACGCTACTTCAAAACCTTCACCCAGTAAGTTAAAACACAGTAGCACCAAAAGTGTGATGAACATCGTTGCATGAGCGTATTGCACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCA  >  NC_000913/1022209‑1022499

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