Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,542,564 C→T H720Y (CAT→TAT)  feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,542,5640CT100.0% 221.2 / NA 56H720Y (CAT→TAT) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/29);  total (27/29)

TCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTTGCTGCGCAGCCAGCACCA  >  NC_000913/3542418‑3542711
                                                                                                                                                  |                                                                                                                                                   
tCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATcc                                                                                                                                                 >  2:739549/1‑151 (MQ=255)
 ccTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATcc                                                                                                                                                 >  1:69027/1‑150 (MQ=255)
   tgctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAAc                                                                                                                                              >  1:168526/1‑151 (MQ=255)
     ctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACtt                                                                                                                                            >  2:757214/1‑151 (MQ=255)
       gtatgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACtt                                                                                                                                            <  1:448529/149‑1 (MQ=255)
         atgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACa                                                                                                                                         <  1:1100334/150‑1 (MQ=255)
         atgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACa                                                                                                                                         <  1:705804/150‑1 (MQ=255)
         atgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAg                                                                                                                                        <  2:183052/151‑1 (MQ=255)
             aCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTg                                                                                                                                    >  2:588348/1‑151 (MQ=255)
             aCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTg                                                                                                                                    <  1:398813/151‑1 (MQ=255)
              ccATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTgg                                                                                                                                   <  2:833514/151‑1 (MQ=255)
                aTGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGgtg                                                                                                                                 <  1:12704/151‑1 (MQ=255)
                  gTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGtgtg                                                                                                                               <  2:1156740/151‑1 (MQ=255)
                      ctcGGTGATGGGGGCAATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCAtt                                                                                                                           >  2:1003524/1‑151 (MQ=255)
                        cGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCAtt                                                                                                                           <  1:413229/149‑1 (MQ=255)
                        cGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCt                                                                                                                         <  1:1057047/151‑1 (MQ=255)
                          gTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTgg                                                                                                                       >  1:814979/1‑151 (MQ=255)
                              tGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGt                                                                                                                   >  1:984545/1‑151 (MQ=255)
                                ggggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGtt                                                                                                                  >  1:1070439/1‑150 (MQ=255)
                                 gggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTAt                                                                                                                >  1:803707/1‑151 (MQ=255)
                                     tATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTg                                                                                                            >  2:345907/1‑151 (MQ=255)
                                         gCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGttt                                                                                                         >  1:201780/1‑150 (MQ=255)
                                         gCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTa                                                                                                        <  1:679758/151‑1 (MQ=255)
                                          cccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTa                                                                                                        >  2:984118/1‑150 (MQ=255)
                                           ccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAAc                                                                                                      >  1:988101/1‑151 (MQ=255)
                                                   cAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTgg                                                                                               <  1:356267/150‑1 (MQ=255)
                                                   cAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTgg                                                                                               <  1:35349/150‑1 (MQ=255)
                                                    aaGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGtt                                                                                             >  1:486840/1‑151 (MQ=255)
                                                    aaGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGtt                                                                                             >  2:733389/1‑151 (MQ=255)
                                                     aGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTa                                                                                            <  1:466635/151‑1 (MQ=255)
                                                     aGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTa                                                                                            <  1:70757/151‑1 (MQ=255)
                                                         gTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCgg                                                                                        <  2:1175753/151‑1 (MQ=255)
                                                               ggATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTctgct                                                                                  <  1:831418/151‑1 (MQ=255)
                                                               ggATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGCTATCGGTctgc                                                                                   >  1:486862/1‑150 (MQ=255)
                                                                gATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTctgctg                                                                                 <  2:494089/151‑1 (MQ=255)
                                                                   gggCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTgcgc                                                                              <  2:1187224/151‑1 (MQ=255)
                                                                     gCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGcgccg                                                                            <  2:300969/151‑1 (MQ=255)
                                                                           ccATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCcgcgcgcg                                                                      >  2:1169437/1‑151 (MQ=255)
                                                                           ccATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCcgcgcgc                                                                       <  1:7804/150‑1 (MQ=255)
                                                                               ccTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGc                                                                  <  1:891837/151‑1 (MQ=255)
                                                                                      gggCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTgg                                                           >  1:685628/1‑151 (MQ=255)
                                                                                        gCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGat                                                         <  2:693497/151‑1 (MQ=255)
                                                                                                  gCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAActgc                                               >  1:1090747/1‑151 (MQ=255)
                                                                                                         cACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCaa                                         <  2:1021678/150‑1 (MQ=255)
                                                                                                             ggCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCg                                     <  1:1043339/150‑1 (MQ=255)
                                                                                                                   aTTGTTCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACGTCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCg                              >  2:908299/1‑151 (MQ=255)
                                                                                                                       ttCTATCAGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTa                           >  1:124528/1‑150 (MQ=255)
                                                                                                                              aGGTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTtg                   <  1:299453/151‑1 (MQ=255)
                                                                                                                                gTCGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTtgct                 >  1:1064279/1‑151 (MQ=255)
                                                                                                                                  cGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTtgctgc               >  2:1121347/1‑151 (MQ=255)
                                                                                                                                  cGCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTtgctgc               >  2:314489/1‑151 (MQ=255)
                                                                                                                                   gCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCACGTCGGTAAGCAGCtgc                  >  1:558433/1‑147 (MQ=255)
                                                                                                                                   gCCAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTTGCTgcg              <  2:1108987/151‑1 (MQ=255)
                                                                                                                                     cAGCTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTTGCTGCGca            <  2:885703/151‑1 (MQ=255)
                                                                                                                                        cTACAGTCAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTTGCTGCGcagcc         >  1:148168/1‑151 (MQ=255)
                                                                                                                                           cagtcagtATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTTGCTGCGcagccagc      >  2:870621/1‑151 (MQ=255)
                                                                                                                                               cAGTATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTTGCTGCGCAGCCAGcacca  <  1:972691/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTTATCCTGTTTAACATCGTGGTTATCGGTCTGCTGCGCCGCGCGCGTAGCCGGGTGGATATCGAACTGCTGGCAACCCGCAAGTCGGTAAGCAGTTGCTGCGCAGCCAGCACCA  >  NC_000913/3542418‑3542711

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: