Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,619,957 G→A *128* (TAG→TAA marA → multiple antibiotic resistance transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,619,9570GA100.0% 239.8 / NA 73*128* (TAG→TAAmarAmultiple antibiotic resistance transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/34);  total (39/34)

ATGGCTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATG  >  NC_000913/1619824‑1620091
                                                                                                                                     |                                                                                                                                      
atGGCTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGaa                                                                                                                                   <  2:600359/139‑1 (MQ=255)
          gTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         <  1:948760/139‑1 (MQ=255)
             gCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGt                                                                                                                      >  2:1613690/1‑139 (MQ=255)
                 caaACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACt                                                                                                                  <  2:1455013/139‑1 (MQ=255)
                         gACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAAt                                                                                                          >  1:1477714/1‑139 (MQ=255)
                         gACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAAt                                                                                                          >  1:365234/1‑139 (MQ=255)
                               aaCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCAt                                                                                                       <  2:1328494/136‑1 (MQ=255)
                               aaCCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCAt                                                                                                       >  1:1328494/1‑136 (MQ=255)
                                           ttACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATc                                                                                        <  2:1477714/139‑1 (MQ=255)
                                             aCTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCg                                                                                      >  1:94193/1‑139 (MQ=255)
                                                  gATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATa                                                                                 >  2:1299347/1‑139 (MQ=255)
                                                   aTGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATag                                                                                >  2:364862/1‑139 (MQ=255)
                                                      ttCCGCCGCATAAATACCCGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATcgcag                                                                             >  2:1184568/1‑139 (MQ=255)
                                                       tCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                            <  1:1106713/139‑1 (MQ=255)
                                                       tCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                            <  2:365234/139‑1 (MQ=255)
                                                         cgccgcATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTg                                                                          <  1:364862/139‑1 (MQ=255)
                                                             gcATAAATACCGGATGACCAATATGCAGGGCGAACCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGt                                                                                                                              >  1:646358/1‑83 (MQ=255)
                                                             gcATAAATACCGGATGACCAATATGCAGGGCGAACCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGt                                                                                                                              <  2:646358/83‑1 (MQ=255)
                                                              cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                            <  1:1548011/132‑1 (MQ=255)
                                                              cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                            >  2:1548011/1‑132 (MQ=255)
                                                              cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  2:928457/139‑1 (MQ=255)
                                                              cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     >  1:753576/1‑139 (MQ=255)
                                                              cATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  2:602273/139‑1 (MQ=255)
                                                                   aTACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATtct                                                                <  1:1299347/139‑1 (MQ=255)
                                                                        ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         >  2:944215/1‑77 (MQ=255)
                                                                        ggATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         <  1:944215/77‑1 (MQ=255)
                                                                         gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTc                                                          >  2:340444/1‑139 (MQ=255)
                                                                         gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTc                                                          >  2:1659301/1‑139 (MQ=255)
                                                                         gATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTc                                                          >  2:481338/1‑139 (MQ=255)
                                                                          aTGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTcc                                                         >  2:401184/1‑139 (MQ=255)
                                                                           tGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCcg                                                        >  1:1676987/1‑139 (MQ=255)
                                                                                      cAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGtt                                             >  1:75161/1‑139 (MQ=255)
                                                                                       aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         <  2:465896/62‑1 (MQ=255)
                                                                                       aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         >  1:465896/1‑62 (MQ=255)
                                                                                       aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                            >  1:649124/1‑107 (MQ=255)
                                                                                       aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATagcagc                                                                            <  2:649124/107‑1 (MQ=255)
                                                                                       aGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTg                                            >  1:4366/1‑139 (MQ=255)
                                                                                         ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            >  1:520877/1‑73 (MQ=255)
                                                                                         ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            <  2:520877/73‑1 (MQ=255)
                                                                                         ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                   <  2:1208454/82‑1 (MQ=255)
                                                                                         ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                   >  1:1208454/1‑82 (MQ=255)
                                                                                         ggCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCg                                          >  2:1517358/1‑139 (MQ=255)
                                                                                          gCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            >  1:1348725/1‑72 (MQ=255)
                                                                                          gCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCa                                                                                                            <  2:1348725/72‑1 (MQ=255)
                                                                                           cGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGa                                        >  1:945365/1‑139 (MQ=255)
                                                                                              aTCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACa                                     >  1:1618491/1‑139 (MQ=255)
                                                                                                    cTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         >  2:1341367/1‑49 (MQ=255)
                                                                                                    cTTTTTACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACaa                                                                                                                         <  1:1341367/49‑1 (MQ=255)
                                                                                                        ttACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAg                           <  2:945365/139‑1 (MQ=255)
                                                                                                         tACATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGc                          <  2:753576/139‑1 (MQ=255)
                                                                                                          aCATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGcc                         <  2:94193/139‑1 (MQ=255)
                                                                                                           cATCCATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCa                        <  1:1517358/139‑1 (MQ=255)
                                                                                                              ccATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCGGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgtt                     <  1:54360/139‑1 (MQ=255)
                                                                                                               cATTAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttg                    >  2:52622/1‑139 (MQ=255)
                                                                                                                 ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACt                                                                                              <  2:1395071/63‑1 (MQ=255)
                                                                                                                 ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACt                                                                                              >  1:1395071/1‑63 (MQ=255)
                                                                                                                 ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                  >  2:770087/1‑139 (MQ=255)
                                                                                                                 ttAAATCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAgttgtt                  >  1:1535244/1‑139 (MQ=255)
                                                                                                                     aTCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGtt                                             >  1:1011205/1‑108 (MQ=255)
                                                                                                                     aTCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGtt                                             <  2:1011205/108‑1 (MQ=255)
                                                                                                                      tCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTActtc             >  2:1760825/1‑139 (MQ=255)
                                                                                                                      tCATTACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTActtc             <  2:1618491/139‑1 (MQ=255)
                                                                                                                         ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                   <  1:123358/50‑1 (MQ=255)
                                                                                                                         ttACAACAGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGaaa                                                                                                   >  2:123358/1‑50 (MQ=255)
                                                                                                                             aacaGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGcc      <  1:401184/139‑1 (MQ=255)
                                                                                                                             aacaGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGcc      <  2:4366/139‑1 (MQ=255)
                                                                                                                             aacaGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGcc      <  2:1535244/139‑1 (MQ=255)
                                                                                                                               caGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     >  1:1319361/1‑74 (MQ=255)
                                                                                                                               caGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCtt                                                                     <  2:1319361/74‑1 (MQ=255)
                                                                                                                                aGCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAAt   >  2:38336/1‑139 (MQ=255)
                                                                                                                                 gCTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATg  >  1:474673/1‑139 (MQ=255)
                                                                                                                                  cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGcc      <  1:1686881/134‑1 (MQ=255)
                                                                                                                                  cTAATTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGcc      >  2:1686881/1‑134 (MQ=255)
                                                                                                                                     |                                                                                                                                      
ATGGCTTCGAGTCGCAACAAACTCTGACCCGAACCTTCAAAAATTACTTTGATGTTCCGCCGCATAAATACCGGATGACCAATATGCAGGGCGAATCGCGCTTTTTACATCCATTAAATCATTACAACAGCTAGTTGAAAACGTGACAACGTCACTGAGGCAATCATGAAACCACTTTCATCCGCAATAGCAGCTGCGCTTATTCTCTTTTCCGCGCAGGGCGTTGCGGAACAAACCACGCAGCCAGTTGTTACTTCTTGTGCCAATG  >  NC_000913/1619824‑1620091

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: