Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,930,568 +G intergenic (+66/‑373) kup → / → rbsD low affinity potassium transport system/RbsD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,930,5671.G98.2% 304.5 / ‑0.7 114intergenic (+65/‑374)kup/rbsDlow affinity potassium transport system/RbsD
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base G (46/66);  total (48/66)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 120 total aligned reads displayed.

TATGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATG  >  NC_000913/3930492‑3930637
                                                                              |                                                                      
tatGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGtt                                                                         <  1:1512529/76‑1 (MQ=255)
 atGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTg                                                                       <  1:3243817/76‑1 (MQ=255)
  tGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGGGCAAGGTTAGGGTTgg                                                                      >  2:3316025/1‑75 (MQ=255)
   gTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCCCTTTCTGGCGCGTGGTGGCAAGTTATGGTTggc                                                                     >  2:412867/1‑74 (MQ=255)
        tAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTNGCAAGTTATGGTTGGCGAGaa                                                                 <  1:2345067/76‑1 (MQ=255)
        tAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGaa                                                                 >  1:2535923/1‑76 (MQ=255)
        tAACGTCCGTCGGGCACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGGGGCAAGCTCTGGTTGGCGAGaa                                                                 >  1:3183119/1‑76 (MQ=255)
         aaCGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGaaa                                                                <  2:2986870/76‑1 (MQ=255)
           cgtccgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACg                                                              <  1:2847975/76‑1 (MQ=255)
           cgtccgtcGGGTACAGAGTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACg                                                              <  1:3532547/76‑1 (MQ=255)
              ccgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCa                                                           <  2:4695341/76‑1 (MQ=255)
               cgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaa                                                          >  2:637156/1‑76 (MQ=255)
               cgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaa                                                          <  1:571915/76‑1 (MQ=255)
               cgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaa                                                          <  2:2327679/76‑1 (MQ=255)
               cgtcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaa                                                          >  2:2809998/1‑76 (MQ=255)
               cgtcGGGTACAGATTGAAAAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaa                                                          <  1:1512337/76‑1 (MQ=255)
                 gcGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAAc                                                        <  2:812109/75‑1 (MQ=255)
                  cGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACg                                                       <  1:1211521/76‑1 (MQ=255)
                  cGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACg                                                       >  1:1644616/1‑76 (MQ=255)
                   gggTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGt                                                      <  1:2184757/76‑1 (MQ=255)
                   gggTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGt                                                      <  2:2897149/76‑1 (MQ=255)
                    ggTACAGATTGAACAACTGTCGCCGTCTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTa                                                     <  2:2677371/76‑1 (MQ=255)
                    ggTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTa                                                     <  2:3869442/76‑1 (MQ=255)
                     gTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTag                                                    >  1:2821828/1‑76 (MQ=255)
                      tACAGATTGAACGACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTaga                                                   <  1:4551122/76‑1 (MQ=255)
                      tACAGATTGAACAACTGTCGACCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAGACGCCAAACGCaga                                                   >  1:4402820/1‑76 (MQ=255)
                       aCAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGACACGCCAAACGTagaa                                                  >  2:1851487/1‑76 (MQ=255)
                       aCAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTagaa                                                  >  1:3166723/1‑76 (MQ=255)
                        cAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTagaag                                                 >  2:677402/1‑76 (MQ=255)
                         aGATTGAGCAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTagaaga                                                <  2:1397549/76‑1 (MQ=255)
                          gATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAgaagaa                                               <  2:4552049/76‑1 (MQ=255)
                          gATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAgaagaa                                               >  1:3154303/1‑76 (MQ=255)
                          gATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAgaagaa                                               >  1:2867383/1‑76 (MQ=255)
                           aTTGGACAACTGTGGCAGACGTGCTGGCGCGTGGTGGCAAGTGAGGGTGGGAGAGAAATGCGAAACGTcgaaggag                                              >  2:3083387/1‑73 (MQ=255)
                           aTTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGaagaag                                              <  2:4167508/76‑1 (MQ=255)
                            ttGAACAACTGTCGCTCTCTTTTTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGCAGAAGt                                             <  1:2147299/76‑1 (MQ=255)
                            ttGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGt                                             >  2:4196503/1‑76 (MQ=255)
                            ttGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGt                                             <  1:2972264/76‑1 (MQ=255)
                            ttGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGt                                             >  2:2517725/1‑76 (MQ=255)
                            ttGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGt                                             >  1:3866881/1‑76 (MQ=255)
                             tGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGtt                                            >  2:2039764/1‑76 (MQ=255)
                                acaccTGTCGCACACATTCTTGCGCGTTGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTc                                         <  1:3785874/76‑1 (MQ=255)
                                acaacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTc                                         <  1:3746559/76‑1 (MQ=255)
                                acaacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTc                                         <  2:2014130/76‑1 (MQ=255)
                                  aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaaa                                                         >  1:349184/1‑58 (MQ=255)
                                  aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCaaa                                                         <  2:349184/58‑1 (MQ=255)
                                  aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAGGAAGTTTTCCa                                       >  1:406113/1‑76 (MQ=255)
                                  aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCa                                       <  2:4407751/76‑1 (MQ=255)
                                  aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCa                                       <  2:707418/76‑1 (MQ=255)
                                  aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCa                                       <  2:1650375/76‑1 (MQ=255)
                                  aacTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCa                                       >  1:1900900/1‑76 (MQ=255)
                                   acTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCAc                                      <  2:4284187/76‑1 (MQ=255)
                                     tGTCGCGCACTTGCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCg                                    <  2:434828/76‑1 (MQ=255)
                                     tGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCg                                    >  2:2838583/1‑76 (MQ=255)
                                      gTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGc                                   <  1:2111033/76‑1 (MQ=255)
                                      gTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGc                                   >  1:1240110/1‑76 (MQ=255)
                                       tCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCt                                  <  1:2039765/76‑1 (MQ=255)
                                       ccGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCt                                  <  2:2120149/75‑1 (MQ=255)
                                        cGCCCAGTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTg                                 <  1:4328651/76‑1 (MQ=255)
                                         gCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTAGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGc                                >  1:2809202/1‑76 (MQ=255)
                                         gCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGc                                <  1:4287667/76‑1 (MQ=255)
                                         gCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGc                                <  2:2382629/76‑1 (MQ=255)
                                          cccGCTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCg                               <  2:461689/76‑1 (MQ=255)
                                          cccACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCg                               <  2:4465552/76‑1 (MQ=255)
                                          cccACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCg                               >  1:3284369/1‑76 (MQ=255)
                                           ccAGTTTCGGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCgg                              <  2:1018396/76‑1 (MQ=255)
                                             aCTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTc                            >  2:2017606/1‑76 (MQ=255)
                                              cTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCt                           >  1:3339877/1‑76 (MQ=255)
                                              cTTTCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCt                           >  1:1734165/1‑76 (MQ=255)
                                               tttCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTg                          <  1:1398775/76‑1 (MQ=255)
                                               tttCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTg                          >  2:1379107/1‑76 (MQ=255)
                                                ttCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTgg                         >  2:1530486/1‑76 (MQ=255)
                                                 tCTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGa                        >  1:654067/1‑76 (MQ=255)
                                                  cTGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGaa                       >  1:3580419/1‑76 (MQ=255)
                                                   tGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAg                      <  2:3496414/76‑1 (MQ=255)
                                                   tGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAg                      >  1:1345004/1‑76 (MQ=255)
                                                   tGGCGCGTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCGGCGGTCTGGCAg                      >  2:1658052/1‑76 (MQ=255)
                                                    ggCGCGGGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAgg                     <  2:33605/76‑1 (MQ=255)
                                                     gcgcGTGGTGGCAAGTTATGGTTGGCGCGAAACTCCACGCGTAGCATAGGTTTTCCGCCGCTGCGGTCTGGAATGa                    <  1:122372/76‑1 (MQ=255)
                                                      cgcgTGGTGGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGAt                   >  1:3237132/1‑76 (MQ=255)
                                                         gtggtgGCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTaa                >  2:4021995/1‑76 (MQ=255)
                                                         gtggtgGCAAGTTATGGCTGGCGAGACCCGCCACACGTAGAAGAAGTTTTCCACCGCTGCGGTGTGGCAGGATTac                >  2:4030724/1‑75 (MQ=255)
                                                            gtggCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGtt             <  2:670477/76‑1 (MQ=255)
                                                             tggCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTg            <  1:1165428/76‑1 (MQ=255)
                                                              ggCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGc           <  1:103705/76‑1 (MQ=255)
                                                               gCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGcc          <  1:1890864/76‑1 (MQ=255)
                                                               gCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGcc          <  2:590608/76‑1 (MQ=255)
                                                               gCAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGcc          >  1:1374616/1‑76 (MQ=255)
                                                                cAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCg         <  2:2637031/76‑1 (MQ=255)
                                                                cAAGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCg         <  1:817377/76‑1 (MQ=255)
                                                                 aaGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGt              <  1:3843748/70‑1 (MQ=255)
                                                                 aaGTTATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGt              >  2:3843748/1‑70 (MQ=255)
                                                                    gTTATGGTTGGCGAGCAACGCGAAAGGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGgat      <  2:347713/76‑1 (MQ=255)
                                                                      tATGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGgatga    >  1:1627463/1‑76 (MQ=255)
                                                                       aTGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGgatgat   <  1:4196503/76‑1 (MQ=255)
                                                                       aTGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGgatgat   <  1:504274/76‑1 (MQ=255)
                                                                       aTGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGgatgat   <  1:903371/76‑1 (MQ=255)
                                                                       aTGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTCGAAGGATTAAGTTGCCGgatgat   >  1:2627950/1‑76 (MQ=255)
                                                                        tGGTTGGCGAGAAACGCCAAACGTAGATGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTACGTTGCCGgatgatt  <  1:334428/76‑2 (MQ=255)
                                                                        tGGTTGGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGGatgatg  <  1:3583887/76‑1 (MQ=255)
                                                                              |                                                                      
TATGTCCATAACGTCCGTCGGGTACAGATTGAACAACTGTCGCCCACTTTCTGGCGCGTGGTGGCAAGTTATGGTTGCGAGAAACGCCAAACGTAGAAGAAGTTTTCCACCGCTGCGGTCTGGAAGGATTAAGTTGCCGGATGATG  >  NC_000913/3930492‑3930637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: