Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,831,451 A→G E332G (GAA→GGA)  yicH → putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,831,4510AG100.0% 208.3 / NA 57E332G (GAA→GGA) yicHputative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/39);  total (18/39)

AAGATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTA  >  NC_000913/3831307‑3831593
                                                                                                                                                |                                                                                                                                              
aaGATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAggggg                                                                                                                                          >  1:780420/1‑151 (MQ=255)
 aGATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAggggg                                                                                                                                          >  2:428584/1‑150 (MQ=255)
  gATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTa                                                                                                                                        <  1:408665/151‑1 (MQ=255)
  gATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTa                                                                                                                                        <  1:414370/151‑1 (MQ=255)
      aCTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGGTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATg                                                                                                                                      >  2:38232/1‑149 (MQ=255)
              acacAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGt                                                                                                                            <  1:243007/151‑1 (MQ=255)
                acaGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCa                                                                                                                          <  1:791808/151‑1 (MQ=255)
                acaGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCa                                                                                                                          <  2:157806/151‑1 (MQ=255)
                   ggAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGTTGGGGAGGGGGTATGGTCAGAACGTCAgg                                                                                                                        <  1:576620/150‑1 (MQ=255)
                    gAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGa                                                                                                                      <  1:112249/151‑1 (MQ=255)
                      aGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAAc                                                                                                                    <  1:869617/151‑1 (MQ=255)
                       ggCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACt                                                                                                                   <  2:206117/151‑1 (MQ=255)
                         cAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTgg                                                                                                                 <  1:744426/151‑1 (MQ=255)
                               gTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGt                                                                                                           <  1:953853/151‑1 (MQ=255)
                                    tGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACg                                                                                                       <  2:260386/150‑1 (MQ=255)
                                     gAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACgg                                                                                                      <  1:427860/150‑1 (MQ=255)
                                      aaCGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGa                                                                                                    <  2:332994/151‑1 (MQ=255)
                                         gCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGaaaa                                                                                                 <  2:30791/151‑1 (MQ=255)
                                          cTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAAc                                                                                                >  1:256727/1‑151 (MQ=255)
                                            aGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCTCCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACg                                                                                               <  1:815785/150‑1 (MQ=255)
                                              cGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGtt                                                                                             <  2:244816/150‑1 (MQ=255)
                                                    cATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCtt                                                                                      <  1:928890/151‑1 (MQ=255)
                                                     aTTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTtg                                                                                     <  2:936660/151‑1 (MQ=255)
                                                       ttATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTtga                                                                                    <  1:709379/150‑1 (MQ=255)
                                                              tCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCg                                                                             <  1:522578/150‑1 (MQ=255)
                                                               cGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCgg                                                                            <  2:814313/150‑1 (MQ=255)
                                                                  tGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAAt                                                                        >  2:596393/1‑151 (MQ=255)
                                                                   gCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAAtt                                                                       <  1:571219/151‑1 (MQ=255)
                                                                        tatttGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCgg                                                                  >  1:887304/1‑151 (MQ=255)
                                                                         atttGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCggg                                                                 <  1:87768/151‑1 (MQ=255)
                                                                              aCCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTgg                                                             >  2:341231/1‑150 (MQ=255)
                                                                                     tataAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATAcc                                                     >  1:297247/1‑151 (MQ=255)
                                                                                        aaaCACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCtt                                                   <  2:353189/150‑1 (MQ=255)
                                                                                           cacGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCg                                               >  2:734671/1‑151 (MQ=255)
                                                                                           cacGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCg                                               <  1:595894/151‑1 (MQ=255)
                                                                                           cacGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCg                                               <  2:799444/151‑1 (MQ=255)
                                                                                             cGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGaa                                             <  1:998675/151‑1 (MQ=255)
                                                                                               gAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGaaaa                                           <  1:535148/151‑1 (MQ=255)
                                                                                                     tCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGc                                     >  1:465292/1‑151 (MQ=255)
                                                                                                      ccccGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGc                                     <  1:576484/150‑1 (MQ=255)
                                                                                                         cGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACa                                 <  1:883290/151‑1 (MQ=255)
                                                                                                          gCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAAc                                  >  1:64042/1‑149 (MQ=255)
                                                                                                           cAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGt                               <  1:177638/151‑1 (MQ=255)
                                                                                                             gggCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTtg                             >  1:655267/1‑151 (MQ=255)
                                                                                                              ggCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTtgt                            <  2:389790/151‑1 (MQ=255)
                                                                                                               gCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTtgtg                           >  1:517044/1‑151 (MQ=255)
                                                                                                                cGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGtgg                          >  2:729358/1‑151 (MQ=255)
                                                                                                                cGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGtgg                          >  1:832442/1‑151 (MQ=255)
                                                                                                                 gTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGtgga                         <  2:832366/151‑1 (MQ=255)
                                                                                                                 gTAGCGCTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGtgg                          >  2:753529/1‑150 (MQ=255)
                                                                                                                     cgcTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGtggatgga                     <  1:821364/151‑1 (MQ=255)
                                                                                                                       cTGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGaa                    <  1:150972/150‑1 (MQ=255)
                                                                                                                        tGCGCCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAAc                  >  1:349294/1‑151 (MQ=255)
                                                                                                                         gcgcCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGaaa                   <  1:179545/149‑1 (MQ=255)
                                                                                                                          cgcCAGTTCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGa                <  2:780339/151‑1 (MQ=255)
                                                                                                                                ttCACCAGCCGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCg          >  1:306972/1‑151 (MQ=255)
                                                                                                                                        ccGCTGGGGAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTa  <  1:560926/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
AAGATGACTGGCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCGGTTGGTTA  >  NC_000913/3831307‑3831593

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