Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 516,951 G→A V123V (GTG→GTA fetB → iron export ABC transporter permease; peroxide resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913516,9510GA100.0% 191.1 / NA 58V123V (GTG→GTAfetBiron export ABC transporter permease; peroxide resistance protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/32);  total (26/32)

AATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAG  >  NC_000913/516808‑517093
                                                                                                                                               |                                                                                                                                              
aaTGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCt                                                                                                                                          <  2:879045/150‑1 (MQ=255)
        gTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCggg                                                                                                                                 <  1:372227/151‑1 (MQ=255)
          ggAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGgat                                                                                                                               <  1:569834/151‑1 (MQ=255)
            aaCGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTTTCTCCGGGGCGGTTGAATTTATCCCGATGCAGGGAATCCCTATCGCCggg                                                                                                                                 >  2:688501/1‑147 (MQ=255)
                   aaaaaCGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCgg                                                                                                                      <  2:577696/151‑1 (MQ=255)
                    aaaaCGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGt                                                                                                                     <  1:583909/151‑1 (MQ=255)
                             tAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGGTTATTTTCTTCGGGGCGGTTGTATTTTTTCCGATGCAGGGAATCCCTCTCGCCGGGATGGTTGCCGGTAACGccct                                                                                                             >  2:251729/1‑148 (MQ=255)
                              aaaTATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGGTGATTCTCTCAGGGGCGATTGAATTTATTCCGGTGCCGGGAATCCCTATCGCCGGGGTGATTGCCGGGAAAGCCATg                                                                                                            >  2:508785/1‑150 (MQ=255)
                               aaTATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTCGCGGGGCTTATTTTCTCCGGGTCGCTTGGATTTTTTCCGATGCCGGGAATCCCCATCCCCGGGGTGATTTGCGGTAACGcc                                                                                                               >  2:719616/1‑146 (MQ=255)
                                 tatatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTTTCTCAGGGTAGGTTGAAATTTTCCCGATGGAGGGAATCCCTATCGCCGGGGTGATTGCCGGGAACGGCATgg                                                                                                           >  2:313046/1‑148 (MQ=255)
                                   tatTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTTCCCTGGCGGTGCTGATTTTCTCTGGGTCGATTTGATTTATCCCGATGCAGGTAATTCCTATCGCCCGGATGGTTGCCGGTAACGCCATGGTAGCg                                                                                                      >  2:899550/1‑151 (MQ=255)
                                      tGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCGCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGt                                                                                                    <  2:18090/150‑1 (MQ=255)
                                        cTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGGGCTGATTTTCTCCGGGGCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGGTGATTGCCGGGAACGCCATGGGAGCGGTggg                                                                                                 >  2:642798/1‑151 (MQ=255)
                                        cTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGGGCTGATTTTCTCAGGGGCGATTGAATTTATTCCGATGCAGGGAATCCCTATTGCCGGGATGATTTCCGGTAACGCCATGGTAGCgggggg                                                                                                 >  2:601654/1‑151 (MQ=255)
                                          aaaGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGt                                                                                               <  2:53450/151‑1 (MQ=255)
                                          aaaGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGGCGATTGAATTTATCCCGATGCAGGGAATCCCTATCGCCGGGGTGATTGCCGGGAACGCCATGGTAGCGgggggg                                                                                                >  1:347571/1‑150 (MQ=255)
                                          aaaGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGGGCTTATTTTCTCAGGGGCGATTGAATTTATCCCGCTGCAGGGAATCCCTATCGCCGGGGTGATTTCCGGTAACCCCCTGGGAGCGGTgggg                                                                                                >  2:185245/1‑150 (MQ=255)
                                                 ttATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTTATTTTCTCAGGGTCGGTTGGATTTATCCCGATGCAGGTAATCCCTTTCGCCGGGGTGATTGCCGGTAACGCCATGGGAGCGGTGgggttgggtt                                                                                         >  2:309868/1‑150 (MQ=255)
                                                 ttATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGtt                                                                                         <  1:864059/150‑1 (MQ=255)
                                                 ttATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGGCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGGAGCGGGGGGGGTGTGtt                                                                                         >  1:536139/1‑150 (MQ=255)
                                                 ttATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGGGCTGATTCTCTCAGGGGCGATTGAATTTATCCCGATGCAGGGAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGGGGGGTTGTGtt                                                                                         >  1:381911/1‑150 (MQ=255)
                                                     ctcATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTTTCTCAGGGTCGATTGAATTTATTCCGATGGAGGGAATCCCTATCGCCGGGATGATTTCCGGTAACGCCATGGGAGCGGGGGGGGTGTGGTacaa                                                                                     >  2:647002/1‑150 (MQ=255)
                                                      tcATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATTCCGCTGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGGAGCGGGGGGGTTGTGTTacaa                                                                                     >  2:230815/1‑149 (MQ=255)
                                                        aTCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCCGGCGGGGCTGATTTTTTCAGGGGCGATTGAATTTTTTCCGATGCAGGGAATCCCTATTGCCGGGGTGATTTGCGGGAAAGCCCTGGGAGCGGGGGGGGTGGGTTTtaacaa                                                                                  >  2:477550/1‑150 (MQ=255)
                                                             ttATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTTATTTTCTCAGGGGCGATTGAATTTATTCCGATGCAGGGAATCCCTATCGCCGGGATGGTTGCCGGTAACGCCCTTGTAGCGGTGGGGTTGTGTTACAACAATTTAgg                                                                            >  2:784285/1‑151 (MQ=255)
                                                                  gCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCaa                                                                        <  1:598947/150‑1 (MQ=255)
                                                                    cATTACGGTCGGGGCGGGAATTACCCTGGCGGGGCTGATTCTCTCAGGGGTGATTGAATTTATCCCGATGCAGGGAATCCCTTTCGCCGGGGTGTTTGGCGGTAAAGCCATGGTAGCGGGGGGGGTGTGTTTCAAAAAATTTGGGGAAAgg                                                                     >  2:796200/1‑151 (MQ=255)
                                                                    cATTACGGTCGGGGCGGGAATTAACCTGGCGGTGTTGGTTCTCTTAGGGGTGGGTGGATTTATATCGATGCAGGGAATCCCTTTAGCCGGGGAGAGTGGCAGGAAAACCATGGGAGCGGTGGGGGTGTGTTTTAAAAATTTAAGGCAACgg                                                                     >  2:628450/1‑151 (MQ=255)
                                                                       tACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGt                                                                   <  1:688550/150‑1 (MQ=255)
                                                                        aCGGTCGGGGCGGGAATTACCCTGGCGGGGCTGATTCTCTCCGGGGCGGTTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGGTGATTTCCGGTAACGCCATGGTAGCGGTGGGGGTGTGTTTCAACAATTTAGGGCAACggggc                                                                  >  2:346376/1‑148 (MQ=255)
                                                                               gggCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGa                                                          <  2:462138/151‑1 (MQ=255)
                                                                                 gCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAAc                                                        <  2:694643/151‑1 (MQ=255)
                                                                                  cGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACa                                                       <  1:185270/151‑1 (MQ=255)
                                                                                     gAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCa                                                    <  2:347545/151‑1 (MQ=255)
                                                                                       aTTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAAc                                                  <  2:181635/151‑1 (MQ=255)
                                                                                       aTTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAAc                                                  <  1:855873/151‑1 (MQ=255)
                                                                                        ttACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAAc                                                  <  1:860370/150‑1 (MQ=255)
                                                                                         tACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACa                                                 <  2:376633/150‑1 (MQ=255)
                                                                                            ccTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCATATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGAtt                                             >  1:299597/1‑151 (MQ=255)
                                                                                               ggCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCa                                           <  2:536089/150‑1 (MQ=255)
                                                                                                  ggTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGTAACGGGTCATTAGCGAACAGCAACAGATTCAGgag                                       <  1:435143/151‑1 (MQ=255)
                                                                                                  ggTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAgg                                         >  2:130863/1‑149 (MQ=255)
                                                                                                     gCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGaaa                                    <  2:111839/151‑1 (MQ=255)
                                                                                                     gCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGaa                                     <  1:97200/150‑1 (MQ=255)
                                                                                                        gATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACt                                  <  1:434600/150‑1 (MQ=255)
                                                                                                                cAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTgg                         <  1:3130/151‑1 (MQ=255)
                                                                                                                cAGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTgg                         <  1:893439/151‑1 (MQ=255)
                                                                                                                 aGGGTCGATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGt                        <  2:276923/151‑1 (MQ=255)
                                                                                                                       gATTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACg                  <  2:371532/151‑1 (MQ=255)
                                                                                                                        aTTGAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGc                 >  1:649482/1‑151 (MQ=255)
                                                                                                                           gAATTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGa              <  1:754594/151‑1 (MQ=255)
                                                                                                                            aaTTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAGCTGAGTCTTGGTGCGACGCCGa              <  2:792876/150‑1 (MQ=255)
                                                                                                                            aaTTTATCCCGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGa              >  2:169868/1‑150 (MQ=255)
                                                                                                                                    ccGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCtt     >  2:307701/1‑151 (MQ=255)
                                                                                                                                    ccGATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCtt     <  2:446017/151‑1 (MQ=255)
                                                                                                                                      gATGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTc    >  1:689503/1‑150 (MQ=255)
                                                                                                                                      gATGCAGGTAATCCATATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCa   <  2:299571/151‑1 (MQ=255)
                                                                                                                                       aTGCAGGTAATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAg  <  2:156158/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
AATGCGGCGTGGAACGCGCAAAAACGCAGTAAATATATTGCTAAAGCTTTTATCTCATCGTTTATTGCCATTACGGTCGGGGCGGGAATTACCCTGGCGGTGCTGATTCTCTCAGGGTCGATTGAATTTATCCCGATGCAGGTGATCCCTATCGCCGGGATGATTGCCGGTAACGCCATGGTAGCGGTGGGGTTGTGTTACAACAATTTAGGGCAACGGGTCATTAGCGAACAGCAACAGATTCAGGAGAAACTGAGTCTTGGTGCGACGCCGAAGCAGGCTTCAG  >  NC_000913/516808‑517093

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